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Summary 

y 

The  objective  of  the  research  was  to  provide  theoretical  and 
experimental  backing  for  three  concepts  that  will  facilitate  the 
design  of  structures  capable  of  arresting  cracks.  The  first 
concept  is  that  the  conventional  static  toughness  measures,  JIc  or 
the  JR  curve,  are  bona  fide  measures  of  the  fracture  arrest 
capability  of  tough  ductile  steels.  The  second  is  that  the 
J-values  represent  the  crack  driving  force  when  a  crack  arrests  in 
a  tough  material,  provided  the  bulk  of  the  propagation  just  prior 
to  arrest  occurs  in  a  relatively  brittle  material. 

The  work  includes : 

(1)  Verification  of  J-integral  evaluation  of  arrest  properties.' 

(a)  Analysis  of  the  crack  velocity  dependence  of  toughness  for 
the  fibrous  mode  of  crack  extension. 

(b)  A  new  test  method  for  studying  fast  fracture  and  arrest  in 
tough  steels. 

(c)  Measurements  of  fast  fracture  and  crack  arrest. 

(d)  Fractographic  studies.’ 

(2)  Application  of  J-resistance  curve  concept  to  crack  arrest.' 

(a)  Development  of  a  finite  element  model  for  analysis  of  a 
rapidly  loaded  stationary  crack  of  a  compact  specimen  of 
a  tough,  rate  dependent  steel. 

(b)  Development  of  a  general  method  for  quantifying  temperature 
and  rate-sensitive  constitutive  relations  for  steel. 

(c)  Development  of  crack  extension  capability  using  "local 
control"  in  the  finite  element  program. 

(d)  Finite  element  modeling  of  rapid  crack  extension  in  ductile 
viscoplastic  materials. 


Proceedings  of  US-Japan  Cooperative  Seminar 
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DEFINITION  OF  CRACK  ARREST  PERFORMANCE  OF  TOUGH  ALLOYS 


V.  Dantatn  and  G.T.  Hahn 

Mechanical  and  Materials  Engineering, 

Vanderbilt  University,  Nashville,  TN 

1.  Introduction 

Tough  alloys  are  inserted  into  structures  to  arrest  unstable 
fractures  but  the  performance  of  such  materials  in  terms  of  the  size  of 
cracks  that  will  be  stopped  at  different  stress  levels  is  not  well 
defined.  The  conventional  wedge-loaded  compact  specimen  crack  arrest 
test  is  difficult,  to  perform  on  tough  materials  for  which  O.ZSm*/^ 

<  KIa/oc  <  0.5nr'^.  The  test  calls  for  large  stress  intensity  values  at 
the  onset  of  cracking,  i.e.  Kj  >  2Kxc,  and  prohibitive  specimen 
size  requirements,  '  i.e.,  W  >  8(KIa/o-)  >  lm.  For  this  reason, 

the  possibility  of  evaluating  Kja — the  J-equivalent  of  Kja  expressed  in 
stress  intensity  units — without  recourse  to  a  crack  arrest  test  deserves 
attention.  The  basis  for  this  idea  is:  (1)  the  fact  that  tough  alloys 
normally  display  flat  fibrous  fractures  with  shear  lips  or  the  (fibrous) 
full  shear  mode,  and  (2)  evidence  that  the  resistance  to  penetration  by 
a  fibrous  fracture  normally  increases  continuously  with  crack  velocity 
(see  Figure  1).  In  this  case,  KIa  (the  minimum  of  the  Kjo-veloclty 
curve)  Is  found  at  zero  velocity  (see  Figure  lb)  and  corresponds  with  KjC. 
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Figure  1.  Schematic  representation  of  the  dependence  of  the  fast  fracture 
toughness  (K^)  on  crack  velocity:  (a)  cleavage  fracture,  and 
(b)  fibrous  (ductile)  fracture. 


This  paper  briefly  summarizes  the  evidence  that  the  fibrous  fracture 
resistance  increases  with  crack  velocity.  It  offers  two  'working'' 
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definitions  for  Kja  for  tough  alloys  based  on  the  J  jc-  anc*  (  Aa  “l*50mm)_ 
values  obtained  with  conventional  slow  loading.  Finally,  it  describes  an 
experimental  approach  utilizing  deep  side  grooves  and  preliminary  results 
for  a  7075-T651  Aluminum  plate. 

2.  Previous  Work 

2. 1  Experimental  Measurements 

Evidence  that  the  resistance  to  fibrous  fracture  displays  a 
positive  crack  velocity  dependence  is  reviewed  in  Reference  2  and  is 
drawn  from  a  number  of  sources: 

(1)  CVN-energy  values  for  the  fibrous  mode  (on  the  ductile  shelf)  for 
ABS-C,  A302B,  HY-80,  A  517F,  HY-130  and  a  10NiCr-Mo-Co  steel  are  higher 
for  standard  impact  rates  of  loading  than  for  slow  bending. 

(ii)  Kjc~values  for  the  fibrous  mode  of  crack  extension  in  steels  and 
high  rates  of  loading  exceed  the  values  obtained  with  conventional  slow 
rates  of  loading.  The  only  exception  is  the  results  of  Klepaczko^*^ 
for  PA6  and  UTD  5U2A  aluminum  which  reveal  a  40%  decrease  in  Kjq  with 
increasing  loading  rate  in  the  range  10^  <  X  <  10^  MPam^^g-l. 

(iii)  Kja-values  for  relatively  high  strength  AIS1  4140  and  4340  steel 
on  the  ductile  shelf  fall  close  to  or  slightly  above  the  corresponding 
Kic-values. 

(iv)  Actual  measurements  of  Kjq,  the  resistance  to  fast  fracture,  for 
relatively  high  strength  steels  all  show  an  increase  in  Kjd  with  crack 
velocity  for  the  fibrous  mode. 

(v)  The  J-resistance  curve  of  A533  B  on  the  ductile  shelf  for  high 
rates  of  loading  that  produce  a  crack  velocity  of  V  •  0.25ms-1  is 
elevated  by  210KJm-^  (1200  lbin-1)  relative  to  conventional  J -measurements 
(V~10~5  ms-1). 

2.2  Theoretical  Calculations 

Freund  and  Douglas^***)  have  estimated  the  crack  velocity  dependence 
of  the  ductile  (fibrous)  toughness  on  the  basis  of  dynamic,  elastic  plastic, 
Mode  III  (antiplane  shear  analysis)  which  predicts  a  plastic  zone 
contraction  and  steeper  strain  gradients  directly  in  front  of  the  crack 
with  increasing  crack  velocities  above  V/Cs  >  0.2  (Cs  is  the  shear  wave 
velocity).  The  present  authors^)  have  incorporated  this  result  in  a 
simplified  model  of  the  Mode  I  plastic  strain  gradient  in  an  effort  to 
evaluate: 

(i)  The  plastic's  strain  rates  generated  by  a  fast  propagating  crack 
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and  the  resulting  elevation  of  the  yield  stress, 

(ii)  The  adiabatic  heating^)  and  corresponding  reduction  in  the  stress 
field  intensity,  and 

(iii)  The  increase  in  applied  stress  intensity  (fracture  toughness) 
required  to  offset  the  net  reduction  of  the  size  of  the  "process  zone". 

Figure  2  illustrates  that  the  strain  gradient  of  the  model  —  a  stationary 
crack  —  is  similar  to  the  near  tip  gradient  accompanying  steady  state 
crack  extension  proposed  by  Hermann  and  Rice.(®)  Figure  3  illustrates 
substantial  increases  of  Kjp  with  crack  velocity  calculated  for  fibrous 
mode  fracture  in  a  number  of  materials.  The  calculations  are  based  on 
assumed  critical  plastic  zone  strain  -  distance  values  which  are  independent 
of  crack  velocity.  However,  this  assumption  may  not  be  valid  when  the 
number  of  voids  nucleated  ahead  of  the  advancing  crack  tip  is  sensitive  to 
the  peak  normal  stress  and  increases  with  velocity  (or  loading  rate)  as  a 
result  of  the  strain  rate  induced  elevation  of  the  yield  stress.  This 
would  have  the  effect  of  reducing  the  critical  strain-distance  value  and 
the  toughness  at  higher  velocities  (and  loading  rates)  along  the  lines 
reported  by  Klepaczko. Consequently,  the  calculations  in  Figure 
3  need  the  support  of  the  f ractographic  observations  that  confirm  that  the 
dimple  population  does  not  increase  excessively  with  crack  velocity. 
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Figure  2.  Strain  gradient  in 
advance  of  crack  (p  is  extent 
of  plastic  zone). 

3.  General  Approach 


Figure  3.  Calculate  variation  of 
Kn)/Kic  with  crack  velocity. 


In  view  of  the  foregoing,  the  value  of  Kla  for  fibrous  fracture 
can  be  equated  with  Kjq  or  Kjq  (the  stress  intensity  equivalent  of  Jjc) 
or  more  correctly,  with  Kjr,  the  stress  intensity  equivalent  of  the  point 
on  the  J-resistance  curve  coresponding  with  crack  instability.  For  tough 


materials  the  J-res istance  curve  itself  comes  closest  to  a  geometry  indepen¬ 
dent  material  property  that  reflects  the  crack  penetration  resistance. 
However,  two  points  on  the  J-resistance  curve,  expressed  as  stress  intensity 
parameters,  can  serve  as  convenient  crack  arrest  design  parameters: 

KJaU)  =  Kj  -  [JIC  E/Cl-v2)  ] 1/2  O) 

KJaU)  =  Kjr«  IJr  E/(1-v2)]1/2  (2) 

where  Jr  corresponds  with  a  Aa  ■  1.50mm  offset  from  the  blunting 
line,  for  a  specimen  with  10%  to  20%  deep  side  grooves.  The  quantity 
Kja(l)  is  an  appropriate  measure  of  arrest  toughness  for  relatively 
short  cracks,  i.e. ,  a  <  100mm,  stopping  near  a  free  surface  where  large 
values  of  dJ/da  are  encountered.  The  quantity  Kja(2)  is  appropriate  for 
the  arrest  of  long  cracks,  i.e.,  a  >  500mm,  remote  from  a  free  surface, 
conditions  for  which  unstable  tearing*  are  not  likely  to  be  met  in  the 
first  Aa  ■  1.50mm  of  penetration.  The  Aa  ■  1.50mm  offset  is  selected, 
because  this  JR-value  is  evaluated  as  part  of  the  ASTM  (tentative) 
standard  Jic“determination  and,  because  the  reliability  of  JR-values 
for  larger  crack  extensions  is  currently  open  to  question.  Values  of 
Kja  (1)  and  Kja(2)  for  a  number  of  structural  alloys  based  on  measure¬ 
ments  by  the  authors,  Gudas  and  Co-workers^) ,  ancj  Joyce  and  Hanson (  »  ' 

and  estimates  of  the  crack  arrest  capabilities  are  summarized  in  Table  1. 

4  Experimental  Verification 

Work  is  in  progress  to  test  the  ideas  introduced  in  the  preceeding 
sections,  specifically  to: 

(i)  Devise  techniques  for  producing  fast,  unstable  fracture  in  tough 
arrester  grade  materials, 

(ii)  Measure  the  fibrous  mode  KjQ-velocity  curve, 

(ili)  Examine  fractographs  for  changes  in  the  fracture  appearance  with 
crack  velocity  and 

(iv)  Test  the  usefulness  of  Kja(l)  and  Kja(2)  as  structural  design 
parameters. 

As  a  first  step,  the  possibility  of  employing  deep  face  grooves  to 
reduce  the  stored  elastic  energy  requirement  for  fast  fracture  has 
been  explored.  Conventional  J-  curve  measurements  have  been  performed 
on  7075-Aluminum  in  different  heat  treated  conditions  and  on  HY-130 
steel,  with  45°  face  grooves  occupying  0%,  25%  and  80%  of  the  cross  section. 
The  unloading  compliance  was  used  to  evaluate  the  effective  crack  length, 
and  the  J-values  were  corrected  for  face-groove  depth:  J  ■  J(measured) 
(B/B^)**.  Some  of  the  results  of  this  study  are  summarized  in  Figure  4 

*  dJ/da  >  dJR/da 

**  B  is  the  ungrooved  plate  thickness  and  Br,  the  net  thickness  after 
grooving. 


TMLE  1.  MACK  ARREST  TOUGHNESS  VALUES  FOR  SELECTED  STEELS  FOR  SERVICE  TEMPERATURES  CORRESPONDING  TO  THE  UPPER  SHELF 
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CVN.J 
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CRACK  ARREST 
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REFERENCE 

HV-130 

RT 
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112,  116 
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0 

HY-130 

2S*C 

002 

80 

189 

— 

10 

HV-130 

RT 

082 
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2 

HT-00 

RT 

550 
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259 
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RT 
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— 

11 

A533-BI02) 

140 

413 

176 

207 

329 

1620 

0.  12 

A533-BI03) 

140 

401 

87 

202 

274 

1190 

0 

A51C-G70 

140 

303 

137 

180 

259 

1860 

0 

7075  A1 

RT 

280 

— 

78 

85 

235 

PrpSMt  Study 

707S-A1 

RT 

552 

— 

35 

56 

26 

Prtssnt  Study 

•Lsngth  of  •  contrail/  located  crock  2a  tbst  can  be  (topped  by  long,  noarinallj.tlatilc  arrottort,  «  •  O.S«B 
(o0  Is  the  /laid  (trass),  located  on  either  side  of  the  crack:  2a  •  */•£*)*  Vo]  . 


which  compares  the  Kjr/Kjr(25%  grooved)-ratios  obtained  at  different 
relative  toughness  levels.  The  results  show  that  side  grooves 
occupying  up  to  80%  of  the  cross  section  do  not  alter  either  the  JiC“  or 
Jg-values  substantially.  It  appears  that  the  use  of  deep  face  grooves  to 
reduce  stored  energy  requirements  for  a  fast  fracture  can  be  justified. 

The  effects  of  face  grooves  on  fibrous  mode  crack  arrest  measurements 
are  being  studied.  So  far,  run-arrest  events  have  been  successfully 
produced  in  25mm-thick,  150mm  by  150mm  compact  specimens  of  7075-T651 
with  blunt,  2mm-radius  starting  slots.  The  relatively  low  touahness-to- 
yield  ratio  of  this  heat  treated  condition:  Kj q/  <jq  m  0.07  m1^  permits 
an  LEFM  interpretation  of  the  arrest  event  and  the  evaluation  of  Kja. 
Results  are  summarized  in  Table  2  and  show  little  influence  of  side 
groove  depth.  In  addition,  the  results  seem  to  be  consistant  with  the 
proposed  features  of  fibrous  mode  crack  arrest: 

(i)  The  Kia-values  are  close  to  and  somewhat  larger  than  Kja(2), 

(ii)  The  Kjp-values  (the  average  values  during  run-arrest-event) 

are  larger  than  Kjg  reflecting  a  rising  K^-crack  velocity  curve. 

Propagation  and  arrest  in  the  7075-T651  test  specimens  (see  Figure  5) 
is  complicated  by  the  relatively  weak  grain  boundaries  parallel  to 
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Figure  4.  Effect  of  face  grooves 
on  Kjr  (Aa“l.  50mm)  values  ob¬ 
tained  from  stable  crack  growth 
measurements. 
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Figure  5.  Appearance  of  fractured 
7075-T651  aluminum,  compact  crack 
arrest  test  specimens. 
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Figure  b.  Appearance  of  fracture  surfaces  7075-T651  Aluminum  produced 
by  (a):  slow  stable  crack  extension  in  a  conventional  J-test  and  (b) 
fast  fracture  in  a  compact  crack  arrest  test. 


TABLE  2.  RESULTS  OF  STABLE  CRACK  GROWTH  AND  CRACK  ARREST  MEASUREMENTS  IN  707S-T6S1<«>  ALUMINUM 
COMPACT  SPECIMENS  WITH  0%,  251  AND  801  FACE  GROOVES. 


STABLE  GROWTH 


CRACK  ARREST 


FACE  GROOVE 

t 

KJC,MP*  M1/2 

Kjr.MP. 

■1/Z  Kq,HP*  «172 

At.  m  Ku 

MPt  m1'* 

r1D,MPt  i 

0 

38 

50 

139 

58 

56 

86 

0 

38 

50 

118 

53 

54 

79 

25 

39 

49 

139 

51 

64 

94 

25 

39 

49 

139 

52 

57 

89 

25 

39 

49 

174 

>88 

<84 

<121 

80 

40 

— 

118 

56 

60 

83 

(a)  25»- thick.  707S-T651  plate:  o0  *552  MPa,  HRB  »  91;  with  cracks  extending  i  to  the 
LT-dl recti on. 

(b)  Percent  of  cross  section  reaioved  by  face  grooves. 

(c)  -  [Kq  .  A,.]"* 

the  rolling  plane.  Tensile  stresses  In  the  short  transverse  direction 
Initiate  cracks  in  these  boundaries  (delaminations)  at  the  root  of  the 
face  grooves  and  these  permit  the  fracture  path  to  deviate  from  the 
specimen  symmetry  plane.  In  addition,  the  high  speed  fracture  surfaces 
become  inclined  to  the  tensile  axis.  Fractographs  of  these  surfaces 
display  a  step-like  appearance  (Figure  6a  and  6b)  with  the  vertical 
segments  of  the  steps  corresponding  with  the  delamination  of  the  weak 
boundaries.  Consequently,  the  fractographs  contain  evidence  of  changes 
in  the  details  of  the  fibrous  mode  mechanism  of  7075-T651  aluminum 
that  could  reduce  the  toughness  values  for  high  crack  velocities  in 
Figure  3. 
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Optimal  Discretization  of  Power 
Stress-Strain  Law  Curves 

A  criterion  for  optima I  discretization  of  power  stress-strain  law  curies  is  proposed. 
The  criterion  is  based  on  the  assumption  that  it  is  desirable  to  have  the  fewest 
possible  line  segments  without  exceeding  some  predetermined  bound  on  the  error. 
The  formulation  produces  a  system  of  simultaneous  nonlinear  equations  which  are 
solved  using  an  iterative  search  technique.  Solutions  are  presented  in  both  graphical 
and  tabular  form  for  a  wide  range  of  strain  hardening  exponents  and  acceptable 
error  bounds.  It  is  shown  that  stress  and  energy  density  can  be  accurately  and  ef¬ 
ficiently  modeled  using  the  optimal  discretization. 


Introduction 

The  power  stress-strain  law  is  frequently  used  in  the 
mathematical  analysis  of  plasticity  problems  since  it  is  a 
reasonable  description  of  the  behavior  of  many  materials 
beyond  first  yield.  One  of  the  most  common  applications  of 
the  power-law  is  in  elastic-plastic  fracture  mechanics  since  the 
stresses  and  strains  at  the  tip  of  a  stationary  crack  were 
determined  by  Hutchinson  [1],  and  Rice  and  Rosengren  [2]  on 
the  basis  of  a  power-law  maierial. 

The  finite  element  method  (FEN1)  is  widely  applied  in  the 
study  of  elasto-plasiic  fracture  mechanics.  Whenever 
possible,  the  user  should  code  the  actual  exact  analytic 
representation  of  the  constitutive  relation  into  the  FEM  code 
as  a  part  of  a  material  library.  Neverthe'ess,  there  may  be 
occasions  when  this  is  impossible  or  impractical,  notably 
when  certain  commercial  FEM  software  is  used.  A  common 
approach  when  this  occurs  is  to  approximate  the  stress-strain 
law  with  a  series  of  discrete  line  segments.  Until  now.  how  to 
perform  the  discretization  has  been  largely  a  matter  ot  art  and 
experience.  Many  strategies  offer  plausible  solutions.  Ex¬ 
perience  suggests  that  the  greatest  success  can  be  obtained  by 
using  more  segments  on  the  portion  of  the  curve  with  greatest 
cirvature.  In  this  note  we  suggest  an  objective  criterion  for 
the  discretization  based  on  the  premise  that  it  is  desirable  to 
have  as  few  segments  as  possible  without  exceeding  some 
predetermined  bound  on  the  error. 

In  order  to  accurately  describe  a  smooth  curve  such  as  the 
power  stress-strain  law,  a  large  number  of  line  segments  is 
typically  required.  However,  a  large  number  of  line  segments 
in  the  discretized  stress-strain  curve  is  computationally  un¬ 
desirable  since  this  increases  the  CPU  time  required  to  sort 
and  interpolate  in  order  to  find  the  stress  corresponding  to  a 
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Fig.  1  Discretized  power-law  stress-strain  curve 


particular  value  of  strain.  This  note  describes  a  technique 
whereby  the  power-law  stress-strain  curve  is  discretized  using 
the  minimum  number  of  line  segments  possible  so  that  a 
predetermined  error  criterion  is  not  violated.  Results  are 
presented  for  a  large  range  of  values  of  the  strain  hardening 
exponent,  rt,  for  cases  where  the  stresses  are  held  to  within 
±0.1,  ±0.5  or  ±1.0  percent  of  the  exact  value. 

Problem  Formulation 

Consider  a  power-law  stress-strain  curve  of  the  form 

( /f,  =  o/o,  as  a,  ( lu  1 
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e/«0  =(a/o0Y  a>aa,n>\  (lb) 

Let  e  and  s  be  normalized  quantities  so  that  e  =  t/ta  and 
s  =  o/o0.  Equations  (1)  can  then  be  rewritten 

e  =  s  j£l  (2a) 

e  =  s "  s>l  (2b) 

These  relations  can  be  more  conveniently  addressed  if  we 
define  N=  l/n  so  that 

s  =  e  esl  (3a) 

s=ev  e>\,  0<A/<  1.  (3b) 

If  the  allowable  relative  error  is  denoted  by  6,  the  line 
segments  which  approximate  equation  (3b),  s',  must  be 
within  the  limits 

(1 -6)s<s*<(l +6)5.  (4) 

Clearly  the  longest  line  segment  which  will  be  bounded  by 
equation  (4),  as  shown  in  Fig.  1,  will  intersect  the  (l+6)s 
curve  at  points  /,  and  /  +  1 ,  and  be  tangent  to  the  ( 1  -  6)s  curve 


at  some  intermediate  point  /  +  £.  i.e., 
at  point  /:  5*=(l+6)e;v  (5a) 

at  point /+1:  s% ,  =  ( 1  +  &)e\ i  (5b) 

at  point /+£:  s'.i  =  (l  -  6)e,\(  (5c) 

=  (l-6)[(l-|)e, 

But  since  points  i,  t  +  £,  and  i  +  1  all  lie  on  a  straight  line, 

j:_s = (i <6a> 
which,  upon  substitution  fors'ands*,,,  becomes 

5\;  =(1  +a)[(l  +?e, v,1].  (6b) 


Since  /  +  £  is  a  point  of  tangency,  the  value  of  s'..  as 
calculated  from  the  (1  -6)s  curve  (5c)  must  equal  the  value 
calculated  from  the  line  segment  (6b): 

( 1  +  A)[(  1  -  £)*,'  +  ie\ ,}-(!-  6)l(  1  -  t)e, 

+  £<?,.,]' =0.  (7) 

The  slope  of  the  (1-6)5  curve,  calculated  at  the  point  of 

tangency  is 

ET=d[(\-b)es\/de\,.i  (8) 

=  M1 -«)[(! -«)*,  +  $*, 

Likewise,  the  slope  of  the  line  segment  in  this  region  is 

E'r  =  (s’. , -sy(e,.l-e,).  (9) 

Since  the  slopes  must  be  equal  at  the  point  of  tangency, 

(5%,  -sy(e,.  i -e,)-N(\  -6)[(l-i)e, 

+  =0,  (10a) 

or, 

( 1  +  6)(e\ ,  -  e,s )  -  M 1  -  &)(e, . ,  -  e, )(( 1  -  { )e, 

+  ?e,.,r-|=0.  (10b) 
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Fig.  2  Values  of  strain  at  the  end  of  the  first  line  segment,  a, .  for 
discretized  power-law  stress-strain  curve 
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Fig.  3  Strain  ratio,  «  =  e(„  t  le,,  for  discretized  power-law  stress-strain 
curve 


In  order  for  the  discrete  approximation  to  exhibit  the  same  three  unknowns  e ] ,  e:  and  f  (if  /  =  1 )-  Unfortunately,  these 
yield  point  as  the  power-law  curve,  and  in  order  to  mitigate  equations  are  nonlinear  and  no  closed  form  solution  is  ap- 

errors  for  small  plastic  strains,  we  suggest  that  the  first  line  parent.  These  equations  can  nonetheless  be  solved  using  an 

segment  should  begin  tangent  to  the  power  law  curve.  iterative  search  strategy.  An  objective  function  is  constructed 

Therefore,  the  first  line  segment  is  defined  by  rules  that  differ  by  summing  the  squares  of  equations  (7),  (10b),  and  (13): 


from  the  approximation  strategy  used  for  other  points. 

The  slope  of  thes  =  es  curve  at  the  yield  point  is 

ET  =  d(es)/de\t^=N.  (II) 

The  slope  of  the  first  line  segment  is 

£*T=(s;-l)/(e, -1).  (12) 

Since  these  slopes  are  equal, 

(1 +6)^-1 -<V(e,  -  1)  =  0.  (13) 


/(el,e:.{)=l(l+6)[(l-0t’,'  +  i<?:'l 

-(l-6)[(l-{)e, +{e2Jv|:  +  ((l+«)(e2v 
-e,  V)-M1  -«X<?2 -<?,)((!  - €)e,  +£e2]'-'  l: 

+  i  (1  +  6)c,  v  —  1  -  Af(e,  -  1)J!-  (14) 

This  function  is  equal  to  zero  at  the  solution  to  equations  (7), 
(10b),  and  (13).  The  function, /(e, ,  e:.  0,  was  minimized  for 
several  values  of  N  and  6  using  the  BFGS-Armijo  variable 


Equations  (7),  (10b),  and  (13)  are  three  equations  in  the  metric  algorithm  [3-7]  as  implemented  by  Santi  et  al.  [8-9). 
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Fig.  5  Error*  in  *n*rgy  d*mity  lor  *  *  »°  ’  with  t «  0.5  p«rc*nt 

Convergence  was  assumed  when  the  values  for  e,.  e:,  and  £ 
did  not  charge  in  the  eighth  decimal  place.  Execution  times 
were  typically  5  seconds  on  a  DEC-1099  computer. 

It  is  a  simple  matter  to  reformulate  equations  (7)  and  <10/>) 
into  independent  relations  of  the  form/(©,£)  =  0  where  ©  = 
e,.s/e,.  Since  these  relations  admit  but  one  solution  for  ©  and 
once  e,  and  e2  are  known,  all  other  points  on  the  discrete 
stress-strain  curve  can  be  determined  using  the  recurrence 
relationship 

e:/e{  =  e,/e:  =  ...  =  e,.,/e,  =  ...  =©  (15) 

Other  formulations  of  this  problem  are  possible.  The 
present  approach  is  straightforward  and  computationally 
tractable. 

Results 

Equation  (14)  was  minimized  for  15  different  values  of  A', 
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and  for  values  of  6  =  0.1,  0.5,  and  1.0  percent.  Table  1 
presents  results  for  c. ,  the  ratio  o  =  e..t/e,,  and  £.  All  data 
have  been  verified  as  correct  and  none  produce  errors  greater 
than  10  '  percent  outside  of  the  (1  ±6)s  range.  Values  of  e,. 
o,  and  £  are  plotted  in  Figs.  2-4,  respectively. 

Attempts  were  made  to  determine  if  a  linear  relationship 
could  be  found  between  any  of  [e, ,  e:,  O.  £,  log  e, ,  log  e:,  log 
O,  log  £1  and  any  of  In,  <V,  \  ,V,  log  ,V|.  without  success.  If 
results  are  desired  for  values  of  .V  not  listed  in  Table  I,  then 
an  appropriate  interpolation  scheme  should  give  quite  good 
values  of  *|  and  o.  Stress-strain  curves  produced  in  this 
manner  will  still  be  optimal,  but  the  error  term,  6,  may  not  be 
exactly  the  same  as  stated  in  Table  1 . 

Another  frequently  overlooked  aspect  of  stress-strain  curve 
discretization  is  that  the  energy  density  is  poorly  represented  if 
a  number  of  points  on  the  stress  strain  curve  are  used  to 
discretize  the  stress-strain  curve.  The  energy  density  is  the  sum 
of  the  strain  energy  density  and  the  plastic  dissipation  at  a 
given  material  point: 

W'=  \  o  dt  ..  (16) 

Jo 

The  optimal  stress-strain  discretization  technique  gives  a 
much  better  representation  of  the  energy  density  than  would 
be  obtained  if  a  number  of  equal  increments  were  used  to 
discretize  the  stress-strain  curve.  Figure  5  depicts  a  rather 
extreme  example  of  this.  The  curve  s  =e"  1  Is  discretized  using 
six  line  segments  in  the  range  l<e<IOO.  In  one  case  the 
optimal  discretization  technique  is  used  with  six  intervals  (i.e.. 
6  =  0.5  percent),  and  in  the  other  case  six  equal  intervals  of  e 
are  used.  The  errors  in  energy  density  are  plotted  in  Fig.  5. 
The  optimal  discretization  technique  produces  errors  in 
energy  density  of  roughly  -0.2  percent,  whereas  the  equal 
interval  technique  produces  errors  of  as  high  as  -  6.6  percent . 
Indeed,  the  equal  interval  technique  also  produces  a  very  poor 
representation  of  the  stress  as  a  function  of  strain,  where 
errors  as  large  as  -8.1  percent  occur  at  small  values  of  strain. 
This  is  depicted  in  Fig.  6. 


Conclusions 

An  objective  technique  was  described  whereby  the  power- 
law  stress-strain  curve  was  represented  by  an  optimal  com¬ 
bination  of  line  segments.  Results  were  presented  for  a  large 
variety  of  strain  hardening  exponent  values,  n,  at  each  of 
three  allowable  deviation  values,  6.  A  detailed  investigation  of 
one  case  revealed  that  stress  and  energy  density  are  accurately 
modeled  using  the  optimal  discretization  technique. 
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INTRODUCTION 

As  part  of  a  study  of  the  crack  arrest  capabilities  of  tough 
steels[l],  efforts  are  underway  to  simulate  rapid  crack  extension 
and  arrest  in  elastic-plastic  finite  element  models.  As  a  first 
step,  stationary  cracks  in  compact  tension  specimens  have  been 
modelled  and  the  effects  of  loading  rate,  strain  rate  sensitivity 
and  inertia  on  Jj  have  been  examined.  The  aim  of  this  work  is  to 
examine  those  features  of  the  plastic  zone  influential  in  deter¬ 
mining  the  toughness,  namely,  the  size  of  the  process  zone,  and 
the  crack  tip  opening  displacement. 

Plasticity  associated  with  a  stationary  crack  has  been  charac¬ 
terized  as  occurring  in  two  separate  zones,  as  shown  in  Fig.  1.  The 
larger  region,  called  the  plastic  zone,  features  small  plastic 
strains  in  the  range  0  <  Ep  <  0.1.  The  size  of  the  plastic  zone 
(at  0  *  0°)  is  given  by  a  characteristic  dimension,  rG.  Levy,  et 
a  1 - [ 2 ]  have  determined  r0  for  a  nonhardening  material  with  a 
semi-infinite  crack  in  an  infinite  plate.  In  cases  where  the 
plastic  zone  is  small  compared  to  the  specimen  dimensions 

rQ  -  0.036  EJI/oo2(l-  v2).  (1) 

Even  closer  to  the  crack  tip  is  the  heavily-strained  process 
zone,  where  plastic  strains  range  roughly  from  0.1  to  1.0.  This  in¬ 
tensely  non-linear  zone  can  be  given  a  characteristic  dimension,  w, 
which  can  be  related  to  the  crack  tip  opening  displacement,  6,  by 
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Rice[3]  suggests  1.9  as  a  possible  value  of  B,  and  Paris[4]  states 
It  Is  of  the  order  of  2.  Since  the  process  zone  Tias  no  “obvious" 
boundary,  the  value  of  B  depends  on  how  the  process  zone  is  defined; 
for  instance,  a  critical  strain  value  may  be  chosen  to  define  the 
process  zone  boundary.  For  convenience,  B  m  1  has  been  used  in 
this  analysis. 

The  crack  tip  opening  displacement,  6,  can  be  expressed  as  a 
function  of  the  J-integral[5]  by: 

6  -  dnJ/oc,  (3) 

for  a  nonlinear  elastic,  power  law  hardening  material.  The  constant, 
dn,  is  a  function  of  material  properties,  a,  oQ,  E  and  n  (Ramberg- 
0sgoodI6]),  where 

e/c0  -  a(o/oo)n.  (4) 

For  small  values  of  strain,  e  can  be  interpreted  as  the  total 
strain;  for  large  values  of  strain,  e  can  be  interpreted  as  the 
plastic  strain.  Using  values  for  a  steel  such  as  A533B,  a  value  of 
dn  can  be  determined [5 ]  as  dn  ■  0.52.  Equation  (3)  can  be  substituted 
into  (2)  and  then  divided  by  (1)  giving  the  relative  size  of  the 
process  and  the  plastic  zones: 

w/rQ  -  13.14  o0/E .  (5) 

If  E  *  197  GPa  and  o0  *  415  }Pa  as  in  the  case  of  A533B  steel  at 
93°C,  then 

w/rQ  *  0.03.  (6) 

The  plastic  strain  rates  are  related  to  the  strain  gradients 
and  the  extent  of  the  plastic  and  process  zones.  If  it  is  assumed 
that  the  equivalent  plastic  strain  is  0.0  at  r  ■  r0  (the  plastic  zone 
boundary),  and  varies  linearly  up  to  a  value  of  0.005  at  r  *  w  (the 
process  zone  boundary),  then  the  plastic  strain  can  be  expressed 
as: 


ep  *  0.0Q5(ro  -  r)/(rD  -  w).  (7) 

After  making  appropriate  substitution  from  (6)  and  (1),  and  dif¬ 
ferentiating  with  respect  to  time: 

ep  -  0.005(r/ro)(JI/J1).  (8) 

An  "average"  value  of  plastic  strain  rate  could  be  calculated  at 
the  middle  of  the  plastic  zone,  i.e.  where  r  «  0.5ro. 
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This  gives 

I  -  0.0025(iI/JI)  -  O.OOSCKj/Kj) 


(9) 


for  the  plastic  zone,  and  is  equivalent  to  the  ep  equation  of 
Wllson[14] . 

Another  expression  for  the  strain  rate  in  the  plastic  zone 
can  be  derived  from  the  J-field  solutions  of  Hutchison{7] ,  Rice 
and  Rosengren[8] .  Hie  equivalent  plastic  strain  is 


e  “  (o0/E)[(EJ)/(oo02Inr)]n/n+1  efen(n,6). 


(10) 


(Note  that  equation  (10)  assumes  elastic  strains  to  be  negligible.) 
The  constants  In  and  a  are  given  in  Refs.  {8,9],  and  the  remain¬ 
ing  values  are  defined  by  the  Ramberg-Osgood  model  of  eq.  (4). 
Differentiating  (10)  with  respect  to  time,  and  substituting  the 
appropriate  values  for  A533B  (a  *  1.12,  n  ■  9.71) [10]  at  0  ■  0, 
yields 


ep  -  3.24  x  IQ”11  r-0.9066J-0.09337i> 


(11) 


Again,  an  "average"  value  of  plastic  strain  can  be  calculated  in 
the  middle  of  plastic  zone  (r  ■  0.5ro)  giving 


m 


ep  -  0 ,00024( j/J)  -  0.00048  (fcj/Kj).  (12) 

Hie  plastic  strain  rate  predicted  by  (9)  is  10  times  larger  than 
that  predicted  by  (12).  Hiis  is  due  to  (i)  the  nature  of  the  linear 
approximation  used  in  (9),  and  (ii)  the  fact  that  the  J-fleld  is 
not  a  good  approximation  where  plastic  Btrains  are  small  (such  as 
in  most  of  the  plastic  zone).  It  is  not  clear  which  model  is  better 
or  more  reliable. 

A  similar  linear  approximation  could  be  used  to  determine  the 
order  of  magnitude  of  plastic  strains  in  the  process  zone.  Here 
the  assumptions,  that  ep  ■  0.005  at  r  ■  w  and  ep  ■  0.25  at  r  «*  0, 
are  suggested  by  finite  element  results  of  M:teeking[ll]  for  a 
material  where  n  •  10.  Typical  strain  rates  in  the  center  of  the 
process  zone  for  8*0  are 

ep  -  O.SStfj/Jj)  -  O^Kj/Kj).  (13) 

The  plastic  strain  rates  can  directly  be  calculated  for  8*0  using 
the  J-field  solution  of  (11).  Hie  result  given  below  is  not  reliable 
since  r  ”  w/2  {18] . 

Ep  -  0.007( jj/Jj)  -  O.OHtfj/Kj).  (14) 

In  view  of  equations  (9)  and  (13),  plastic  strain  rate  can  be 


determined  at  the  onset  of  crack  extension  of  ductile  materials  such 
as  A533B.  Table  1  gives  typical  values  of  plastic  .strain  rates  for 
"slow"  and  "fast"  loading  rates.  Later,  these  shall  be  compared 
with  finite  element  computations. 

Ductile  crack  extension  proceeds  by  void  nucleatlon  and  growth 
in  the  process  zone.  Toughness  is  controlled  by  void  spacing  which 
determines  wc,  the  critical  size  of  the  process  zone.  From  equa¬ 
tions  (2)  and  (3)  one  can  deduce  Jjc  in  terms  of  the  critical  process 
zone  size 

Jlc  ■  wco0/6dn  -  2o0wc.  (15) 

Changes  in  the  yield  stress,  for  constant  values  of  wc,  will  alter 
Jjc  according  to 

■  «ob'«o.-  (16> 
The  subscripts  a  and  b  in  (16)  refer  to  two  different  loading  rates. 

Experimental  data  is  available  from  a  number  of  sources [12-16] 
which  relate  the  variation  in  Kjc  to  the  loading  rate  fcj.  Some  of 
these  data  are  plotted  in  Figure  2.  Using  the  strain  rate  expression 
(9),  and  the  yield  stress  data  given  in  Ref s . [12-16 ] ,  comparisons 
are  made  between  the  changes  in  yield  stress  and  the  changes  in 
Kj,  as  a  function  of  Kj-rate.  These  comparisons  are  given  in 
Table  2,  Equation  (16)  suggests  that  the  last  two  columns  in 
Table  2  should  have  equal  values.  In  fact,  good  agreement  is 
achieved  only  for  AISI  1018.  The  lack  of  agreement  may  have  2 
sources.  The  first  one  is  connected  with  the  relative  contribution 
of  the  flow  stress  in  the  plastic  and  process  zones.  The  yield 
stress  ratios  in  Table  2  are  based  on  the  strain  rates  in  the 
plastic  zone.  Since  the  strain  rates  in  the  process  zone  are  10^ 
to  10-*  times  the  values  for  the  plastic  zone,  the  yield  stress 
ratios  appropriate  for  the  process  zone  would  be  larger  for  mate¬ 
rials  whose  rate  sensitivity  increases  with  strain  rate  (see  Figure 
Ab).  Consequently,  a  significant  contribution  by  the  resistance 
to  flow  in  the  process  zone  could  account  for  the  greater  rate 
sensitivity  of  the  Kjc-values  of  A533B  and  AISI  1020.  The  second 
source  is  the  rate  sensitivity  of  wc,  which  Is  neglected  by  Equa¬ 
tion  (16).  The  increases  in  the  flow  stresses  and  normal  stresses 
associated  with  the  higher  strain  rates  can  facilitate  void  nucle- 
ation  and  reduce  the  void  spacing  and  wc.  A  large  reduction  in 
wc  in  the  face  of  the  increase  in  o0  could  account  for  the 
relatively  low  Kjc-ratios  displayed  by  AISI  1018  and  the  PA-6 
aluminum  alloy. 

The  existing  analyses  do  not  clearly  distinguish  the  relative 
contributions  of  the  changes  in  the  plastic  and  process  zones.  The 
fact  that  the  rate  of  strain  hardening  (the  value  of  n)  influences 


both  5  and  w[5,ll],  when  the  large  strains  and  hardening  are  largely 
confined  to  the  process  zone,  is  evidence  that  the  .rate  sensitivity 
of  flow  in  the  process  zone  cannot  be  Ignored.  The  finite  element 
calculations  described  in  the  following  sections  were  undertaken  to 
illuminate  this  issue. 

FINITE  ELEFENT  MODELLING 

A  stationary  crack  in  a  compact  tension  specimen  has  been 
modelled  using  finite  elements.  Figure  3  shows  a  plot  of  the  mesh. 
To  limit  the  complexity  of  the  model,  the  elements  nearest  the  crack 
tip  were  sized  to  1  mm.  The  mesh  was  composed  of  8-noded  isopara¬ 
metric  elements.  Numerical  studies  by  deLorenzi,  Sh ih [19],  Hoff 
and  Byrne [20]  have  shown  that  these  elements  are  suitable  for  use 
in  the  fully  plastic  range. 

The  model  shown  in  Fig.  3  has  3500  degrees-of-f reedom,  and  con¬ 
sequently  the  analysis  times  for  the  elastic-plastic  problem  are 
very  long.  As  a  preliminary  step,  the  same  mesh  configuration  was 
used  in  an  analysis  employing  4-noded  isoparametric  elements, 
which  resulted  in  execution  times  of  about  one-tenth  of  the  8-noded 
analysis.  The  4-noded  isoparametric  elements  tend  to  be  too  stiff 
in  the  fully  plastic  range,  as  discussed  in  Ref.  [19],  and  this 
will  be  borne  out  in  the  results. 

The  mesh  was  constructed  so  that  the  stationary  crack  analysis 
could  later  be  extended  into  the  advancing  crack  regime.  The  method 
of  crack  extension,  as  implemented  by  Kanninen  et  al.[17],  involves 
the  release  of  constraints  at  the  crack  tip.  This  technique  gives 
no  "obvious"  crack  tip  opening  displacement  for  the  stationary 

crack  problem.  The  authors  have  selected  the  separation  of  the 
crack  faces  1  mm  away  from  the  crack  tip  to  be  the  CTOD,  as  shown 
in  Figure  6.  A  more  refined  approach  could  be  effected  by  imposing 
displacements  on  the  boundary  of  a  region  very  close  to  the 
crack  tip,  as  has  been  done  by  Sorensen[22] . 

An  important  feature  of  the  analysis  is  the  incorporation  of 

the  strain  rate  dependence  of  the  flow  stress.  Considerable 

experimental  data  is  available  for  A533B  steel[23]  as  a  result  of 
testing  for  nuclear  applications.  The  uniaxial  6tress-strain  curve 
is  given  in  Figure  Aa [17],  and  corresponds  to  8  nominal  strain 

rate  of  The  relationship  between  the  instantaneous  flow 

stress,  o,  and  the  plastic  strain  rate,  ip,  is  given  by  the 
Ffelvern[24]  equation 

a/a  -  [c  /D] i/p  +  1.  (17) 

o  p 

The  parameters  D  and  p  must  be  determined  experimentally,  and  c0 
is  the  flow  stress  at  zero  strain  rate.  Using  the  data  in  Ref. 
[23]  at  93C,  the  rate  dependence  relationship  is  plotted  in  Fig.  4b. 


muHHiuniunii/  imnmrRinnnmiK  mil)  VIV.  V^ATLWV 


RESULTS 

9 

The  mesh  configuration  and  material  model,  described  in  the 
previous  section,  were  incorporated  in  6  analyses;  3  different 
loading  rates  were  employed  with  4-noded  and  8-noded  elements. 
Table  3  summarizes  the  analyses  which  have  been  conducted. 

Analysis  1  is  static  using  o-e  relations  for  A533B  for  an 
infinitesimally  slow  rate  of  loading.  A  displacement  is  applied 
to  the  load  point  in  Increments  of  0.2  mm  until  a  maximum  displace¬ 
ment  of  2.4  mm  is  reached.  The  maximum  equivalent  stress  in¬ 
tensity  factor,  Kj,  is  calculated  from 

Kj  -  [JE/(1-v2)]1/2.  (18) 

Analysis  2  is  a  dynamic  analysis.  A  constant  velocity  of 
0.008  m/s  is  applied  to  the  load  point,  until  a  displacement  of 
2.4  mm  is  reached.  This  results  in  an  almost  constant  K-rate,  but 
not  a  constant  J-rate. 

Analysis  3  is  also  a  dynamic  analysis.  A  constant  velocity  of 
0.8  m/s  (100  times  faster  than  Analysis  2)  is  applied  to  the  load 
point.  Previous  studies  using  8.0  m/s  as  a  loading  rate  had  been 
attempted,  but  the  inertia  of  the  model  caused  unwanted  vibra¬ 
tions  which  made  the  Interpretation  of  results  very  difficult. 

Results  for  the  static  analysis  are  plotted  in  Figure  5. 
Fig.  5a  shews  a  dimensionless  J-integral  as  a  function  of  a  dimen¬ 
sionless  load  line  displacement.  A  comparison  with  experimental 
results  by  Andrews  and  Shih[25]  reveal  that  the  numerically  deter¬ 
mined  J-integral  values  are  low.  Figure  5b  shows  a  dimensionless 
crack  tip  opening  displacement  as  a  function  of  a  different  dimen¬ 
sionless  load  line  displacement.  Results  from  the  present  study 
are  again  smaller  than  the  experimental  results  of  [25] .  However, 
the  finite  element  results  of  this  study  are  in  good  agreement 
with  the  finite  element  results  of  Shih[18].  Equation  (3)  suggests 
that  6  varies  linearly  with  the  J-integral.  Fig.  7  plots  this 
relationship  and  shows  remarkable  agreement  with  the  static  finite 
element  results.  The  degree  of  agreement  could  be  fortuitous 
since  the  definition  of  6  is  somewhat  arbitrary,  but  is  does 
suggest  that  the  definition  of  6  in  Fig.  6  is  reasonable. 

The  plastic  strains  in  the  crack  tip  vicinity  are  of  interest 
since  they  essentially  determine  the  plastic  strain  rates  in  a 
dynamic  analysis.  The  magnitude  of  the  plastic  strain,  as  expressed 
by  tensor  product 

ep  -  [(3/2)ePlj:ePlj]1/2,  (19) 

is  shown  in  Figure  8,  as  well  as  a  comparison  with  the  J-field 


solution  of  equation  (10)  at  an  angle  of  45°.  Also  shown  in 
Figure  8  are  results  from  the  slip  line  solution  of  Rice  and 
Johnson[26],  which  were  later  duplicated  by  hfc Ffeekingl 1 1 ]  using 
finite  elements,  for  a  blunt  notch.  A  comparison  shows  that  the 
present  finite  element  results  underpredict  the  plastic  strains 
within  a  region  of  about  36  from  the  crack  tip.  This  is  not 
surprising  since  the  mesh  is  not  sufficiently  fine  to  capture  this 
large  strain  gradient.  It  follows  that  the  present  finite  element 
analysis  only  shows  the  effects  of  plastic  strain  rate  in  the 
plastic  zone,  and  it  is  not  refined  enough  to  show  the  effects  of 
plastic  strain  rate  in  the  process  zone. 

The  effect  of  loading  rate  is  apparent  in  Fig.  9.  If  a  criti¬ 
cal  value  of  crack  tip  opening  displacement  is  chosen  as  a  failure 
criterion,  then  higher  values  of  J  can  be  tolerated  at  higher  load¬ 
ing  rates.  In  fact,  for  the  range  of  strain  rates  occurring  in  the 
plastic  zone,  as  given  in  equation  (9),  the  "average"  increases  in 
yield  stress  would  be  112  and  212,  as  given  by  the  >felvern  equation 
(17).  In  view  of  equation  (16)  one  would  expect  equal  increases  in 
J-integral  over  this  range.  Fig.  9  shows  that  the  J-integral  values 
increase  by  142  and  262  for  a  CTOD  of  0.28  mm.  Although  the  agree¬ 
ment  between  the  two  predictions  of  yield  stress  Increase  looks 
good,  it  is  possible  that,  even  at  these  relatively  low  loading 
rates  the  inertia  of  the  specimen  may  have  some  effect  on  the  CTOD. 
Therefore,  it  may  be  desirable,  in  future  analysis,  to  attempt  to 
separate  the  effects  of  inertia  and  strain  rate  sensitivity  on  the 
crack  tip  opening  displacement. 

The  equivalent  plastic  strain  rate  is  given  in  Fig. 10.  Results 
are  plotted  for  6  *  45°,  since  at  this  angle,  the  plastic  strains 
are  larger  and  comparisons  are  more  obvious.  Equation  (10)  can 
be  differentiated  with  respect  to  time  (at  0  ■  43°)  in  the  same 
sense  that  equation  (11)  was  derived.  Values  of  the  equivalent 
strain  rate  from  the  J-field,  in  Fig.  10,  are  in  reasonable  agreement 
with  dynamic  finite  element  results.  Differences  Increase  as 
the  plastic  strains  become  small,  since  the  J-field  is  only  valid 
where  the  elastic  strains  can  be  neglected.  Here  again,  the  J- 
field  solution  is  only  valid  over  a  very  small  region. 

CONCLUSIONS 

1)  The  loading  rate  sensitivity  of  the  fracture  toughness  for 
the  ductile,  fibrous  mode,  as  expressed  by  Jjc,  should  correspond 
with  the  rate  sensitivity  of  the  flow  stress  in  the  plastic  and 
process  zones,  when  the  critical  crack  tip  opening  displacement  is 
1  ixed. 

2)  The  1  mm,  8-noded  finite  element  mesh  used  to  model  the 
region  near  the  crack  tip  describes  the  plastic  zone,  but  is  not 


sufficiently  refined  for  an  adequate  representation  of  the  process 
zone  of  A533B  steel  at  the  Jj  ■  230  KJ/m  level. 

3)  Numerical  results  obtained  with  the  finite  element  model 
confirm  that  in  the  absence  of  a  process  zone,  the  increase  of  Jj 
with  loading  rate  at  constant  CTOD  corresponds  closely  with  the 
increase  in  the  flow  stress  in  the  plastic  zone. 

4)  Experimental  measurements  reported  in  the  literature, 
contain  examples  of  Jjc-values  with  a  greater  rate  sensitivity,  and 
values  which  are  less  rate  sensitive  than  the  flow  stress  in  the 
plastic  zone.  These  cases  are  associated  with  either  the  103-  to 
10^-fold  higher  strain  rates  generated  in  the  process  zone,  or  the 
possible  negative  rate  depending  on  the  void  spacing. 

5)  Values  of  the  equivalent  (total)  strain  rate  in  the  plastic 
zone  derived  from  the  J-field  agree  with  equivalent  plastic  strain 
rate  values  from  the  finite  element  model  in  the  range  36<r<0.5  r0. 
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TABLE  1:  Estimates  of  Plastic  Strain  Rates 


Loading 

J*I 

“i 

Vs'1! 

CpIs'1] 

Rate 

[kJ*m”2*  S~1 ] 

[  M’a'm^/^*  s"1  ] 

Plastic  Zone 

Process  Zone 

Slow 

2 

1 

3xl0-5 

4x1 0~3 

Fast 

2xl05 

105 

3 

400 

Tabulated  values  are  for  a  ductile  steel  such  as  A533B. 
JIc  -  180kJ/m2;  Klc  -  200  M>a*mI/2 


TABLE  2:  Comparison  of  the  Strain  Rate  Dependence  of  Yield 
Strength  and  the  Kj-Rate  Dependence  of  Kjc  for  Fibrous 
M>de  Crack  Extension  in  Several  Alloys 


Alloy 

Ref. 

Kj-Range 

[  M^m*  s- 

• 

c  p-Range 

l)  ts"1] 

(o  ob^ 
°oa^ 

1/2  ^leb^Ica 

A533B(0177C) 

12,13 

2-2xl05 

5x10-5-5 

1.09 

1.60 

A1SI  1020(HR) 

14 

l-2xl06 

6x10-5-80 

1.36 

1.65 

AlSI  1018(CR) 

15 

l-2xl06 

4x10-5-80 

1.14 

1.08 

PA6  Aluminum 

16 

0.5-2xl06 

6x10-5-350 

1.14 

0.66 

Subscript  a  refers  to  the  lowest  strain  rate  and  b  refers  to  the 
highest  rate. 


TABLE  3:  Summary  of  Finite  Element  Analyses 


Analysis  Analysis 

*L 

''Lmax 

JImax 

^Jmax 

Kj 

Number 

Type 

[m/s  ] 

[mm] 

[kJ /m^] 

[  hPa*  m*  ] 

[M>a»m^^/s] 

1 

Static 

0.0 

2.4 

245 

230 

0 

2 

Dynamic 

0.008 

2.4 

254 

235 

689 

3 

Dynamic 

0.8 

2.4 

260 

238 

71000 
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Equivalent  stress  [M’s) 


Figure  6:  Details  of  crack  tip  mesh,  showing  crack 
tip  opening  displacement,  6. 


Figure  7:  Relationship  between  crack  tip  opening 

displacement  and  J-integral  for  a  static 
finite  element  analysis  for  A533B  @  93C. 
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Figure  8:  Variation  of  plastic  strains  with  distance 
from  the  crack  tip  @  6  ■  45°, 


Figure  9:  Relationship  between  crack  tip  opening 

displacement  and  J-integral  for  statically 
and  dynamically  loaded  analyses. 
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Figure  10:  Relationship  between  plastic  strain  rate 
and  distance  from  crack  tip  for  0  »  45°, 
vL  -  0.8  m/s  and  Jj  ■  245  kJ/m. 
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The  Influence  of  Microstructure  on  Brittle  Fracture  Toughness 


G  T  HAHN 


This  lecture  presents  viewpoints  of  the  role  of  microstructure  in  brittle  fracture  that  have  emerged  in 
the  past  tw  o  decades  The  fracture  mechanics  concept  of  crack  arrest  is  inserted  into  the  Griffith  theory 
of  microcracks  to  describe  the  resistance  of  microstructural  boundaries.  The  implications  of  crack 
blunting  are  related  to  the  essential  role  of  carbide  particles  and  other  brittle  phases  in  the  steel 
Dislocation  pile-ups  and  the  femte  grain  size  are  accorded  a  diminished  role.  The  brittle  fracture  stress 
is  related  to  the  dimension  of  the  largest  "eligible"  parent  particle  and  the  resistance  experienced  w  hen 
the  advancing  microcrack  crosses  the  boundary.  “Eligibility"  is  connected  with  the  probabilities  of 
nucleating  a  crack  and  finding  a  boundary  with  minimal  resistance.  The  boundaries  of  carbide 
particles,  carbide  films  at  the  grain  boundaries,  brittle  inclusion  particles,  ferrite  grains,  pearlite 
colonies,  and  bainite  packets  are  possible  barriers  The  use  of  the  grow  ing  data  base  of  Griffith  energy 
values  —  or  equivalently,  the  local  arrest  toughness  value  —  to  identify  controlling  microstructural 
features  is  demonstrated  The  lecture  touches  on  the  relation  between  the  transition  temperature  and 
the  brittle  fracture  stress  and  draws  attention  to  the  latter's  dependence  on  the  size  of  the  stressed 
volume.  Finally,  the  role  of  the  stressed  volume  in  analyses  of  the  brittle  fracture  A\ -value  and  a 
statistical  treatment  that  clarifies  its  origin  are  discussed.  These  analyses  indicate  that  a  small  number 
of  the  largest  particles  may  have  disproportionate  influence  on  .  and  that  the  microstructural 
features  that  effect  "eligibility"  may  have  a  modest  effect  on  Kv . 
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outstanding  scientific  contributions  to  the  metallurgical  profession  b\ 
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and  engineering 
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I.  INTRODUCTION 

I  want  to  thank  the  ASM  for  the  privilege  of  presenting  the 
Edward  DeMille  Campbell  Lecture  This  lecture  provides  a 
few  moments  to  dwell  on  a  life  that  is  conventional  and 
inspiring.1 21  Campbell  was  a  professor  of  metallurgy, 
chemistry,  and  chemical  engineering  at  the  University 
of  Michigan  from  1890  to  1925.  He  was  director  of  the 
University’s  Chemical  Laboratory.  He  was  a  productive 
researcher — the  author  of  75  technical  papers'  —  with  pio¬ 
neering  work  on  the  interpretation  of  the  X-ray  diffraction 
patterns  of  martensite  to  his  credit  At  the  same  time,  he  and 
his  wife  Jennie  Ives  raised  and  educated  six  children.  In 
these  respects,  Campbell's  life  is  not  unlike  that  of  many 
other  distinguished  members  of  the  society  . 

What  makes  his  life  entirely  unique  is  that  Campbell  was 
blind  —  blinded  at  the  age  of  28  in  a  laboratory  accident. 
Campbell  fathered  children,  taught  metallurgy  and  chem¬ 
istry,  administered  a  laboratory,  unraveled  X-ray  patterns, 
and  published  all  but  three  of  his  papers  after  he  became 
blind.  Campbell’s  productive  life  reveals  the  enormous  re¬ 
serve  of  inner  strength  —  the  "toughness"  of  the  human 
spirit  —  that  resides  in  us.  Fortified  with  that  spirit  we  can 
overcome  our  share  of  tragedy,  jettison  our  disappoint¬ 
ments.  and  meet  the  challenges  of  our  time. 

Like  the  "toughness”  of  the  human  spirit,  brittle  fracture 
and  the  toughness  of  steel  are  fundamental  issues  that  have 
engaged  many  people  for  a  long  time.  Reports  of  brittle 
fracture  in  service  appeared  as  early  as  1879/  the  Charpy 
test  in  1905.  Current  explanations  of  brittle  fracture  and  the 
ductile-to-brittle  transition  have  their  roots  in  papers  pub¬ 
lished  by  Griffith'  and  Ludw  ik"  in  1920  and  192 1 .  The  actual 
observations  of  Griffith  cracks  i;  led  to  the  development  of 
detailed  mechanistic  treatments IJ  I"  in  the  decade  from  about 
1953  to  1963. 

In  the  two  decades  that  have  followed,  mechanistic 
studies  have  been  extended  to  spherodized  carbide.1 
pearlitic,2"24  bainitic,2'2'*  and  martensiti c:< *'  micro¬ 
structures.  McMahon  and  his  collaborators'  12  have  shown 
that  except  for  the  crack  path,  the  mechanisms  of  inter¬ 
granular  fracture  of  temper  embrittled  steels  and  trans- 
crystalline  (cleavage)  fracture  are  nearly  indistinguishable 
(Figure  2.3  and  Figure  2.6)  Knott  and  his  collaborators” 
have  linked  the  micromechanism  and  the  statistics  of 
particles2"  with  the  macroscopic  fracture  toughness.  In  addi¬ 
tion,  progress  has  been  made  in  analyzing  toughness  on  the 
atomic  scale  14 "  * 

Today,  brittle  fracture  problems  still  exist  and  important 
questions  remain  unanswered  For  example,  ways  of  modi¬ 
fying  the  operation  of  existing  nuclear  plants  are  currently 
being  sought  to  avoid  brittle  fracture  Toughness  specifica¬ 
tions  for  bridge  steels  are  controversial  '  The  tw  o  common 
measures  of  toughness  —  the  CVN-energy  used  in  material 
specifications  and  the  Aj, -value  used  in  design  —  resist  cor¬ 
relation.1'  The  microstructural  features  governing  brittle 
fracture  are  not  well  defined  for  many  steels  Even  the 
inherent  toughness  of  unalloyed  ferrite  in  the  absence  of 
carbides  and  other  particles  is  an  open  question 

I  acquired  an  enduring  interest  in  this  subject  in  1953  from 
Maxwell  Gensamer"*  who  faced  ship  plate  fracture  problems 
in  World  War  II  Since  that  time,  graduate  studies  at  MIT 
and  a  succession  of  fracture  prevention  projects  for  gas 
pipes,  rocket  casings,  and  nuclear  components  have  enabled 
me  to  view  the  field  from  time  to  time  The  invitation  to 
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prepare  this  lecture  was  an  opportunity  for  me  to  examine 
recent  work  and  to  reflect  on  how  the  views  of  micro- 
structural  involvement  have  changed  since  my  introduction 
in  1953.  What  follows  reflects  my  limited  comprehension  of 
the  vast  literature  of  this  field  —  for  example.  Kotilaimnen4 
alone  cites  over  400  references  in  a  recent  dissertation  on 
brittle  fracture  of  bainitic  steels.  I  apologize  to  both  readers 
and  contributors  for  the  omissions 


II.  MECHANISMS  OF  BRITTLE  FRACTURE 

Metallographic  observations  of  cleavage  microcracks 
in  mild  steels  were  the  key  that  led  to  detailed  mecha¬ 
nistic  treatments  in  the  period  1953  to  1963  The  work  ot 
Allen.  Rees.  Hopkins,  and  Tipler.  Low.1"’  and  the  group 
at  MIT10  "  12  revealed  Griffith-like  cracks  form  after  the 
onset  of  yielding1* 111  and  span  individual  ferrite  grains 
and  carbide  films  at  the  gram  boundaries  (see  Figure  I  > 
The  association  with  yielding  made  it  possible  to  connect 
nucleation  of  the  cracks  with  dislocation  pile-ups  as  pro¬ 
posed  by  Zener  in  1948"  and  elaborated  on  by  Stroh  ‘  and 
Petch1'  (see  Figure  2.1).  All  that  was  needed  to  make  the 
Griffith  theory-  meaningful  was  to  suppose  that  more  (rap¬ 
ture  energy  is  spent  disrupting  the  lattice  when  the  micro- 
crack  crosses  a  grain  boundary  than  in  the  grain  interior 
Consistent  with  this,  estimates  of  the  Griffith  fracture 
energy  term  y  are  ~  10*  larger  than  y  .  the  specific  surface 
energy  ,  which  is  a  rough  measure  of  y  in  the  interior  ot  a 
crystal  (see  Table  I).  In  fact,  grain  boundaries  musi  otter  a 
range  of  y-values  depending  on  the  misonentation  Sm-e 
fracture  is  associated  with  the  first  boundary  crossed  suc¬ 
cessfully.  the  weakest  link  and.  consequently,  a  boundary 
with  a  relatively  low  y-value  will  be  controlling 

Subtleties  of  the  problem  and  the  application  ot  the 
Griffith  Theory  are  more  easily  discussed  w  ith  the  language 
of  LEFM  (Linear  Elastic  Fracture  Mechamcsi  Since  the 
grain  boundary  region  is  presented  with  an  advancing 
(rather  than  a  stationary  )  microcrack,  fracture  occurs  when 
the  microcrack  fails  to  arrest  For  this  reason.  1  believe 
the  Griffith  energy  term  y  is  properly  viewed  as  an  arrest 
toughness  parameter44  in  this  case  G*  =  2y  and  A;  - 
|2£y/(l  -  i’:))"\  where  F  and  r  are  the  elastic  modulus 
and  Poisson's  ratio  (see  Table  h  The  superscript  H  is  added 
to  emphasize  that  G ?  and  At.  apply  to  a  relatively  smooth, 
microscopic  crack  crossing  a  boundary  offering  minin', a! 
resistance  As  illustrated  in  Table  1.  At.  values  are  mudi 
lower  thanjhe  conventional  crack  arrest  toughness.  A, 

50  MPa\  m.  which  applies  to  a  rough  macrocrack  that  o 
overcoming  all  the  boundaries  in  its  path  as  well  as  ductile 
ligaments 44  The  LEFM-form  of  the  Griffith  theory  e\ 
presses  cr,.  the  brittle,  cleavage,  or  transgranular-traciuie 
stress  in  terms  of  2a.  the  (micro)  crack  size.  a  crack  shape 
factor,  and  A’*: 

<7,  =  A/.  </3\  Tti/I  |  l.ij 

and  more  specifically 

it,  =  I  25  At  <!  '  |  lh| 

a.  =  0.798  At<  '  |k| 

where  d  is  the  diameter  of  a  penny -shaped  crack  and  i  the 
total  length  ot  a  crack  whose  depth  is  much  greater  than 
its  length  Since  the  microcrack  is  generated  attei  the 
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Table  1.  Summary  of  Microcrack  and  Macrocrack  Resistance  Values  in  Terms  of  y,  Af„.  Au.  and  A t. 
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specific  surface  energy.  n  is  the  shear  modulus  and  h  is  the  Burgers 
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'  onset  of  yielding,  the  probity  of  applying  LEFM  must  be 

»  questioned  The  use  of  Eq  (1|  can  be  justified  provided 

k  the  nucleated  crack  advances  rapidly  and  test--  the  condi- 

ition  for  extension  in  a  time  duration  too  short  for  significant 
yielding  of  the  ferrite  Dynamic  contributions  to  Eq  |l| 
«  may  also  be  negligible  since  little  kinetic  energy  is  reflected 

:  back  to  the  microcrack  from  the  far  boundaries 


The  theory  connects  fracture  resistance  with  the  micro¬ 
structure  because  the  critical  microcrack  dimensions  corre¬ 
spond  w  ith  either  the  carbide  film  thickness  try  I  or  the  ferrite 
grain  diameter  id  l.  This  was  demonstrated  by  Allen  and  co- 
workers  and  by  Low  4  The  1953  Allen  ci  al.  paper  shows 
that  the  refinement  of  carbide  films  by  manganese  additions 
is  accompanied  by  a  systematic  lowering  of  the  transition 
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MECHANISMS  OF  BRITTLE  FRACTURE 


Fig  2  —  Schematic  representations  of  the  mechanisms  of  brittle  frac¬ 
ture  according  to  various  authors  "  14 1,24 11  The  letter  B  identifies 
barriers  to  crack  extension 


temperature:  Low  demonstrated  that  the  brittle  fracture 
stress  varies  linearly  with  d ~°\ 

At  first,  the  dislocation  pile-up  was  considered  only  in 
the  nucleation  stage  and  not  in  the  analysis  of  the  fracture 
condition,  i.e..  microcrack  extension  across  a  boundary  as 
in  Figures  2.2  and  2.3.  This  was  changed  by  the  painstaking 
metallographic  studies  of  McMahon  and  Coheni:  which 
revealed  that  the  large  microcracks  observed  in  the  ferrite 
grains  were  invariably  connected  with  a  cleaved  carbide 
particle  located  somewhere  in  the  grain  or  on  the  surround¬ 
ing  boundary  .  Their  findings  imply  that  the  rupture  of  the 
carbide  is  an  essential,  intermediate  event  between  the 
formation  of  the  pile-up  and  the  cleavage  of  the  ferrite. 
Accordingly.  Smith1*’  and  Almond.  Timbres  and  Embury1 
formulated  models  in  which  the  pile-up  continues  to  con¬ 
tribute  to  the  crack  driving  force  after  the  carbide  crack 
nucleates  (Figure  2.4).  Accordingly,  the  expression  for  a,. 
the  fracture  stress  for  microcrack  extension  across  the  car¬ 
bide  ferrite  boundary,  contains  a  term  arising  from  the 
carbide  film  and  one  from  the  pile-up:1*' 


2  ( Kty 
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carbide  film  pile-up 

The  terms  rr„  and  r,  are  the  effective  shear  stress  and  the 
lattice  friction  stress,  respectively  Equation  |2]  involves 
two  microstructural  dimensions:  (i)  the  dislocation  pile-up 
size  which  is  equated  with  d.  the  ferrite  grain  diameter;  and 


(ii)  the  carbide  film  thickness.  c0.  These  analyses  seem  to 
ascribe  a  more  central  role  to  the  carbide  and  reinforce 
the  then  prevalent  view  that  the  refinement  of  both  the  grain 
size  and  carbide  dimensions  enhance  toughness.  However, 
neither  Smith  nor  Almond  et  al.  really  explain  why  the 
carbide  is  essential  to  the  nucleation  process.* 

•The  carbide  particle  is  not  an  essential  pan  of  the  fracture  condition 
since  <■  is  replaced  by  d  when  the  grain-size  femte  microcracks  are 
observed 


III.  CRACK  BLUNTING  AND  THE  ESSENTIAL 
ROLE  OF  HARD  PARTICLES 

The  distinction  between  the  static  and  dynamic  aspects  of 
the  fracture  mechanism  also  has  a  bearing  on  the  transition 
temperature  and  the  role  of  carbides  Rice  and  Thomson'4 
have  treated  the  origins  of  toughness  in  crystals  and  con¬ 
clude  that  iron  is  brittle  at  temperatures  close  to  0  K  and 
crack  resistant  at  higher  temperatures  The  brittle  behavior 
arises  when  bonds  rupture  before  dislocations  can  form  to 
blunt  the  crack;  the  tough  response  when  the  crack  blunts 
before  the  bonds  rupture.  Results  of  a  recent  computer  simu 
lation  of  a  crack  tip  in  an  atomistic  iron  lattice  at  0  K  and 
400  K  performed  by  Mullins  and  Dokainish'*’  are  in  accord 
with  this  view.  The  temperature  at  which  an  iron  crystal 
containing  a  stationary  crack  changes  its  behavior  from  brit¬ 
tle  to  tough  —  TT(K i,).  the  A",, -transition  temperature  —  is 
not  well-defined  by  these  analyses  and,  to  my  knowledge, 
has  not  been  measured  experimentally  The  TT(A',  l  for  Fe- 
3  Si  single  crystals,  about  - 100  °C.  can  be  deduced  from 
the  work  of  Pilkington  and  Hull*  (see  Figure  3).  Since  Fe-3 
Si  is  in  many  respects  more  brittle  than  unalloyed  ferrite, 
the  TT(Ki)  for  ferrite  crystals  must  fall  below  -100  C 
and  probably  below  -200  °C 

The  response  of  a  crystal  to  a  rapidly  advancing  crack, 
characterized  by  Ku.  differs  from  that  for  a  stationary  cravk 
because  there  is  less  time  for  thermally  assisted  generation 
of  dislocations  near  the  crack  tip  The  time  available  for 
blunting  a  slowly  loaded  stationary  crack  is  r  —  10  second' 


Fig  3 —  The  At,, -transition  for  Fe-3  Si  single  crystals  with  stanonarv 
cracks  after  Pilkington  and  Hull4’  and  an  estimate  of  the  corresponding 
A, .-transition  for  cracks  propagating  with  a  velocity  a  -  100  ms  based 
on  the  theoretical  model  of  Jokl.  Vitek,  and  McMahon  ”  The  dashed 
portions  are  based  on  the  extrapolation  of  the  measured  K,  -curve 
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In  contrast,  the  time  available  to  blunt  a  brittle  crack  ad¬ 
vancing  with  a  velocity  of  a  =  100  ms  is  /  -  10 b/a  = 
2  •  10  "  s  (10f>  represents  the  extent  of  the  highly  stressed 
tip  region  capable  of  generating  dislocations  and  b  = 
0.2  nm  is  the  Burgers  vector).  The  approximate  10i:-fo!d 
reduction  in  the  available  time  produces  a  ~  10'-fold  re¬ 
duction  in  the  dislocation  contribution  to  Ku  according  to 
the  theoretical  model  of  Jokl.  Vitek,  and  McMahon,"  This 
means  that  A|U-values  for  Fe-3  Si  crystals  remain  at  the 
1  MPaVm  level  displayed  by  Klt  below  the  transition 
for  temperatures  up  to  about  100  °C  and  possibly  higher 
(see  Figure  3).  In  other  words,  the  7T(Alu)  for  Fe-3  Si 
crystals  occurs  above  100  °C.  The  same  reasoning  applies 
to  microcracks  in  the  vicinity  of  grain  boundaries  although 
the  relevant  A*-values  are  likely  to  be  higher  than  Alu  for 
the  single  crystal 

Several  other  observations  support  the  view  that  A*  - 
values  for  both  Fe-3  Si  and  unalloyed  ferrite  reflect  tem¬ 
perature  independent  brittle  behavior  to  temperatures  above 
~  100  °C.  For  example,  the  constancy  of  the  brittle  fracture 
stress  of  Fe-3  Si  in  the  range  -200  °C  to  at  least  -‘-50  CC. 
reported  by  Griffith  and  Owen4"  (see  Figure  4).  taken  to¬ 
gether  w  ith  Eq.  1 1  ]  shows  that  A®  must  also  be  constant  over 
this  range  since  the  microstructura!  dimension,  d  (or  C0l. 
does  not  vary  with  test  temperature.  The  photomicrograph 
of  a  fast  propagating  crack  that  arrested  in  Fe-3  Si  at  22  CC. 
in  Figure  5.  shows  that  no  dislocations  are  produced  while 
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Fig  4  —  The  temperature  dependence  of  the  brittle  fracture  stress  of  Fe- 
3  Si  steel  measured  with  CVN-specimens  after  Griffith  and  Owen  " 


Fig  5  —  Micrograph  of  the  tip  of  a  fast  propagating  crack  in  Fe-3  Si  steel 
that  arrested  at  22  °C  after  Hoagland.  Rosenfield.  and  Hahn  *  The  section 
has  been  etched  to  reveal  dislocations  The  dark  etching  region  at  the  crack 
tip  reflects  dislocation  generation  and  crack  blunting  which  were  absent 
while  the  crack  is  propagating 
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the  crack  is  propagating,  but  that  the  crack  does  blunt  after 
it  arrests.4-'  Another  piece  of  evidence  is  that  At.- values  for 
steel  correspond  quite  closely  with  A,, -values  for  brittle  ce¬ 
ramics  bicrystals  and  polycrystals  (see  Table  1).  Lange4-  has 
identified  processes  responsible  for  fracture  energy  dissi¬ 
pation  at  the  grain  boundaries  of  the  MgO  bicrystals.  They 
include  the  formation  of  multifaceted  fracture  surfaces  with 
overlapping  cracks  with  some  plastic  flow  in  the  region  of 
overlap.  Similar  processes  are  probably  involved  in  steel 
The  TT(Kt)  for  ferrite  is  not  established,  but  it  seems  likely 
that  it  corresponds  closely  with  the  7T(A|J  for  commer¬ 
cial  grades  of  steel  which  occurs  in  the  range  50  CC  to 
100  °C.4'  49  As  noted  by  Rice  and  Thomson,'4  the  transition 
temperature  depends  sensitively  on  the  resistance  to  bond 
rupture.  Thus,  temper  embrittled  steels  with  weak  bonds  at 
the  grain  boundaries  can  display  brittle  behavior  to  tem¬ 
peratures  above  400  °C.': 

These  circumstances  provide  possible  explanations  for 
the  behavior  of  microcracks  and  the  essential  role  of  a  car¬ 
bide  or  other  hard,  brittle  particles  or  inclusions  in  the  brittle 
fracture  of  steel. 

(i)  Dislocation  pile-ups  against  ferrite  grains  cannot  ordi¬ 
narily  produce  microcracks  because  the  crack  nucleus  (in  the 
ferrite)  blunts  when  the  temperature  is  above  7Y(A|,  i; 
i.e..  >  -200  °C. 

(ii)  Dislocation  pile-ups  can  nucleate  cracks  in  a  carbide 
particle  because  the  7T(A]t)  of  Fe>C  is  above  300  CC  ' 
The  same  is  true  for  other  hard  brittle  particles  below 
their  7T(AI(>. 

(iii)  The  fracture  of  a  brittle  carbide  particle  presents  the 
adjacent  ferrite  with  a  rapidly  advancing  crack.  This  crack 
can  penetrate  the  ferrite  because  A® -values  remain  at  the  brit¬ 
tle  level  in  the  temperature  range  TltA®)  <  T  <  7T(A'*i. 

(iv)  The  microcrack  will  anest  at  the  grain  boundary  when 
the  Griffith  condition  (Eq.  [1])  is  not  satisfied.  Ordinarily, 
such  an  arrested  crack  will  blunt  above  7T(A,, )  and  can¬ 
not  reinitiate. 


IV.  THE  DIMINISHED  ROLE  OF 
THE  DISLOCATION  PILE-l  P  AND  THE 
ORIGINS  OF  GRAIN  SIZE  EFFECTS 

Curry  and  Knott19  have  pointed  out  that  the  Smith  pile-up 
theory  of  fracture  (Eq.  [2] )  does  not  really  predict  a  femte 
grain  size  dependence.  They  note  that  the  value  of  re(,  is 
limited  by  the  operation  of  dislocation  sources  in  neigh¬ 
boring  grains  and  cannot  exceed  the_value  at  yielding 
Ttff  =  <r  5A’  where  k\  =  0.33  MPa  \  m  is  the  Hall-Petch 
parameter.15  51  When  (dQik[  is  substituted  for  rcff  in  Eq.  [2). 
the  grain  size  dependence  of  <j  is  lost  Secondly,  they  note 
that  the  contribution  of  the  pile-up  term  in  Eq.  [2]  is  neg¬ 
ligible  for  reasonable  values,  i.e..  d  =  20  ^m,^0.1  < 
c0  <  1  n,  r,  =  100  MPa,  and  A®  =  2.5  MPa  V  m.  and 
Eq  [2]  reduces  to  Eq.  [ lc] .  Curry  and  Knott  propose  that 
the  strong  carbide  size  dependence  of  Eq.  [lc]  has  been 
incorrectly  assigned  to  ferrite  grain  size  because  c0  and  d 
frequently  go  hand-in-hand  as  shown  in  Figure  6(a).  To  test 
this  point,  they  formulate  the  apparent  grain  size  depen¬ 
dence  (replacing  r0  in  Eq  [lc]  with  its  dependence  on  d 
as  specified  in  Figure  6(a)),  and  compare  it  with  mea¬ 
surements  in  Figure  6(b)  for  different  values  of  the  fracture 
energy  term.  It  can  be  seen  that  the  curve  for  y  =  14  Jm  •’ 
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a  ferrite  matrix  is  cr*  -1500  MPa.  Since  this  value  is  well 
below  the  theoretical  strength,  the  fracture  of  the  carbide 
particle  must  be  assisted  by  a  local  stress  raiser  such  as  a 
defect  in  the  panicle  or  a  pile-up  Defects  cannot  be  the 
answer  for  small  carbide  particles,  d,  c(l  ^  1  /xm.  since  the 
minimum  Griffith  crack  diameter  is  d  =  0.7  /xm  (for  an 
assumed  value  K\,  =  1  MPa  /xm).  For  this  reason  pile-ups 
cannot  be  ruled  out. 


V.  THE  EVALUATION  OF  Kt  AND  THE 
IDENTITY  OF  THE  CONTROLLING 
MICROSTRUCTURAL  FEATURES 


(a) 


( b ) 


Fig  6  — Relations  among  the  bnttle  fracture  stress,  femle  grain  diameter, 
and  the  coarsest  observed  carbide  film  thickness  after  Currv  and  Knott 
(a)  Experimentally  determined  relation  between  c„.  the  coarsest  observed 
carbide  film  thickness,  and  d.  the  ferrite  grain  diameter  c,.  =  fid l.  and 
(i)  the  “apparent  '  </  "“ '  dependence  of  a  based  on  reported  measurements 
The  solid  lines  were  obtained  by  replacing  c(1  in  the  expression  <J,  = 
0  798  Ki  t\'  with  fid  \  for  different  values  of  Ki  The  fact  that  the 
data  are  consistent  with  a  linear  dependence  of  a  on  rf '' ' .  but  not  on  <T‘  ' 
supports  the  view  that  r<.  is  the  controlling  variable  The  value  Xf.  = 
2  5  MPa  m" '  ( y  =  14  Jm‘;)  gives  the  best  fit  The  data  in  (bl  are  derived 
from  Refs  18.  52-61 


(A®  = 2.5  MPa  Vm)  agrees  with  the  experiments.  The 
scatter  is  not  surprising  since  the  relation  between  c„  and  d 
(Figure  6(a))  can  hardly  be  precise.  The  fact  that  the 
dependence  of  a  on  d "°5  is  not  linear  also  supports  the 
view  that  c0  is  the  controlling  variable  in  these  cases  A 
bonafide  d~° 5  ferrite  grain-size  dependence  would  still  be 
displayed  by  steels  with  ferrite  grain-size  microcracks,  as 
shown  by  Low.” 10 

These  findings  do  not  alter  the  view  that  microcrack  nu- 
cleation  is  assisted  by  slip  and,  possibly,  even  by  pile-ups. 
This  has  been  clarified  by  Argon,  Im.  and  Safoglu62  and  by 
Lindley,  Oates,  and  Richards6'  who  show  that  cr*.  the  peak 
tensile  stress  acting  locally  at  the  boundary  of  a  non-yielding 
particle  in  a  plastic  matrix.  <r*  =  ( 1  /d)  aQ  am.  is  the  sum 
of  terms  involving  1/d.  the  particle  aspect  ratio  and  yield 
stress  and  the  mean  stress.  In  other  words,  the  tensile  stress 
on  the  particle  is  amplified  by  the  onset  of  yielding  particu¬ 
larly  in  the  presence  of  a  larger  hydrostatic  component. 
Studies  by  Argon  and  Im64  and  Gurland6'  indicate  that  the 
local  tensile  stress  needed  to  fracture  cementite  particles  in 
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The  value  of  Ky  is  defined  by  measurements  of  the  brittle 
fracture  stress  and  the  dimension  of  the  particle  whose 
boundary  limits  the  microcrack  size  together  with  the  appro¬ 
priate  form  of  the  Griffith  Equation  (Eq.  [lb]  or  [  1  c] )  In 
addition  to  the  carbide  particle,  film,  or  lamella,  other  hard 
particles  or  inclusions  (Figure  7(a)),  fracture  may  be  con¬ 
trolled  by  one  of  the  following  features  depending  on  the 
microstructure  of  the  steel  (Figure  7(b)): 

(i)  ferrite  grain9 11 

(ii)  pearlite  colony21"24 

(iii)  bainite  packet25"28 

(iv)  martensite  packet. 25,29  30 

Particle-size  microcracks  have  been  observed  in  all  of  these 
microstructural  features  except  in  martensite  packets.  The 
microcracks  in  the  pearlite  nucleate  in  individual  carbide 
lamellae  in  linear  arrays.  They  are  joined  by  the  failure  of 
the  intervening  ferrite  and  form  colony-  and  prior  austenite 
grain-size  cracks23  24  (see  Figures  7(c)  and  7(d)).  It  is 
not  clear  whether  the  critical  dimension  of  the  microcrack 
that  determines  a-  is  the  colony  or  austenite  grain  di¬ 
ameter  as  assumed  by  Rosenfield.  Votava.  and  Hahn.21  or 
the  lamella  thickness,  as  assumed  by  Park  and  Bernstein" 
although  correlations  with  the  7T  (CVN)  favor  the  former. 
The  Af0-value  consistent  jxith  lamella  thickness.  A®  =  1.5 
MPa  vm.  is  smaller  than  the  other  values  in  Figure  8. 

The  interpretation  of  brittle  fracture  stress  values  of  dif¬ 
ferent  microstructures  involves  two  complications.  One  is 
that  the  particle  (or  grain)  generating  the  microcrack  is  rep¬ 
resented  by  a  distribution  of  sizes.  The  other  is  that  not  all 
particles  are  "eligible".  "Eligibility"  is  here  connected  with 
the  probability  a  given  particle  can  produce  a  crack:  (a)  that 
it  is  associated  with  a  suitable  pile-up  or  stress  raiser,  (b)  that 
it  has  an  orientation  favorable  for  nucleating  a  crack,  and 
(c)  that  the  misorientations  at  the  particle  boundary  produce 
a  low  value  of  A&.  For  example,  if  the  probability  of  satis¬ 
fying  each  of  these  requirements  is  p  =  0.1,  only  1  of  10' 
particles  is  “eligible”.  It  is  the  dimension  of  the  largest 
"eligible"  particle  within  the  stressed  volume  that  enters  into 
the  expression  for  the  brittle  fracture  stress  (Eq.  [1]).  Knott 
and  Curry19  select  the  “largest  observed"  carbide  thickness 
for  the  value  of  cn  in  Figure  6.  They  also  select  the  95th 
percentile  diameter  to  analyze  results  for  spheroidized  car¬ 
bides  because  this  diameter  leads  to  a  A'£-value  identical  to 
the  one  derived  from  the  "largest  observed"  film.  Other 
workers  have  used  the  average  particle  diameter,  and  while 
this  leads  to  an  incorrect  estimate  of  A'* .  the  ard~ 1  ;  relation 
serves  as  a  useful  signature. 

Brittle  fracture  stress-particle  size  relations  have  been 
reported  for  a  number  of  microstructural  elements  and  are 
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Fig  7  —  Barriers.  (Bl.  lo  microcrack  extension  and  corresponding  dimensions  < u i  grain  boundan  carbide  and  film  thickness  t  .  ibi  femic  gram  boundary 
and  grain  diameter  d.  and  tel  and  id  I  pearlite  lamella  or  colons  boundan  10  and  idi  are  photographs  of  the  surfaces  of  tensile  bars  of  mild  steei  •  steel  L' 
tested  at  - 140  °C  " 


summarized  in  the  convenient  <j-a~' 2  form  of  Eq  [la]  in 
Figure  8.  The  values  of  Kfc  can  be  deduced  from  the  slopes 
of  the  lines,  and  these  show  that  different  boundaries  hav  e 
characteristic  values  of  Kfa .  Fracture  and  the  identity  of  the 
controlling  microstructural  feature  are  determined  by  the 
particle  with  the  (K*  ■  d'  ^-combination  that  produces  the 
largest  value  of  fracture  stress.  The  presentation  in  Figure  6 
facilitates  extrapolation  and  may  be  useful  for  identify  ing 
the  controlling  particle  when  the  microstructural  dimensions 
are  known.  For  example,  microcrack  studies  reported  by 
Hahn  et  al."  included  work  on  a  coarse  grained  ferrite: 
d  =  40  fim.  with  coarse  carbide  films:  c0  =  15  m™  Fig¬ 
ure  8  shows  that  the  fracture  stress  associated  with  micro¬ 
cracks  in  the  ferrite  grains  of  this  size  is  somewhat  larger 
than  the  value  for  microcracks  in  the  carbide  This  accounts 
for  the  observation  of  both  ferrite  and  carbide  microcracks 
in  this  steel.  Ritchie,  Server,  and  Wullaert'*  report  a  brittle 
fracture  stress  a  =  1800  MPa  for  a  heat  of  grade  A533B 
steel  with  bainite  packets:  d  =  10  gun  to  20  gtm.  cementite 
particles:  d  -  0.5  /xm  and  oxide  particle  inclusions,  d  = 
5  Mm  to  10  Mm.  Figure  8  indicates  that  the  brittle  fracture 
stress  associated  with  spherical  carbide  particles  of  this  size. 
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cr  =  4400  MPa.  is  much  larger  than  the  reported  fracture 
stress.  It  is  therefore  unlikely  that  carbide  particles  are  par¬ 
ticipating  in  the  nucleation  process  On  the  other  hand, 
microcracks  nucleated  in  the  oxide  particles  would  extend  at 
a  stress  that  is  lower  than  the  reported  cr-values.*  and  the 

•Presided  Ki.  for  the  oxide  is  close  to  the  saluc  for  the  carbiJe  in 
Figure  K 

fr  icture  stress  of  the  bainite  packets  is  close  to  the  value  for 
the  steel.  The  different  fracture  stress  values  are  consistent 
with  the  view  that  the  mechanism  of  brittle  fracture  of 
A533B  involves  nucleation  of  cracks  in  oxide  inclusions  but 
is  controlled  by  the  bainite  packet  size  The  interpretation 
provides  indirect  evidence  that  hard  particles  other  than  the 
carbide  can  initiate  cracks  in  steel 


VI.  THE  BRITTLE  FRACTURE  STRESS  AND 
THE  TRANSITION  TEMPERATURE 

The  views  expressed  so  far  imply  that  a.  the  brittle  frac¬ 
ture  stress,  is  a  temperature-  .  loading  rate-  .  and  specimen 
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Pig  8  —  Summary  of  brittle  fracture  stress-particle  dimension  measurements  for  different  microstructures, 
showing  the  charactenstic  values  of  A’f  Rets  v.||.  |v.  24.  26-28 


geometry-independent  material  property  because  it  is  a 
function  of  two  quantities,  e  g  .  A'f„  and  t  .  that  are  both 
only  functions  of  the  microstructure  The  transition  tem¬ 
perature,  TT,  is  the  temperature  below  which  the  peak 
normal  stress  exceeds  a,.  and  is  defined  by  the  condition: 

O’.  =  Pc T  =  a  I  31 

where  P  is  the  plastic  constraint  factor  and  cr,  is  the  yield 

The  quantity  P  depends  on  the  mean  stress  and  strain  hardening. 
P  *  1.3.  =2  7.  *3.5.  for  the  tensile-  .  Charps-  .  and  A’ ,  -specimen  tor 
typical  strain  hardening  rates 

strength.  Accordingly,  the  various  transition  temperatures: 
TT  (tensile  bar).  TT  (CVN),  7T  (A\).  and  TT  can  be 
estimated  when  cr  and  a„  at  the  appropriate  temperature 
and  strain  rates  are  known.  This  is  illustrated  in  Figure  9. 
which  shows  effects  of  microstructure  and  the  test  method 
on  the  TT  of  a  hypothetical  medium  strength  steel  that  are 
qualitatively  in  accord  with  experiments.7*"*" 

However,  the  construction  in  Figure  9(b)  that  assumes 
oy  is  independent  of  the  test  piece  configuration  is  question¬ 
able.  Estimates  of  a,  derived  from  tensile-  .  CVN-  ,  and 
K\{ -specimens  of  the  same  material,  using  Eq.  1 3]  and  the 
appropriate  values  of  P  and  cr,,,  just  below  the  TT.  sum¬ 
marized  in  Table  II.  increase  substantially  as  the  stressed 
volume  of  the  test  piece  decreases.  This  is  consistent  with 
the  view,  developed  in  the  preceding  section,  that  the  frac¬ 
ture  stress  is  determined  by  the  size  of  the  largest  "eligible" 


particle  in  the  stressed  volume,  and  the  fact  that  this  size 
increases  with  the  stressed  volume  for  a  materia!  with  the 
usual  histogram  of  particle  sizes  It  seems  likely  that  the 
distribution  of  particle  sizes  depends  on  composition  and 
processing  history,  and  it  follows  from  this  that  the  effect' 
of  changing  the  stressed  volume  vary  from  one  material  to 
another  Some  evidence  of  this  can  be  found  in  Table  11 


VII.  THE  RELATION  BETWEEN 
FRACTURE  TOUGHNESS  AND  THE 
PARTICLE  SIZE  DISTRIBUTION 

The  size  of  the  stressed  volume  enters  in  the  analyri' 
of  brittle  crack  extension  and  A,,  because  it  is  the  extent  ot 


Table  II.  Influence  of  Test  Specimen  Geometry 
and  the  Relative  Size  of  the  Stressed 
Volume  on  the  Brittle  Fracture  Stress 


Relative  Stressed  Volume 

(7.. 

MPa 

Tensile  Bar 
—400 

CVN 

-20 

A, 

1 

Steel  K.  Kotilainen4" 

995 

— 

255U 

Steel  M.  Kotilainen4' 

1505 

— 

3650 

Steel  T.  Kotilainen4" 

1230 

2500 

3015 

Steel  E.  Hahn7" 

550 

810 
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Fig.  9  —  Schematic  representation  of  transition  temperatures  of  a  hypothetical  medium  strength  steel:  (a)  influence  of  the  microstructural  dimension  d  on 
the  TT(CVN)  and  (fc)  influence  of  test  piece  geometry  on  the  TT  for  a  constant  value  of  a,.  Note;  in  practice  o,  varies  with  test  piece  geometry. 


the  highly  stressed  crack  tip  region,  rather  than  the  magni¬ 
tude  of  the  stresses,  that  is  affected  by  A",  and  increases 
with  the  applied  load.  Ritchie.  Knott,  and  Rice'*  therefore 
postulated  the  additional  condition  for  cleavage  crack  ex¬ 
tension.  namely:  Eq.  [3]  must  be  satisfied  a  characteristic 
distance  X(>  in  front  of  the  crack  (see  Figure  10).  A  similar 
view  of  brittle  intergranular  crack  extension  has  been 
presented  by  Jokl,  Kameda.  McMahon,  and  Vitek.”  The 
value  of  A1(  below  the  TT  (A,,  1  is  then  fully  determined  by : 
(1)  the  character  of  the  stress  field.  (2)  the  values  of  ay.  the 
rate  of  strain  hardening,  and  (3)  a,  and  A'„.  Figure  11(a) 
illustrates  that  a  constant  value  of  A,  =  3 d.  where  d  is  the 
average  grain  diameter,  provides  a  rough  description  of 
the  variation  of  Atl  for  A533B  steel  below  TT  (A|, ). 

At  low  temperature,  the  AVrequirement  reduces  to  the 
following  expression  for  toughness: 

A,.  =  <7„\  25AV  [4J 

implying  that  the  particle  spacing  rather  than  its  size  is  the 
controlling  dimension.  The  correlation  between  X(,  and  Ar. 
the  mean  free  path  between  carbide  particles,  illustrated  by 
Rawal  and  Gurland’s71  measurements  in  Table  III.  supports 
this  concept.  While  Ai.-values  below  the  7T(A'i.)  are  dom  - 
nated  by  a  particle  spacing  (Eq.  (4]),  CVN- values  may  be 
more  strongly  influenced  by  the  largest  eligible  particle 
size  (Eq.  [1]).  This  notion  and  the  fact  that  the  spacing 
and  size  can  vary  independently  account  for  the  absence  of 
a  single,  well-defined  relation  between  the  two  standard 
measures  of  toughness  valid  for  all  steels. ,k  A  classical  ex¬ 
ample  of  non-correlation,  reported  by  Ritchie.  Francis,  and 
Server.73  is  illustrated  in  Table  IV. 

The  characteristic  distance  was  initially  connected  with 
the  need  to  accommodate  a  dislocation  pile-up  and  an 
“eligible"  carbide  in  the  highly  stressed  region  and  esti¬ 
mated  to  be  a  small  multiple  of  the  average  grain  diameter. 
d  (see  Figure  11(a)).  Subsequent  studies  that  reveal  A',  is 


independent  of  d  (see  Figure  11(b))  have  led  Curry  and 
Knott54  to  the  view  that  X0  is  connected  with  the  volume  of 
material  needed  to  assure  the  presence  of  an  "eligible"  par¬ 
ticle  whose  fracture  condition  is  satisfied.  Curry  and  Knott3' 
have  formulated  this  concept  statistically  and  use  it  to 
analyze  the  fracture  toughness  of  quenched  and  tempered 
microstructures  with  spheroidized  carbides.  Their  analysis 
recognizes  contributions  from  all  the  "eligible"  panicles, 
large  and  small  (see  Figure  10(d)).  It  involves  three  statis¬ 
tical  quantities:  (i)  P(r,),  the  probability  a  panicle  has  a 
radius  r,.  (ii)  P(f.r ,).  the  probability  a  particle  of  radius  r 
will  initiate  fracture,  and  (iii)  8,  the  probability  a  particle  is 
"eligible",  and  the  condition  for  crack  extension  * 

* Pi  f.r ,)  =  BsX:naPu,).  where  X  is  Ihe  distance  from  the  cra.k  up  tor 
which  ihe  fraclure  condtion  (Eqs  |1)  and  |?]i  is  satisfied  tor  t  -parnJes. 
5  is  a  factor  of  order  unity  descnbinp  the  shape  of  ihe  stressed  region. 
sX:  is  area  wilhin  which  the  fracture  condition  is  satisfied,  and  «  is  the 
loia)  number  of  particles  per  unit  area  Equation  |?)  defines  A',  =  A  where 
ZPl/M  =  1  since  X,  =/(r,.A'|) 


2>(/.r,)  =  1  151 

i 

Results  of  the  analysis,  in  Figure  12.  reveal  that  the  -5  pci 
of  the  particle  population  with  the  largest  radii  a-e  the 
most  probable  source  of  fracture  in  the  crack  tip  region 
This  arises  because  the  reduced  stresses  needed  to  fracture 
the  largest  particles  are  obtained  in  a  volume  ~  30  ‘  larger 
than  that  corresponding  to  average-size  particles  The 
analysis  indicates  that  the  contribution  of  the  different  par¬ 
ticle  sizes  vanes  with  the  temperature,  and  that  the  fraction 
of  the  eligible  particles  10:  s  8  s  10'  is  roughly  indepen¬ 
dent  of  temperature.  It  should  be  noted  that  the  value  of  8 
depends  on  the  nucleation  process  and  enters  the  statis¬ 
tical  fracture  condition.  For  example,  the  Curry  and  Knott 
analysis  indicates  that  a  15-fold  decrease  of  8  leads  to  a 
2-fold  increase  of  A|. .  In  other  words,  the  microstructural 
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Fig  10  —  Mechanisms  of  brittle  crack  extension  <al  and  (hi  transgranular  crack  extension  after  Riuhic.  Knott,  and  Rice'  .  u  umergranular  crack  euensn  i 
after  Jokl.  Kameda.  McMahon,  and  Vitek.'  and  id i  statistical  interpretation  of  Curts  and  Knott  s 


Table  III.  Relation  between  the  Characteristic  Distance  (A0>  and  the  Carbide  Particle  Mean  Free  Path  ( 
in  Spheroidized  Carbon  Steels  Tested  at  -196  °C  According  to  Measurements  bv  Rawal  and  C.urland’1 


Steel 

Carbon  Weight.  Pet 

a,..  MPa 

A'i. .  MPa  Vm 

a.  X.i.  jam 

b.  A, .  jim 

RC 

0.13 

766 

18.0 

22.1 

39  3 

YC 

0.43 

912 

14  9 

10  7 

8.05 

BC 

0.83 

841 

13.7 

10.6 

7  88 

GC 

1  46 

813 

113 

7.7 

5.28 

YF 

0  43 

997 

15  6 

9.8 

4  32 

BF 

0.83 

946 

14.6 

9.5 

4  86 

a.  X„  =  25  of  i 

b  Ar  =  ( I  -  f)iS ir.  where  /  is  the  volume  fraction  calculated  from  the  chemical  analyses  and  ,V,r  is  the  number  of  panicles  intercepted  per  unit  lengt’ 
of  random  test  line' 


features  that  control  nucleation  can  make  a  modest  con¬ 
tribution  to  K\, . 

The  Curry  and  Knott  analysis  is  based  on  the  questionable 
assumption  that  the  fracture  condition,  Eq.  [5],  must  be 
satisfied  on  every  cross  section  because  "the  initiation  of 
fracture  has  to  be  a  more  or  less  sympathetic  event  across  a 
large  proportion  of  the  specimen  thickness  to  allow  suf¬ 
ficient  elastic  energy  release  In  contrast,  if  the  fracture 
condition  must  be  satisfied  only  once  in  the  stressed  volume 
rather  than  on  every  cross  section.*  the  relative  contribution 

*Pt/.  M  -  8BSx:  lna/2r,i  Pir, i  where  B  is  the  thickness 


Table  IV.  Example  of  Non-Correlating  CVN- 
and  AVToughness  Measurements  of  AISI  4340 
Steel  after  Ritchie,  Francis,  and  Server72 


Condition* 

CVN-Energy.  J 

At  .  _  A,.,.  _ 

MPa  \  m  MPa  \  m 

A 

6  6 

70  52 

B 

9  5 

34  35 

*A — 1  hour-1200  C.  sail  quench  to  K'u  C.  1  hour-S'ii  C.  oil  iiuer.h 

to  Rt.  a  =  159. 

2  MPa 

B  —  1  hour-8' 

?0  CC.  oil  quench  to  RT  a  =  15*0  MPj 
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Fig  II  — Estimates  of  the  characteristic  distance  A'  to)  after  Ritchie.  Server,  and  Wullaen*  and  (b)  after  Cum  and  Knott 
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Fig  12  —  Contribution  of  the  carbide  panicle  size  distribution  to  the  fracture  toughness  of  a  quenched  and  tempered  steel  with  spheroidized  carbide  particle 
deduced  by  Cum  and  Knott-'  with  an  areal  analysis,  and  the  corresponding  values  for  a  volumetnc  analysis  The  quantity  Ptr, )  is  the  probability  a  panivlt 
has  a  radius  r,\  P(  fr, i  is  the  probability  a  particle  of  radius  r.  will  initiate  fracture 
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of  the  different  particle  sizes  is  not  altered  very  much 
(see  Figure  12).  but  the  estimated  fraction  of  "eligible" 
particles  is  reduced  to  0*  KT".  This  fraction  is  consistent 
with  the  large  number,  N  -  1.1*  lO”  mm‘\  of  carbide  par¬ 
ticles  in  steel,  but  it  is  so  small  as  to  raise  questions  about 
the  identity  of  the  nucleating  particle.  The  general  con¬ 
clusion  to  be  drawn  from  either  the  area!  or  volumetric 
assumption  is  that  the  very  small  number  of  the  largest 
particles  has  a  disproportionate  influence  on  both  TT(K\C) 
and  the  value  of  Ku  below  the  transition  temperature. 


the  different  transition  temperatures  and  can  make  the 
A'k-CVN  energy  relation  microstructure  dependent 
9.  Statistical  analyses  of  the  contribution  of  different  par¬ 
ticle  sizes  indicate  that  the  small  number  of  the  very 
largest  particles  has  a  disproportionate  influence  on  K, 
below  the  transition  temperature.  Microstructura)  fea¬ 
tures  that  influence  the  fraction  of  the  particles  that  is 
“eligible"  have  a  modest  effect  on  Ku . 
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VIII.  CONCLUSIONS 

1 .  The  brittle  fracture  of  steel  is  governed  by  the  extension 
of  fast  propagating  Griffith-type  microcracks  across 
microstructural  boundaries. 

2.  The  microcracks  are  not  preexisting  defects  and  must 
first  nucleate.  Carbides  and  possibly  other  brittle  par¬ 
ticles  have  an  essential,  intermediate  role  because,  unlike 
the  ferrite,  they  permit  the  crack  nucleus  to  form  without 
blunting  it. 

3.  The  brittle  fracture  stress  depends  on  the  values  of  two 
quantities:  (i)  the  size  of  the  microcracks;  and  (ii)  the 
resistance  it  experiences,  relevant  when  the  microcrack 
attempts  to  cross  the  boundary.  The  resistance  is  ex¬ 
pressed  by  the  Griffith  energy  term:  15  JnT:  ^  y  ^ 
100  Jm“\  or  a  corresponding  LEFM  crack  arrest  tough¬ 
ness  value:  25  MPa  Vm  <  Ku  s  7  MPa  Vm,  that 
applies  to  a  rapidly  advancing  microcrack.  The  size  of 
the  microcrack  corresponds  with  the  parent  particle's 
dimension. 

4.  Since  the  parent  particle  or  grain  is  represented  in  the 
microstructure  by  a  distribution  of  sizes,  the  brittle  frac¬ 
ture  stress  is  determined  by  the  largest  “eligible"  particle 
within  the  stressed  volume.  Particles  are  "eligible”  when 
the  conditions  for  nucleating  a  crack  can  be  satisfied  and 
when  the  resistance  their  boundaries  offer  to  microcrack 
extension  is  not  excessive. 

5.  The  crack  nucleation  in  hard  particles  is  assisted  by  plas¬ 
tic  deformation  of  the  surrounding  matrix  but  requires  an 
additional  stress  raiser  such  as  a  dislocation  pile-up  or  a 
defect  in  the  particle.  However,  the  pile-up  dimensions 
do  not  influence  the  brinle  fracture  stress  directly.  The 
microstructural  features  that  affect  nucleation  can  alter 
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the  brittle  fracture  stress  indirectly  by  changing  the 
number  of  “eligible”  particles  and  the  largest  “eligible" 
particle  size. 

6.  The  particles  whose  boundaries  resist  microcrack  exten¬ 
sion  include:  (i)  carbide  particles,  (ii)  carbide  films  at  the 
grain  boundaries,  (iii)  carbide  lamellae  in  the  pearlite, 
(iv)  hard  inclusions,  (v)  ferrite  grains,  (vi)  pearlite  colo¬ 
nies  or  grains,  (vii)  bainite  packets,  and  (viii)  martensite 
packets.  The  different  boundaries  display  characteristic 
/(£- values. 

7.  The  microstructura]  feature  that  controls  the  brittle 
fracture  of  a  particular  steel  is  the  feature  with  the 
“eligible"-particle-size/AT®  combination  that  produces 
the  largest  brittle  fracture  stress  value.  The  existing 
/(^-measurements  provide  insights  to  the  controlling 
microstructural  features  in  different  steels. 

8.  The  differences  in  the  size  of  the  stressed  volume  in  the 
tensile-,  CVN-,  and  AT^-test  complicate  the  analysis  of 
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FOREWORD 


Establishment  of  Microstructural  Spacing 
during  Dendritic  and  Cooperative  Growth 


This  symposium  was  organized  under  joint  sponsorship  by  the  ASM-MSD  Phase  Transformations 
Committee  and  the  TMS-AIME  Solidification  Committee  to  provide  a  forum  for  discussion  of  micro- 
structural  spacing  in  different  morphologies  during  various  transformations.  The  invited  contributions 
were  presented  during  the  March  1983  meeting  of  TMS-AIME. 

The  size  scale  and  phase  spacing  are  among  the  most  important  characteristic  features  of  regular 
microstructural  morphologies.  In  practice  a  number  of  thermal  and  mechanical  treatments  have  been 
developed  to  modify  microstructural  spacing  during  phase  transformations.  For  some  time  there  has  been 
no  generally  established  analysis  to  account  for  the  selection  of  spacing  under  given  conditions.  However, 
recent  advances  in  experiment  have  yielded  accurate  measurements  in  well-characterized  systems  which 
have  stimulated  new  theoretical  work  on  the  analysis  of  the  selection  of  spacing  and  scale. 


Important  theoretical  progress  has  been  made  in  applying  the  principle  of  marginal  stability  to  the  initial 
selection  of  dendrite  spacing  and  with  the  description  of  dendrite  coarsening  to  account  for  structural 
adjustment.  Further  progress  in  analysis  and  new  results  have  been  reported  on  spacings  developed  in 
cellular,  eutectic,  and  monotectic  structures.  For  solid  state  reactions  recent  experimental  work  on  spacing 
selection  has  emphasized  the  importance  of  diffusional  transport  in  the  vicinity  of  the  moving  trans¬ 
formation  interface.  The  written  discussions  in  the  papers  that  were  available  at  the  time  of  publication 
provide  a  valuable  overview  of  recent  work  and  point  out  a  number  of  unresolved  areas  which  should 
stimulate  future  work. 
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REFERENCE:  Hoff,  R.,  Rubin,  C.  A.,  and  Hahn,  G.  T.,  “Strain-Rate  Dependence  of 
the  Deformation  at  the  Tip  of  a  Stationary  Crack,"  Fracture  Mechanics:  Sixteenth  Sym¬ 
posium,  ASTM  STP  868,  M.  F.  Kanninen  and  A.  T.  Hopper,  Eds.,  American  Society 
for  Testing  and  Materials,  Philadelphia,  198S,  pp.  409-430. 

ABSTRACT:  A  new  analytical  technique  is  presented  whereby  the  Hutchinson-Rice- 
Rosengren  singularity  solutions  can  be  used,  in  certain  cases,  to  determine  the  stresses 
and  strains  near  a  crack  tip  for  a  dynamically  loaded,  strain  nte-dependent  material. 

This  technique  involves  a  simple  modification  of  the  constant,  a,  and  the  strain  hard¬ 
ening  exponent,  n.  in  the  power  stress-strain  law.  Elastoplastic,  dynamic  finite  element 
computations  are  in  close  agreement  with  this  analytical  technique. 

Stresses  and  strains  are  generally  higher  in  a  dynamically  loaded  rate  sensitive  mate¬ 
rial  than  in  a  rate  insensitive  material.  Also,7ic  increases  correspond  to  increases  in  the 
flow  stress  in  rate-sensitive  materials,  provided  that  6c  is  rate-insensitive. 

KEY  WORDS:  fracture,  cracks,  plastic  fracture  mechanics,  dynamic,  strain  rate, 
7-integral,  crack  opening  displacement,  finite  element  models,  compact  tension  speci¬ 
men,  steels,  power-law  hardening  materials 

Nomenclature 

Ap  Crack  tip  plastic  zone  area 

dn  Constant  in  crack  tip  opening  displacement  relation  (Eq  2). 

Value  depends  on  the  strain  hardening  exponent  n 
D  Constant  in  Malvern  equation  (Eq  14) 

E  Young’s  modulus 

In  Constant  in  /-field  equations  (Eqs  S  and  6).  Value  depends  on 
the  strain  hardening  exponent,  n 

/,  J\  Path-independent  contour  integral;  energy  release  rate  for  quasi- 
statically  loaded  specimens 
/,  /i  Time  rate  of  change  of  the  /-integral 

1  Graduate  Assistant,  Associate  Professor,  and  Professor,  respectively.  Department  of  Me¬ 
chanical  and  Materials  Engineering,  Vanderbilt  University,  Nashville,  TN  37233.  Dr.  Hoff  is 
currently  Assistant  Professor,  Department  of  Mechanical  Engineering,  University  of  Waterloo, 
Waterloo,  Ontario,  Canada. 
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<v 

«v 

<v 


€P 

<p 

<P(ava) 

<P(av«) 


Path-independent  integral;  energy  release  rate  for  dynamically 
loaded  specimens 

Critical  value  of  7-integral  when  crack  propagation  initiates 
7ic  at  loading  rates  a  and  b,  respectively 
Maximum  value  of  J\ 

Maximum  value  of  J\ 

Value  of  /-integral  at  time  =  0 
Coordinate  directions:  Xk  —  X\ ,  Xj 

Mode  I  stress  intensity  factor  (loosely  defined  as  \jEJ/\  —  v2) 

Time  rate  of  change  of  the  stress  intensity  factor 

K ic  at  loading  rates  a  and  b,  respectively 

Strain  hardening  exponent  in  the  power  stress-strain  law  (Eq  3) 

Strain  hardening  exponent  in  the  modified  power  stress-strain  law 

Direction  cosines 

Reciprocal  of  the  strain  hardening  exponent;  N  =  1  /n 
Reciprocal  of  the  modified  strain  hardening  exponent;  N'  =  1 /n‘ 
Exponent  in  the  Malvern  equation  (Eq  14) 

Radial  distance  from  the  crack  tip  (length  units) 

Plastic  zone  size  at  6  =  0 

Dimensionless  distance  from  crack  tip;  R  =  r/(K\/oof 
Cauchy  deviatoric  stress  tensor 
Arc  length 
Time 

Traction  vector;  /,  =  ayrij 

Kinetic  energy  density;  T  —  (1/2)  pit, it, 

Displacement  vector 
Velocity  vector 
Acceleration  vector 
Volume  f  • 

Strain  energy  density;  W  =  Jn  a, fit,, 

Constant  in  the  power  stress-strain  law  (Eq  3) 

Constant  in  the  modified  power  stress-strain  law  (Eq  18) 

Integration  contour 

Crack  tip  opening  displacement 

Critical  crack  tip  opening  displacement 

Strain  tensor 

Time  rate  of  change  of  the  strain  tensor 

Known  dimensionless  functions  of  the  circumferential 

position,  6,  and  the  strain  hardening  exponent,  n 

Equivalent  plastic  strain 

Equivalent  plastic  strain  rate 

Average  equivalent  plastic  strain 

Average  equivalent  plastic  strain  rate 
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<o  Yield  strain  under  static  tensile  loading;  e0  =  oo/E 
6  Circumferential  position  measured  counterclockwise  from  the 
crack  plane 

X  Loading  rate  parameter 
v  Poisson’s  ratio 
p  Mass  density 

od  Instantaneous  flow  stress  under  dynamic  loading 
ot  (von  Mises)  equivalent  stress 
av  Cauchy  stress  tensor 

b,j  Known  dimensionless  functions  of  the  circumferential  position,  8, 
and  the  strain  hardening  exponent,  n 
oo  Yield  stress  under  static  tensile  loading 
ooa,  oob  Yield  stress  at  loading  rates  a  and  b ,  respectively 
a.  Flow  stress  under  static  loading 


As  a  first  step  in  modelling  rapidly  advancing  cracks  in  tough  steels,  this 
paper  examines  a  stationary  crack  in  a  rapidly  loaded  compact  tension  spec¬ 
imen.  The  plastic  strain  rates  and  therefore  the  strain-rate  sensitivity  effects 
are  smaller  for  a  dynamically  loaded  stationary  crack  than  for  a  rapidly 
propagating  crack.  However,  strain-rate  effects  are  by  no  means  insignifi¬ 
cant  in  the  stationary  crack  case  [/].  Yet,  with  the  exception  of  the  recent  fi¬ 
nite  difference  study  of  elasto-viscoplastic  fracture  by  Aboudi  and  Achen- 
bach  [2]  and  creep  studies  by  Little  et  al  [5],  virtually  no  analytical  or 
numerical  work  has  been  reported  incorporating  strain-rate-dependent  flow 
properties  in  the  determination  of  the  crack  tip  stress  fields.  The  aim  of  this 
work  is  to  examine  the  effects  of  loading  rate  and  rate  sensitivity  on  J\  and 
the  crack  tip  opening  displacement,  features  of  the  plastic  and  process  zones 
that  are  influential  in  determining  the  toughness. 

Plastic  strains,  in  the  region  surrounding  a  stationary  crack,  can  be  char¬ 
acterized  as  occurring  in  two  zones  (Fig.  1),  Small  plastic  strains  in  the  range 
0<«,<0.1  occur  in  the  larger  region  called  the  plastic  zone.  The  size  of  the 
plastic  zone  (at  6  =  0)  is  given  by  a  characteristic  dimension,  ro.  In  a  non¬ 
hardening  material  where  the  plastic  zone  is  small  compared  with  the  speci¬ 
men  dimensions  [4]: 


r0 


0.036  EJ 
a0:  (1  -  v2j 


(1) 


In  the  immediate  vicinity  of  the  crack  tip  is  the  heavily  strained  fracture  proc¬ 
ess  zone,  where  plastic  strains  range  roughly  from  0.1  to  1.0  (/].  The  size  of 
this  intensely  nonlinear  zone  is  on  the  order  of  the  crack  tip  opening  dis- 
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FIG.  1 — Plastic  and  process  zones  in  the  vicinity  of  a  stationary  crack. 


placement,  6.  The  crack  tip  opening  displacement  can  be  expressed  as  a  func¬ 
tion  of  Rice’s  [J]  contour  integral  J  by  [ 6] 


S  —  dn  — 
o0 


(2) 


for  a  nonlinear  elastic,  power  law  hardening  material.  The  constant,  d„,  is  a 
function  of  the  material  properties  a,  oo,  E,  and  n  (in  the  power  stress-strain 
law),  where 


-  =  “(-) 
Co  '  Oo' 


(3) 


The  relative  size  of  the  plastic  and  process  zones  is  approximately 


—  =*  0.03 
To 


(4) 


for  a  material  such  as  A533B  [7].  Consequently  the  volume  of  the  process 
zone  is  only  0.09%  of  the  plastic  zone  volume,  in  two-dimensional,  small-scale 
yielding,  plane-strain  situations. 

Expressions  for  the  stresses  and  strains  in  the  plastic  zone  are  given  by  the 
/-field  solutions  of  Hutchinson  [7]  and  Rice  and  Rosengren  [£]: 


on 


(pj  \ 

- fr-l  b,j{d,  n) 

aoo  Ur) 


(5) 
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and 

(£j  \  [n/n*  1) 

- rr~ )  cv(8,  n)  (6) 

a  oo  hrj 

Also  of  interest  are  expressions  for  the  equivalent  stress  and  equivalent  plas¬ 
tic  strain: 


ot  =  \J (3/2)  SijSij 


and 

ep  =  \J (2/3)  tijty 


0) 

(8) 


Since  the  power  law  relation  (Eq  3)  does  not  accurately  represent  the  linear 
elastic  portion  of  the  stress-strain  curve,  Eqs  5  to  8  are  most  accurate  where 
elastic  strains  are  negligible. 

Strain  rates  can  be  determined  by  differentiating  Eq  6  with  respect  to  time: 


where  the  dot  superscript  denotes  a  time  derivative.  Observing  the  similarity 
between  Eqs  6  and  9,  the  latter  can  be  simplified  to  give 


The  strain  rates  given  in  Eqs  9  and  10  are  not  accurate  in  most  situations. 
The  angular  distributions  of  strain  rates,  represented  by  ev  in  Eq  9,  were  de¬ 
termined  for  static  loading  and  therefore  do  not  account  for  the  effect  of  the 
inertial  loads  on  the  stresses.  The  amplitude  of  the  strain  rate  singularity  is  a 
function  of  a,  n,  a0,  /* ,  /,  and  /.  Many  materials  exhibit  strain-rate-dependent 
material  properties,  in  which  case  the  static  oo,  a,  and  n  values  do  not  repre¬ 
sent  the  stress-strain  relations  under  dynamic  loading. 

The  application  of  Rice’s  /-integral  to  general  dynamic  situations  is  awk¬ 
ward,  since  it  does  not  account  for  kinetic  energy  or  inertia  and  is  no  longer 
path-independent.  Recently,  a  number  of  path-independent  integrals  have 
been  proposed  [9-12]  that  are  applicable  to  general  dynamically  loaded 
structures.  The  most  promising  of  these  is  the  /’-integral  of  Atluri  et  al  [12], 
which  corresponds  physically  to  energy  release  rate  for  stationary  and  prop¬ 
agating  dynamically  loaded  cracks.  The  /’-integral  has  some  computational 
advantages  over  the  others,  since  the  integration  is  performed  along  a  fixed 
contour  in  space  and  not  along  a  contour  that  moves  with  a  propagating 
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crack  such  as  in  Ref  10.  The  T  vector  integral  is  given  by 

Jk~  i  +  f  [pwfc*  -  pu,u,k]dV  (11) 

whose  first  component,  /j,  degenerates  to  Rice’s  /-integral  for  stationary, 
statically  loaded  cracks.  The  domain  of  the  second  integral  in  Eq  1 1  is  the  vol¬ 
ume  contained  within  the  contour  of  the  first  integral  in  the  equation. 

Expressions  such  as  Eqs  5  and  6  for  stresses  and  strains  near  a  crack  tip  are 
not  available  for  most  dynamically  loaded  structures.  However,  Eqs  5  and  6 
can  still  be  used  for  slowly  loaded  specimens  and  for  rapidly  loaded  speci¬ 
mens  with  stationary  cracks  by  using  the  “inertialess”  modelling  technique 
discussed  in  the  next  section. 

Ductile  crack  extension  proceeds  by  void  nucleation,  growth,  and  coales¬ 
cence  in  the  process  zone.  The  critical  crack  tip  opening  displacement,  8C,  has 
been  demonstrated  to  be  a  viable  ductile  fracture  criterion  [13, 14].  From  Eq  2 
one  can  deduce  Jic  in  terms  of  the  critical  crack  opening  displacement: 


Ju  — 


Oo8c 

17 


(12) 


In  rate-sensitive  materials,  the  yield  stress  will  generally  increase  with  in¬ 
creasing  strain  rates.  If  6C  is  rate-insensitive,  then  changes  in  oo  will  alter  Ju 
according  to 

/left  /  K\cb 

/ic a  \  K\ca 

where  the  subscripts  a  and  b  in  Eq  13  refer  to  two  different  loading  rates.  In 
other  words,  J\c  will  increase  with  loading  rate,  provided  that  6C  is 
rate-insensitive. 

The  void  nucleation  and  growth  mechanism,  which  initiates  the  ductile 
fracture  process,  is  enhanced  by  higher  mean  (or  hydrostatic)  stress.  The 
mean  stress  scales  with  static  yield  stress  (see,  for  instance,  the  slip-line  file 
solution  of  Ref  5);  however,  6C  varies  inversely  with  mean  stress  [/5].  There¬ 
fore  a  higher  yield  stress  will  generally  produce  a  lower  critical  crack  tip  open¬ 
ing  displacement,  which  in  turn  would  decrease  Ju.  Clearly,  these  two  com¬ 
peting  effects  are  occurring  simultaneously,  and  it  is  not  known  which  effect 
will  dominate. 

In  the  present  study,  we  consider  only  the  former  case,  where  Ju  increases 
because  of  the  flow  stress  elevation  in  the  dynamically  loaded  plastic  zone. 
An  analytical  model  is  developed  whereby  the  /-field  solutions  are  modified 
to  accommodate  strain  rate  sensitivity.  Finite  element  computations  are  also 
performed  to  verify  the  accuracy  of  this  technique. 


)2  _  OOb 
OOa 


(13) 


HOFF  ET  AL  ON  DEFORMATION  AT  TIP  OF  STATIONARY  CRACK  415 

Experimental  data  from  a  number  of  sources  [16-20]  related  J ic  to  the 
loading  rate,  /  Materials  that  exhibit  an  increase  in  yield  stress  at  high  strain 
rates  often  show  an  increase  in  Ju  at  high  strain  rates  corresponding  to  Eq  13, 
as  was  observed  in  a  previous  paper  [/].  The  relative  contributions  of  the 
strain  rates  in  the  plastic  and  process  zones  are  still  unclear. 

Estimation  Technique  for  Strain-Rate  Effects 

The  stresses  and  strains  near  the  tip  of  a  dynamically  loaded  crack  in  a 
ductile  material  are  affected  by  the  inertia  of  the  structure  and  by  the  strain- 
rate  sensitivity  of  the  yield  or  flow  stress.  It  is  therefore  desirable  to  isolate 
the  strain-rate  sensitivity  effects  in  order  to  ascertain  their  significance.  This 
section  describes  a  technique  whereby  the  parameters  in  the  power  stress- 
strain  law  are  modified  to  account  for  strain-rate  effects. 

Suppose  that  a  body  were  dynamically  loaded,  but  that  inertial  effects 
could  be  ignored.  Physical  situations  that  would  correspond  to  this  are  (l)a 
solid  body  with  negligible  density,  or  (2)  a  body  loaded  at  constant  velocity 
so  that  the  accelerations  are  everywhere  negligible  (perhaps  after  an  initial 
transient).  Differences  between  the  static  behavior  and  the  dynamic  behav¬ 
ior,  in  such  instances  would  be  entirely  due  to  the  strain-rate  effects. 

Experimental  results  have  shown  that  the  flow  stress  typically  increases 
with  increasing  strain  rate  [16-20].  Numerous  models  have  been  proposed  to 
characterize  this  behavior,  some  of  which  are  discussed  in  Ref  21.  We  have 
chosen  a  form  of  the  equation  used  by  Malvern  [22]  because  of  its  computa¬ 
tional  simplicity,  although  any  other  stress-strain-ratc  relation  could  equally 
well  have  been  used.  Malvern  suggested  that  the  instantaneous  flow  stress 
under  dynamic  loading,  od,  is  a  function  of  the  equivalent  plastic  strain  rate, 
ep,  and  the  flow  stress  at  zero  strain  rate,  a„  according  to 

oa  =  o.  [(£)'"+.]  (.4) 

The  constants  o»,  D,  and  p  must  be  determined  experimentally. 

The  power  stress-strain  law  (Eq  3)  can  be  inverted  to  obtain 

N 

(15) 

where  N  —  l/n.  Equation  15  can  be  substituted  into  Eq  14  to  give 

oi=o"(i;)  [(f)  +,1 


(16) 
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Substituting  for  the  equivalent  plastic  strain  rate  from  Eq  10  gives 


-*(±)'[( 


(1  +  N)DJt 


Equation  17  is  too  unwieldy  to  be  used  as  a  stress-strain  relation.  A  more  de¬ 
sirable  dynamic  stress-strain  relation  would  be  a  power  law  of  the  form 


■-  (  <  V 


The  parameters  a'  and  N'  cannot  be  found  by  equating  Eq  18  to  Eq  17  and 
solving  for  a  and  N*.  However,  a  very  good  approximation  to  Eq  17  can  be 
found  by  selecting  a  and  N'  to  minimize  the  deviations  (squared)  between 
Eqs  17  and  18.  That  is,  choose  a  and  N'  to  minimize 


(a,  N')  =  f  (oi~  0d)2  dt 


Or,  making  substitutions  from  above,  to  minimize: 


—/((*)'-(=)■[( 


Jt 

(1  +  N)DJ  < 


+  1 1!  dt  (20) 


(This  minimization  is  best  performed  by  using  a  nonlinear  programming  ap¬ 
proach.  Ref  23  discusses  nonlinear  programming  in  detail.)  The  limits  of  in¬ 
tegration  in  Eq  20  should  be  chosen  to  reflect  the  strain  range  of  interest,  typ¬ 
ically  0  to  0.25.  Note  that  Eq  20  only  applies  to  a  particular  instant  in  time, 
since  J/J  will,  in  general,  not  be  constant  (///  is  constant  only  in  the  special 
case  where  J  varies  with  time  according  to  J  =  Jo  exp{A/}). 

Figure  2  illustrates  the  strain-rate-modified  stress-strain  curve  obtained  by 
using  this  approximation  technique.  The  lower  curve  is  the  zero  strain  rate 
stress-strain  curve  for  a  material  such  as  A533B  ( e/to  =  2.12(a/ao)8  4‘).  The 
upper  curve  is  the  actual  dynamic  stress-strain  curve  produced  according  to 
Eq  17,  with  a.  =  383  MPa,  D  =  62  920/s,  p  =  7.245,  J  =  188.0  kJ/m2,  and 
J  =  176.5  MJ/m2s.  The  curve  reflects  the  variation  of  ep  with  ep  throughout 
the  plastic  zone.  This  curve  can  be  contrasted  with  the  conventional  stress- 
strain  curve,  shown  as  a  dashed  line  in  Fig.  2,  which  assumes  a  constant  av¬ 
erage  strain  rate  throughout  the  plastic  zone.  The  high  strain-rate  estimation 
technique  was  applied  to  parameterize  the  upper  curve.  Values  of  a  =  0.533 
and  n'  =  6.38  correspond  to  the  high  strain  rate  curve  so  well  that  any  devia¬ 
tions  are  not  visible  on  the  scale  of  Fig.  2. 
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FIG.  2 — Power  stress-strain  law  curves  for  A533B  steel  under  static  and  dynamic  loading. 

Finite  Element  Modelling 

A  stationary  crack  in  a  compact  tension  specimen  was  modelled  using  fi¬ 
nite  elements.  Figure  3  is  a  plot  of  the  mesh;  an  enlarged  view  of  the  crack  tip 
region  is  shown  below  a  view  of  the  bulk  of  the  model.  The  model  is  com¬ 
prised  of  646  nodal  points  joined  by  190  eight-noded  isoparametric  finite 
elements.  Numerical  studies  [24,25]  have  shown  these  elements  to  be  suitable 
for  use  in  the  fully  plastic  range. 

The  crack  is  modelled  with  a  sharp  tip;  the  finite  element  analyses  do  not 
account  for  the  continually  changing  tip  geometry  as  progressive  blunting 
takes  place.  This  is,  however,  the  same  way  in  which  the  /-theory  models  the 
crack  tip.  Therefore  the  results  of  the  finite  element  analyses  ought  to  agree 
with  the  /-field  equations  (Eqs  5  to  8).  Another  study,  currently  under  way, 
is  focusing  on  the  more  realistic  problem  where  the  continually  changing 
crack  tip  geometry  is  included  in  the  analysis. 

The  material  properties  are  given  with  an  elastic  strain-hardering  plastic 
model,  based  on  experimental  data  for  A533B  steel  at  93°C  [26,27].  The  zero 
strain  rate  stress-strain  curve  is  given  in  Fig.  4.  Static  and  dynamic  finite 
element  analyses  were  performed  using  this  material  mode.  Dynamic  anal¬ 
yses  were  also  performed  using  strain-rate-sensitive  flow  stress  as  modelled 
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FIG.  3 — Finite  element  model  of  compact  tension  specimen. 


by  Eq  14.  The  three  parameters,  o.  =  383  MPa,  D  =  62920/s,  and p  =  7.245, 
were  determined  by  HofT  ct  al  [28]  for  A533B  steel  at  93°C. 

The  interaction  of  strain  rate  and  inertial  effects  makes  it  difficult  to  iso¬ 
late  the  contribution  of  either.  A  technique  was  devised  to  eliminate  the  iner¬ 
tial  effect  in  order  to  quantitatively  assess  the  importance  of  the  strain  rate 
effect.  A  linear  elastic  static  analysis  was  first  performed.  All  the  resulting 
nodal  displacements  were  scaled  by  a  constant  giving  initial  velocities  at  each 
node.  The  scaling  constant  was  chosen  so  that  the  load  point  velocity  pro¬ 
duced  by  this  scaling  technique  was  equal  to  that  load  point  velocity  required 
to  give  the  desired  A’l-rate,  K\.  A  dynamic  analysis  was  then  performed,  start¬ 
ing  with  an  unloaded  model  (all  displacements  zero)  but  with  an  initial  veloc¬ 
ity  distribution.  The  nodal  velocities  did  not  change  appreciably  during  the 
analysis.  With  almost  constant  nodal  velocities,  the  nodal  accelerations  are 
very  small  and  the  model  is  essentially  “inertialess”.  Further  evidence  of  the 
inertialess  character  of  this  analysis  is  that  the  kinetic  energy  varied  by  only 
±3.4%  throughout  the  analysis. 

Several  finite  element  analyses  were  performed  under  static  and  inertialess 
dynamic  loading.  Table  1  summarizes  the  essential  features  of  these  analyses. 
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FIG.  4—  Zero  strain  rate  stress-strain  curve  for  A533B  steel  at  93°C. 


Results  and  Discussion 

Results  for  the  static  analysis,  S04,  are  given  in  Figs.  5  to  7  at  a  load  level 
of  J\  =  182  kJ/m2.  Figure  5  shows  contours  of  equivalent  stress.  Note  that 
all  distances  are  scaled  by  the  similarity  parameter,  ( Ki/oo)2 ,  so  that 


‘  (Ay  oo)2 

Figure  6  shows  contours  of  equivalent  plastic  strain  in  the  same  region  as 
Fig.  5.  A  comparison  of  the  equivalent  plastic  strain  at  6  =  45°  from  the  fi¬ 
nite  element  computation  with  the  /-field  solution  of  Eq  6  is  shown  in  Fig.  7. 
Note  the  good  agreement  between  these  two  approaches. 

The  average  plastic  strain  can  be  determined  from  an  area  average  of  the 
plastic  strains  in  the  crack  tip  plastic  zone  according  to 


€P(«v«) 


-±/ 


«p  dA\ 


In  spite  of  the  larger  plastic  strains  near  the  crack  tip,  the  average  is  «p<.ve) 
=  0.0016  due  to  the  large  portion  of  the  plastic  zone  that  is  lightly  strained. 

Analysis  S16  features  "fast”  inertialess  dynamic  loading  (see  Table  1  for 
definition  of  fast  and  slow)  without  strain-rate-sensitive  flow  stress.  The 
equivalent  stress  and  equivalent  plastic  strain  contours  for  this  case  are  very 
similar  to  Figs.  5  and  6  and  are  not  reproduced  here.  Plastic  strain  rate  con¬ 
tours  are  plotted  in  Fig.  8  at  a  load  level  of  J\  —  181  kJ/m2.  Plastic  strain 
rates  along  a  line  at  6  —  45°  plotted  in  Fig.  9  show  the  good  agreement  be¬ 
tween  the  finite  element  results  and  the  /-field  solution  of  Eq  10.  Some  scat- 
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DISTANCE  FROM  CRACK  TIP,  R 

FIG.  9 — Plastic  strain  rate  versus  distance  from  crack  tip  at  0  =  45°;  Ji  =  181  kJ/m \ 
Ji  =  1.63  X  10!  kJ/m2s;  rate-insensitive  dynamic  analysis  S16. 

analysis  did  not  account  for  the  possibility,  and  effects,  of  adiabatic  heating 
in  the  very  rapidly  strained  region  at  the  crack  tip.) 

Analysis  S24  is  similar  to  S16,  except  that  a  strain-rate-dependent  flow 
stress  was  used.  The  character  of  the  equivalent  stress,  equivalent  plastic 
strain,  and  plastic  strain  rate  contours  (Figs.  10  to  12)  are  similar  in  form  to 
the  static  and  the  dynamic  (rate-insensitive)  case,  but  the  magnitudes  of  these 
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FIG.  11 —Contours  of  plastic  strain:  Ji  =  188  kJ/m \  Ji  =  1.77  X  10 '  kJ/m:s;  rate-sensitive 
dynamic  analysis  S24. 

quantities  are  quite  different.  It  is  interesting  that  the  magnitudes  of  all  these 
quantities  are  higher  than  for  the  static  case  at  a  comparable  /i-level.  This  is 
largely  due  to  the  change  in  the  “effective”  stress-strain  curve  resulting  from 
the  strain-rate  effects  (Fig.  2).  Figure  13  compares  the  equivalent  stresses  and 
equivalent  plastic  strains  from  the  dynamic  rate-sensitive  analysis  (S24)  with 
the  dynamic  rate-insensitive  analysis  (S16)  along  a  line  at  8  —  45°.  The  plotted 
points  are  from  the  finite  element  results,  and  the  curves  are  from  the  /-field 
solution  (the  rate-sensitive  curve  was  developed  by  using  the  modified  a'  and 
rf  parameters).  These  results  all  show  that  the  effect  of  strain-rate-sensitive 


FIG.  12 — Contours  of  plastic  strain  rate ;  J\  =  188  kJ/m:,  Ji  =  1.77  X  10'  kJ/m:s;  rate-sensi¬ 
tive  dynamic  analysis  S24. 
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FIG.  13 — Comparison  of  equivalent  stress  and  plastic  strain  from  the  fast  dynamic  rate-sensitive 
analysis,  S24,  with  the  fast  dynamic  rate-insensitive  analysis,  St 6. 


flow  stress  is  to  significantly  alter  the  magnitudes  of  the  stresses  and  strains 
at  the  tip  of  a  crack. 

The  constant-rate  stress-strain  curve,  plotted  as  a  dashed  line  in  Fig.  2, 
does  not  map  the  Finite  element  data  in  Fig.  13  when  used  in  the  ./-field  equa¬ 
tions.  The  equivalent  stress  is  underpredicted,  particularly  where  the  plastic 
strains  are  large. 

Analyses  SI 2  and  S20  are  much  more  slowly  loaded  dynamic,  inertialess 
computations.  The  loading  rate  is  K\  =  1.0  MPa\/m/s,  which  is  within  the 
range  of  0.55  <  K\  <  2.75  MPa\/m/s  used  in  ASTM  Test  for  Plane-Strain 
Fracture  Toughness  of  Metallic  Materials  (E  399).  Even  at  this  slow  loading 
rate,  there  is  a  difference  of  roughly  10%  in  equivalent  stress  and  2%  in  equiv¬ 
alent  plastic  strain,  as  shown  in  Fig.  14  along  a  line  at  6  =  45°. 

The  crack  tip  opening  displacement  has  been  used  as  a  fracture  initiation 
parameter.  Figure  15  plots  the  crack  tip  opening  displacement  6  as  a  func- 
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DISTANCE  FROM  CRACK  TIP,  R 

FIG.  14 — Comparison  of  equivalent  stress  and  plastic  strain  from  the  slow  dynamic  rate-sensi¬ 
tive  analysis.  S20.  with  the  slow  dynamic  rate-insensitive  analysis,  SI 2. 

tion  of  J\.  It  is  clear  that  the  strain-rate-dependent  flow  stress  has  a  signifi¬ 
cant  effect  on  6.  Consider  the  rate-insensitive  analysis,  SI 6.  At  the  point  in¬ 
tersected  by  the  dashed  line,  J\  —  182  kJ/mJ  and  6  =  237  pm.  The  average 
plastic  strain  rate  is  1.06/s.  This  strain  rate  would  raise  the  flow  stress  by  a 
factor  of  1.22  according  to  the  Malvern  equation  (Eq  14),  with  D  =  62920/s 
and  p  —  7.245,  and  if  the  material  were  rate-sensitive.  This  increase  in  flow 
stress  would  correspond  to  a  22%  larger  J\  value  according  to  Eq  13,  if  6 
were  constant.  Figure  15  shows  that  the/i  value  for  the  dynamic  rate-sensitive 
analysis,  S24,  is  in  fact  25%  larger  than  the  rate-insensitive  analysis.  A  less 
pronounced  difference  is  found  between  the  slowly  loaded  analyses,  S12  and 
S20,  where  a  4%  change  in  o0  and  J\  is  predicted  by  Eq  13  but  is  not  verified 
by  finite  elements,  since  this  is  within  the  noise  level  of  the  computation.  This 
suggests  that  a  rapidly  loaded  rate-sensitive  ductile  material  could  fail  at  a 
higher  J\c  than  a  slowly  loaded  one,  if  a  critical  crack  tip  opening  displace- 
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J-INTEGRAL  [kJ/m*] 

FIG.  15 — Crack  tip  opening  displacement  versus  3-integral  curves  for  the  fast  dynamic  rate  sen¬ 
sitive  analysis.  S24.  and  fast  dynamic  rate-insensitive  analysis.  S16. 


ment  is  the  failure  criterion.  The  work  of  Little  et  al  [5]  showed  the  same  re¬ 
sult,  albeit  that  computations  were  performed  in  the  10~4  <Ki<  102 
MPa \J m/s  range. 

The  /'-integral,  discussed  earlier  in  this  paper,  was  also  calculated  for  each 
dynamic  analysis.  The  kinetic  energy  density  term  is  very  small  along  the  in¬ 
tegration  contours  close  to  the  crack  tip,  and  the  inertia  term  is  also  very 
small  due  to  the  small  accelerations  produced  by  the  inertialess  loading  tech¬ 
niques.  As  a  result,  the  differences  between  Rice’s  /-integral  and  the  /’-integral 
are  less  than  0.09%  for  the  analyses  considered  in  this  study.  Consequently, 
the  treatment  of  these  dynamic  problems  by  using  the  (static)  /-field  solu¬ 
tions  of  Eqs  5  and  6  is  justified  in  the  case  of  inertialess  loading. 

The  authors  recognize  that  actual  dynamically  loaded  structures  and  lab¬ 
oratory  specimens  do  not  exhibit  the  analytical  tractability  of  inertialess  nu¬ 
merical  models.  Determination  of  the  stress  and  strain  fields  is  complicated 
by  vibrations  that  create  localized  unloading  and  even  reverse  plasticity.  The 
/-field  solutions  are  not  applicable  in  these  cases,  and  Rice’s  /-integral  is  no 
longer  path-independent.  As  a  first  step  in  the  dynamic  analysis  of  cracked 
bodies  where  inertial  effects  are  significant,  the  authors  have  performed  pre¬ 
liminary  elastic-plastic  finite  element  analyses.  The  /'-integral  is  path- 
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independent  in  these  cases  and  shows  promise  as  an  elastoplastic  stress  field 
characterizing  parameter,  since  it  has  already  been  used  to  describe  the  stress 
state  in  elastodynamically  propagating  cracks  [29].  The  details  are  still  unre¬ 
solved  of  how  the  stress  state  can  be  described  using  f  in  the  elastoplastic 
dynamic  case,  providing  a  challenging  topic  for  ongoing  research. 

Conclusions 

A  technique  has  been  described  in  which  finite  element  analyses  were  per¬ 
formed  on  dynamically  loaded  models  without  any  inertial  effects.  This  “in¬ 
ertialess  loading  technique”  permitted  quantitative  determination  of  the  ef¬ 
fects  of  strain-rate-dependent  flow  stress. 

Another  technique  was  developed  that  enabled  the  stresses,  strains,  and 
crack  opening  displacement  to  be  approximated  by  a  simple  modification  of 
the  a  and  n  parameters  in  the  power  stress-strain  law.  Using  the  modified 
parameters,  a'  and  n' ,  in  the  /-field  solution,  stress  and  strains  came  to  within 
10%  of  a  corresponding  finite  element  computation,  in  regions  where 
the  /-field  was  applicable.  This  technique  was  only  usable  in  conjunction 
with  the  inertialess  loading  technique. 

Finite  element  computations  revealed  that  stresses  and  strains  in  the  crack 
tip  plastic  zone  were  higher  in  a  rate-sensitive  material  than  in  a  rate- 
insensitive  material  at  corresponding  Ji-values.  Crack  tip  opening  displace¬ 
ments,  however,  were  lower  in  the  rate-sensitive  material.  The  same  behavior 
was  predicted  using  the  a,  ^'-modified  /-field  solution. 

Finite  element  computations  also  revealed  that  if  a  critical  crack  tip  open¬ 
ing  displacement  were  used  as  a  failure  criterion,  then  the  strain-rate- 
dependent  material  could  withstand  a  higher  /-value  before  this  5C  was 
reached.  Equation  13  showed  that  this  increase  in  /ic  could  be  predicted 
based  on  how  the  (area)  average  plastic  strain  rate  would  alter  the  flow 
stress. 

The  newly  developed  /'-integral  was  determined  from  the  finite  element 
results.  In  the  case  of  a  stationary  crack  under  inertialess  dynamic  loading  it 
is  virtually  identical  to  Rice’s  (static)  /-integral.  The  /’-integral  does  appear 
promising  in  general  dynamic  loading  situations  with  propagating  cracks, 
since  it  remains  path-independent,  whereas  Rice’s  /-integral  does  not. 
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Abstract — Rapid  fracture  in  ductile  viscoplastic  materials  is  discussed  in  this  paper.  Following  a 
review  of  previous  studies  of  rapid  fracture  in  ductile  materials,  a  finite  element  model  is  developed 
for  use  with  a  crack  extension  algorithm  discussed  in  a  companion  paper.  Despite  the  lack  of  an 
established  rapid  ductile  fracture  criterion,  computations  which  assumed  that  the  6-\a  history 
was  the  same  for  rapid  fracture  as  it  was  for  stable  crack  growth  agree  qualitatively  with 
experimental  data  It  is  shown  that  viscoplastic  materials  absorb  considerably  more  energy  during 
rapid  fracture  than  rate  insensitive  materials  and  are  therefore  “tougher”.  A  detailed  examination 
of  the  stresses  in  the  immediate  vicinity  of  the  crack  tip  reveals  a  constantly  changing  stress  field 
that  does  not  reach  a  "steady-state"  as  the  crack  tip  translates. 


NOMENCLATURE 

Ap  crack  tip  plastic  zone  area 

a  crack  length 

a  crack  velocity 

Aa  crack  growth 

C,  shear  wave  speed  of  the  material;  C,*=(p  p)'"1 

D  constant  in  viscoplastic  constitutive  relation  (3) 

£  Young's  modulus 

Ev  undetermined  function  of  0  and  sometimes  v,  m  or  a’C, 

Jt  Rice's  path  independent  contour  integral;  energy  release  rate  under  quasistatic  loading 

Jlc  critical  value  of  ./-integral:  when  crack  propagation  initiates 

J' i  Atluri's  path  independent  coniour  integral;  energy  release  rate  under  dynamic  loading 

K  amplitude  of  stress  or  strain  singularity 

A,  mode  I  stress  intensity  factor 

m  ratio  of  tangent  modulus  to  elastic  modulus  for  a  linear  hardening  material 

n  strain  hardening  exponent 

n,  direction  cosines 

p  exponent  in  viscoplastic  constitutive  relation  (3) 

r,  6  polar  coordinates 

s  arc  length 

s  strength  (exponent!  of  the  stress  and  strain  singularity  for  rapidly  propagating  crack;  given  inequation  (1,2) 

T  kinetic  energy  density;  T=u,u,Z 

t  traction  vector;  l,  =  c,lnJ 

u  dimensionless  function  of  6  and  a 

u  displacement  vector 

u  velocity  vector 

ii  acceleration  vector 

I’  volume 

vL  load  point  displacement 

vL  load  point  velocity 

w  :• 

strain  energy  density;  H  =|  ervdcv 
.0 

cr  crack  tip  opening  angle 

o,  critical  crack  tip  opening  angle 

T  integration  contour 

S  crack  tip  opening  displacement 

c,;  strain  tensor 

t,.r  dimensionless  functions  of  8  and  a 
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4,  yield  strain:  c^=  o„  E 

t,  equivalent  plastic  strain  rate 

if...  average  plastic  strain  rate 

H  shear  modulus 

Y  undetermined  function  of  6  and  sometimes  v.  m  or  a  C, 

Oy  Cauchv  stress  tensor 

Oy  dimensionless  functions  of  6  and  a 

om  mean  or  hydrostatic  stress;  +  a,  +  c.)  '3 

c0  yield  stress  under  static  tensile  loading 

a o,,  a.  normal  stresses  in  the  x.  y  and  r-directions 


INTRODUCTION 

Rapidly  propagating  cracks  in  rate  sensitive  (viscoplastic)  ductile  materials  are  examined  in  this 
paper.  Despite  the  practical  importance  of  this  class  of  fracture  problem,  particularly  in  relation 
to  crack  arrest  phenomena,  relatively  little  analytical  or  numerical  work  has  been  done  in  this 
field.  After  a  brief  review  of  previous  work,  the  fracture  criteria  and  finite  element  model  will  be 
described.  The  results  of  rapid  crack  growth  simulation  with  crack  velocities  up  to  100  m  s  reveal 
that  viscoplastic  materials  are  significantly  more  resistant  to  crack  propagation  than  rate  insensitive 
materials. 

The  analysis  of  dynamic  fracture  in  rate  sensitive  ductile  materials  is  complicated  by  at  least 
three  nonlinearities:  the  material  is  ductile  so  fracture  is  accompanied  by  large  amounts  of  plastic 
flow,  load  rates  and  crack  speeds  are  sufficiently  large  to  necessitate  the  inclusion  of  inertia  effects 
in  the  analysis,  and.  many  ductile  materials  of  industrial  interest  have  strain  rate  dependent  flow 
properties.  Other  nonlinearities  such  as  finite  deformation  effects  and  adiabatic  heating  at  the  tip 
of  a  rapidly  propagating  crack  could  be  investigated,  but  these  will  not  be  included  in  this  study. 

The  stress  and  strain  fields  at  the  tip  of  a  rapidly  propagating  crack  in  a  ductile  rate  insensitive 
material  were  characterized  by  Achenbach  and  Kanninen  [1]  Iot  the  Mode  III  case,  and  later  by 
Achenbach  er  al.  [2]  for  the  Mode  I  case.  For  an  elastic-linear  strain  hardening  material: 

o,i  =  A'  ct0  r’  £„  (ft  v,  m.  a! Q  (I) 

e,j  =  K  a0  r'  E0  ( ft  v.  m.  a  Cs)  .  (2) 

The  order  of  stress  and  strain  singularity  is  in  the  range  —  0.5  <  5  <  0.  The  asymptotic  approach 
used  to  determine  eqs  (1 )  and  (2)  is  unable  to  give  values  for  the  magnitude  factor.  K.  which  must 
be  determined  by  other  means. 

Finite  element  modelling  of  rapid  ductile  fracture  was  performed  by  Ahmad  er  al.  [3]  and 
later  by  Barnes  et  al.  [4]  at  Battelle  Columbus  Laboratories.  Motivated  by  previous  discoveries 
of  Kanninen  er  al.  [5]  that  quasistatic  analyses  or  elastodynamic  analyses  grossly  underestimate 
the  resistance  of  AISI  4340  steel  to  crack  grow  th  under  impact  loading,  the  Battelle  group  performed 
dynamic  finite  element  simulations  of  3-point  bend  specimens  using  experimentally  determined 
elastic-strain-hardening-plastic  materia!  properties.  A  critical  crack  tip  opening  angle  was  used 
as  the  fracture  criterion,  and  crack  extension  was  modelled  using  the  "node  release"  technique 
described  by  various  researchers  [6-10].  Their  numerical  results  for  macroscopic  parameters  such 
as  crack  position  as  a  function  of  time,  and  dynamic  stress  intensity  factor  agreed  very  well  with 
the  corresponding  experiments.  However,  due  to  the  coarseness  of  the  computational  grid,  no 
results  for  near  tip  stresses  or  strain  rates  were  obtained. 

Lo  [1 1]  recently  presented  solutions  to  the  problem  of  rapid  crack  extension  in  ductile  materials 
with  rate  sensitivity.  A  semi-infinite  crack  propagating  in  the  v-direction  under  small  scale  yielding 
conditions  was  modelled.  An  asymptotic  analysis,  in  the  sense  of  Amazigo  and  Hutchinson  [12] 
and  Achenbach  et  al.  [2],  was  performed  on  the  basis  of  crack  velocities  less  than  the  shear  wave 
velocity  of  the  material  The  analysis  uses  one  of  the  elastic-viscoplastic  constitulive  relations 
described  by  Perzyna  [13]  which  does  not  include  strain  hardening  but  models  associated 
viscoplastic  flow.  The  Perzyna  relations  are  stated  in  a  general  three-dimensional  form,  of  which 
the  one-dimensional  case  specializes  to 
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The  parameters  <j0.  D  and  p  must  be  determined  experimentally  [14]. 

The  near  tip  stresses,  strains  and  displacements  were  determined  based  on  the  assumption 
that  they  would  be  singular  in  r: 

o,j  =  r'r-'  du(6,  a)  (4) 

e%  =  r1 e p  (6.  a)  (5) 

1  h  1 r-i 

—  rr  ir  i  u ,(Q.  a) ,  16) 

£  D£ 

where  a  is  the  crack  velocity.  B  =  Dja{,.  D  and  p  are  defined  by  eq.  (3),  £  is  Young’s  modulus,  r 
and  0  are  coordinates  centered  on  the  moving  crack  tip,  and  <ry.  and  c  pv  and  tJ,  are  dimensionless 
functions  of  6  and  a. 

The  possibility  of  using  the  7-integral,  or  any  dynamic  variation  of  the  7-integral  idea,  as  a 
near  tip  characterizing  parameter  is  completely  ruled  out  by  eqs  (4-6).  The  7-integral  is  path 
independent  for  stationary  cracks  since  it  is  comprised  of  energy  density  terms,  order  r~ ',  integrated 
around  a  closed  path,  order  r1.  which  result  in  the  order  r°  integral  7.  (Energy  density  terms  are 
of  the  form  ct„  c„and  are  gi\en  by  the  HRR  singularity  solution  as  being  singular  in  r  to  the  order 
-  1)  In  the  rate  sensitive  rapid  ductile  crack  extension  case,  eq  (4-6)  would  result  in  the  energy 
density,  order  r~2r~l.  being  integrated  around  a  path,  order  r.  giving  a  quantity  of  order  rr~'r~\ 
Therefore  as  r  —  0.  (for  p  >  3)  this  contour  integral  also  tends  to  zero,  not  a  desirable 
property  for  a  path  independent  integral.  Some  form  of  7  may  still  be  used  as  a  (relatively)  far 
field  characterizing  parameter,  providing  it  is  calculated  along  a  contour  where  no  unloading  takes 
place.  That  is  not  unlike  the  justification  for  using  7,  to  predict  tearing  instability  [15]  or  using 
7,  in  the  presence  of  crack  tip  blunting  [16,  17], 

The  Lo  results  also  reveal  that  the  hydrostatic  or  mean  norma]  stresses. 
cr„  =  (cr,,  +  cr,,  +  a..)  3.  ahead  of  the  propagating  crack  tip.  decreases  with  increasing  crack 
velocity.  As  a  C,  increases  from  0  to  0.61.  decreases  in  om  of  40%  and  28%  are  reported  for  p  -  4 
and  p  -  6.  respectively.  It  was  shown  by  McClintock  [18]  and  Rice  and  Tracey  [19]  that  void 
and  hole  growth  is  enhanced  by  high  hydrostatic  tensile  stresses  (neglecting  the  issue  of  whether 
sufficient  time  is  available  for  this  to  occur  in  dynamic  situations).  Since  ductile  fracture  is  the 
end  result  of  a  hole  growth  process,  any  phenomenon  which  affects  hole  growth  will  affect  ductile 
fracture.  However,  it  is  difficult  to  predict  what  affects  the  interrelation  between  hydrostatic  stresses, 
critical  strain  to  fracture,  and  crack  velocity,  will  have  on  a  fracture  parameter  such  as  the  critical 
crack  tip  opening  angle. 

The  only  reported  finite  element  study  of  rapid  crack  propagation  in  rate  sensitive  ductile 
materials  was  performed  by  Brickstad  [20],  The  Brickstad  analysis  makes  use  of  the  Perzyna  [13] 
elastic- viscoplastic  constitutive  model  with  linear  strain  hardening,  for  a  high  strength  carbon 
steel  (ct0  =  1485  MPa).  The  "energy  flow  into  the  crack  tip  region"  was  used  as  the  fracture 
criterion,  and  crack  extension  was  modelled  using  the  "node  release"  technique  [6-10].  The  energy 
flow  into  the  crack  tip  region  is  a  mesh  size  dependent  parameter,  but  previous  finite  element 
calculations  of  this  parameter  by  Rydholm  el  al.  [21]  had  agreed  well  with  theoretical  predictions 
in  elastodynamic  problems. 

Brickstad  [20]  used  a  rather  coarse  mesh,  and  could  not  comment  on  the  near  tip  stress  and 
strain  distributions.  However  he  noted  strain  rates  of  the  order  10’/s  at  the  Gauss  points  nearest 
the  crack  tip  propagating  at  250-750  m  s. 

The  analytical  work  of  Lo  [1 1]  and  the  finite  element  analyses  of  Brickstad  [20]  suggest  that 
the  fracture  resistance  of  viscoplastic  materials  is  enhanced  at  high  crack  velocities.  The  finite 
element  analyses  undertaken  in  this  study  quantitatively  assess  the  increased  fracture  resistance 
of  a  viscoplastic  material:  A533B  steel  at  93  C.  and  elucidate  details  of  the  crack  tip  stress  field 
during  rapid  crack  propagation. 
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FINITE  ELEMENT  MODELLING 

The  finite  element  model  used  in  the  analysis  of  rapid  crack  growth  is  depicted  in  Fig.  1.  The 
material  is  modelled  as  linear  elastic  up  to  the  yield  point,  and  then  with  power  law  hardening 
beyond  first  yield: 
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Bo 
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— ,  a  <  ob 
a0 
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Bo 


<T  >  CT0 


(8) 


with  ct0  =  382.866  MPa,  £  =  197.62  GPa,  %  -  a0/E,  and  n  =  10.  The  rate  dependent  analyses 
employed  the  rate  sensitivity  model  stated  by  Perzyna  in  ref.  [13]. 


£P  —  D 


(9) 


with  D  =  62919.55  and  p  =  7.24545.  Close  correspondence  of  the  analytical  model  with 
measurements  of  the  flow  properties  of  A533B  steel  at  93CC  is  demonstrated  in  a  separate  paper 
[14].  While  it  is  still  debatable  whether  eq.  (9)  is  the  “correct”  representation  of  the  constitutive 
behaviour  of  A533B  steel  at  very  high  strain  rates,  it  does  give  good  agreement  with  experimental 
results  in  the  range  <  104/s. 

Any  crack  growth  simulation  must  incorporate  a  method  for  continuously  changing  the 
boundary  conditions.  A  companion  paper  [22]  presents  a  detailed  discussion  about  how  spring 
and  gap  elements  may  be  used,  in  the  general  purpose  program  ABAQUS,  to  simulate  crack 
extension. 

An  initial  velocity  field  was  imposed  on  the  specimen  at  the  beginning  of  the  analysis.  This 
was  to  prevent  the  occurrance  of  large  load  impulses  at  the  beginning  of  the  analysis.  Impulsive 
loads  induce  large  amplitude  vibrations  in  the  model  which  overwhelm  the  inertial  and  viscoplastic 
effects  at  the  crack  tip.  The  same  technique  was  used  successfully  in  the  analysis  of  rapidly  loaded 
stationary  cracks  [23]. 
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Fig  1  Finite  element  model  of  compact  tension  specimen  used  in  the  investigation  of  propagating  cracks 
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Four  analyses  were  performed:  a  rate  insensitive  one  and  a  viscoplastic  one  at  a  moderate 
loading  rate  of  v*  =  1  m/s,  as  well  as  a  rate  insensitive  analysis  and  a  viscoplastic  analysis  at  a 
higher  loading  rate  of  =  100  m/s.  These  analyses  are  summarized  in  Table  1.  The  computation 
time  required  for  these  simulations  was  typically  12-16  CPU  hr.  Although  these  computations 
were  performed  on  a  Decsystem-10  computer  with  a  Tops- 10  operating  system,  it  was  later  shown 
that  very  similar  times  could  be  expected  on  a  VAX  11/750  running  the  VMS  operating  system. 


Table  1.  Summary  of  finite  element  analyses  for  propagating  cracks 


Analysis 

identification 


Analysis 

type 

Static 

Dynamic 

Dynamic 

Dynamic 

Dvnamic 


Load  line 
velocity 
(m  si 
0 

1.0 

1.0 

100.0 

100.0 


k, 

at  start 
(MPa^ms) 
0 

1.1  x  10s 
1.1  x  10’ 
1.1  x  10’ 
1.1  x  101 


Rate 

dependent? 

No 

No 

Yes 

No 

Yes 


IMPLEMENTATION  OF  FRACTURE  CRITERIA 

The  necessity  of  using  a  particular  fracture  criterion  in  a  finite  element  crack  growth  simulation 
forces  one  to  make  some  a  priori  assumptions  about  the  characterization  of  the  ductile  fracture 
process  under  rapid  loading.  Whereas  stable  crack  growth  can  be  characterized  by  an  experimentally 
obtained  Jrt±a  relationship  just  after  initiation  [7,  24],  very  few  JrAa  results  are  available  for 
rapid  crack  growth.  In  the  absence  of  such  experimental  data,  some  other  fracture  criterion  must 
be  used  at  initiation  and  immediately  following  initiation. 

The  crack  tip  opening  displacement,  S ,  shows  promise  as  a  near  tip  characterizing  parameter 
since  it  can  be  interpreted  as  the  integral  of  strains  around  the  crack  tip  [25].  The  S-&a  history, 
used  as  the  fracture  criterion  for  the  first  4  mm  of  growth  in  the  dynamic  analysis,  was  obtained 
from  the  results  of  the  stable  crack  growth  analysis  [22].  Figure  2  illustrates  that 

S=  -  33.58  (Ao) :  +  0.47082  Ao  +  2.5332  x  I0'4  (10) 


0  J  a  %  8  10  >2  I«  16 

CRACK  GROWTH,  Oo  [mm] 

Fig  2  Crack  tip  opening  displacement  vs  crack  growth 
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(where  S  and  Aa  are  in  meters),  for  Aa  <  4  mm.  A  critical  crack  tip  opening  angle  was  used  as 
the  fracture  criterion  for  Aa  >  4  mm.  The  stable  crack  growth  value  of  ac  —  0.24  radians  was 
used  as  the  fracture  criterion.  The  use  of  rate  invariant  values  of  S  (Aa)  and  ac  is  equivalent  to 
assuming  that  the  local  ductile  fracture  mechanism  is  rate  insensitive.  This  assumption  cannot  be 
defended  at  this  time.  It  is  made  so  that  its  consequences  can  be  tested  against  experiments. 

The  computational  algorithms  which  performed  the  rapidly  growing  crack  simulations  were 
almost  identical  to  those  used  for  stable  crack  growth.  However,  the  first  4  mm  of  crack  growth 
were  ^-controlled  rather  than  ./-controlled,  as  discussed  in  the  preceding  paragraphs. 

RESULTS  AND  DISCUSSION 

The  results  of  the  rapid  crack  growth  simulations  are  presented  in  this  section.  Experimental 
results  [24]  and  numerical  results  [22],  pertaining  to  slow  stable  crack  growth,  are  also  presented 
where  possible  and  appropriate. 

Figure  3  is  a  plot  of  as  a  function  of  crack  growth.  It  is  known  [24]  that  y,c  for  the  slow- 
stable  crack  growth  analysis  is  200  kJ  m  2.  Yet,  fracture  initiated  in  the  rate  insensitive  dynamic 
analysis  at  (J{  =  )  205  kJ  m  :  and  223  kJ/m  2  for  vt  =  1.0  m/s  and  \L  =  100.0  m/s,  respectively; 
and  fracture  initiated  in  the  viscoplastic  analyses  at  257kJ/m-  and  330kJ/m2  for  vL  =  1.0  m/s 
and  vL  =  100.0  m/s,  respectively.  Therefore,  the  effect  of  the  local  material  inertia  is  not  very 
significant  at  initiation,  whereas  the  viscoplastic  effects  are  quite  large. 

The  rate  insensitive  dynamic  results  for  »£  =  1.0  m/s  are  virtually  identical  to  the  slow  stable 
crack  growth  results.  However,  values  of  J{  from  the  viscoplastic  analysis  are  always  about  17% 
larger  than  the  stable  crack  grow  th  values.  The  same  trend  is  evident  in  the  iL  =  100.0  m/s  results: 
the  rate  insensitive  results  average  out  to  about  the  same  as  the  stable  crack  growth  values,  yet 
the  viscoplastic  results  are  about  30%  larger. 

The  observation  that  rapid  crack  propagation  in  ductile  materials  must  overcome  more 
resistance  than  slow  crack  propagation  was  also  observed  by  Kanninen  ei  al.  [5],  Although  the 
Kanninen  et  al.  study  is  inconclusive  as  to  the  source  of  the  difference  in  toughness  between  slow 
and  rapid  loading,  the  finite  element  results  of  this  paper  do  suggest  that  viscoplasticity  could 
account  for  the  difference. 

The  first  2.5  mm  of  crack  growth  are  replotted  in  Fig.  4  and  compared  with  experimentally 
determined  ./-resistance  curves  obtained  by  Joyce  [26].  Joyce  loaded  1  j2T  compact  tension 
specimens  (W  =  25.4  mm  I  of  A533B  steel  at  a  rate  of  600  in  's  whereby  k,  at  the  start  of  the  test 
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Fig  3  J  -integral  vs  crack  growth 
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Fig  4.  Static  and  dynamic  7-resistance  curves 
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is  about  22%  higher  than  for  the  \L  =  lOOm/s  cases  studied  here.  Joyce’s  experimental  results 
show  more  of  a  rate  sensitive  effect  and  this  may  be  due  to  (i)  Joyce  tested  at  150  C,  whereas  the 
finite  element  solutions  assumed  a  temperature  of  93CC,  (ii)  Joyce  tested  1/2"  thick  specimens, 
whereas  the  finite  element  solution  is  based  on  a  plane  strain  model,  (in)  crack  tip  blunting  is 
counted  as  crack  growth  in  the  slow  stable  crack  growth  results  of  Joyce,  and,  (iv)  the  inertia  of 
the  specimen  appears  as  increased  resistance  to  crack  extension  in  the  experimental  results. 

The  iL  =  100.0  m/s  results  have  a  high  frequency  oscillation  superimposed  on  an  otherwise 
monotonically  increasing  response.  This  is  likely  caused  by  the  crack  extension  procedure  which, 
by  incrementally  relaxing  the  boundary  conditions  at  the  crack  tip  nodes,  puts  small  impulses  into 
the  crack  tip  region  at  each  time  increment.  Evidence  that  the  crack  extension  procedure  is  the 
source  of  these  oscillations  is,  that  during  the  first  millimeter  of  crack  growth,  the  oscillations  are 
of  a  higher  “frequency”  than  later  on.  Secondly,  the  “wavelength”  of  oscillations  is  about  1  mm 
(for  Aa  >  1  mm)  where  the  element  size  is  1  mm.  This  crack  extension  technique  may,  therefore, 
be  limited  to  load  point  velocities  below ,  say,  500  m/s.  This  comment  also  applies  to  the  node 
release  technique  mentioned  in  refs  [6-8,  10].  However,  the  moving  singularity  element  of  Atluri 
and  Nishioka  [27-29]  does  show  promise  in  this  area,  if  a  correct  viscoplastic  formulation  can 
be  found. 

If  the  value  of  J,  as  stated  by  Rice. 


J\  =  j.  (H>i,  -  ds. 


(ID 


decreases  as  a  function  of  time  for  the  case  of  a  stationary  crack,  then  the  specimen  must  be 
unloading  and  the  use  of  J,  is  invalid  In  contrast,  the  value  of  J,  as  stated  by  Atluri  [30], 

J'\  =  I  [  (W  +  T)n]  -  r,i/,,]  dy  +  I  -  u,u  ,]dT.  (12) 

r  -  I 


may  decrease  during  the  load  history  without  the  specimen  necessarily  unloading.  Ai  high  loading 
rates,  the  inertia  term  ( f ,  p[ii,.  .  .  .]  d I ' )  of  eq.  (12)  is  comparable  in  magnitude  to  the  strain  energy 
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and  traction  terms.  Small  changes  in  the  material  point  velocities  will  make  the  inertia  term 
negative,  and  reduce  J{,  even  though  the  bulk  of  the  specimen  is  still  loading.  This  is  the 
situation  observed  in  the  =  lOO.Om/s  analyses,  and  is  the  source  of  the  oscillations  in  these 
results.  The  J  -integral  values  quoted  in  this  paper  are  calculated  along  the  contour  shown  by 
the  dashed  line  in  Fig.  1  Values  of  J‘  calculated  at  contours  interior  to  this  one  are  path  independent 
(to  within  +  10%)  until  the  points  at  which  the  integration  is  performed  experience  unloading 
due  to  the  effects  of  the  growing  crack. 

Figure  5  is  a  plot  of  the  crack  tip  opening  angle  as  a  function  of  crack  growth.  The  crack 
tip  opening  angle  is  poorly  defined  during  the  first  millimeter  of  crack  growth,  as  in  the  case  of 
stable  crack  growth.  Between  1  and  4  mm  of  crack  growth,  the  crack  tip  opening  angle  is  very 
similar  for  the  quasistatic  and  dynamic  computations.  This  is  likely  due  to  the  fact  that  the  S-\a 
relationship  from  the  stable  crack  growth  analysis  was  used  as  the  fracture  criterion  for  the  first 
4  mm  of  crack  growth  in  the  dynamic  analyses.  After  4  mm  of  crack  growth,  all  of  the  analyses 
use  the  same  value  of  u  =  0.24  radians  as  the  fracture  criterion.  The  viscoplastic  vL  =  1.0  m/s 
analysis  has  two  anomalous  points  at  A  a  =  0.125  and  2.0  mm.  After  verifying  that  these  data 
points  are  indeed  correct,  we  can  still  offer  no  explanation  for  them. 

Figures  6  and  7  are  plots  of  load  point  displacement  and  load  as  a  function  of  crack  growth, 
respectively.  The  rate  insensitive  iL  —  1.0  m  s  results  are  virtually  the  same  as  the  slow  stable 
crack  growth  results.  Since  the  load  point  displacement  and  load  are  higher  in  the  viscoplastic 
l"i  =  1-0  m/s  analysis  (2%  and  20%.  respectively),  an  analysis  which  does  not  address  viscoplastic 
effects  for  this  case  would  seriously  underestimate  the  energy  dissipated  during  rapid  crack  growth. 
The  iL  =  100.0  m  s  results  have  a  different  character  than  the  slower  loading  rate  results.  The 
general  form  of  the  load  vs  Aa  curves  in  Fig  7  is  that  of  a  sine  wave  with  a  period  of  190  ps  (i.e. 
twice  the  time  for  about  6  mm  of  crack  growth).  This  implies  that  the  natural  frequency  of  the 
overall  specimen  is  about  5.3  kHz,  which  agrees  very  well  with  a  handbook  value  for  an  cantilever 
beam  with  the  same  length  as  the  arms  of  the  compact  tension  specimen.  This  illustrates,  not 
surprisingly,  that  the  bulk  motion  of  the  specimen  is  greatly  affected  by  the  dynamic  response  of 
the  specimen.  The  viscoplastic  \L  —  100.0  m  s  analysis  shows  again  that  considerably  more  energy 
is  dissipated  due  to  viscoplasticity  than  could  be  accounted  for  in  a  rate  insensitive  elastoplastic 
dynamic  analysis 

The  normalized  crack  velocity,  a  \L,  is  plotted  as  a  function  of  crack  growth  in  Fig.  8.  For 
the  four  cases  here  it  appears  that 

a  ~  (o  +  fi  Ao)  \  L  (13j 

=  (0.4  +  1.40  Ao)  ii 
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Fig  5  Crack  up  opening  angle  \s  crack  growth 
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Fig.  6.  Load  point  displacement  vs  crack  growth. 


( Aa  is  in  meters).  The  constants  a  and  /?  are  geometry  dependent:  indeed,  a  and  /?  reflect  the  fact 
that  a  compact  tension  specimen  is  kinematically  equivalent  to  two  cantilever  beams.  It  appears 
that  a  and  /3  are  insensitive  to  inertia,  load  rate  and  strain  rate. 

Stresses  ahead  of  the  propagating  crack  tip  are  shown  in  Fig.  9—1 1 .  Unless  otherwise  stated 
in  this  discussion  the  word  “stress”  shall  refer  to  equivalent  stress,  oe;  “strain”  shall  refer  to 
equivalent  plastic  strain.  eP:  and  “strain  rate”  shall  refer  to  equivalent  plastic  strain  rate,  eP.  At 
initiation,  the  stresses  ahead  of  the  crack  tip  are  not  a  function  of  the  material  inertia  so  the 
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Fig  8  Normalized  crack  velocity  vs  crack  growth. 


vL  =  1.0  m/s  results  are  the  same  as  the  t£  =  100.0  m/s  results  in  Fig.  9.  This  idea  was  proved  by 
Achenbach  ei  al.  [31]  as  r  -*  0,  and  introduced  into  the  finite  element  computation  using  the 
inertialess  loading  technique  discussed  in  ref.  [23].  Figure  9  shows  that  the  stress  is  elevated  by 
about  21%  for  the  vL  =  100.0  m/s  analysis  due  to  viscoplastic  effects. 

Figure  10  is  a  plot  of  the  stresses  for  the  four  dynamic  analyses  after  4  mm  of  crack  growth. 
The  results  are  not  nearly  as  smooth  as  the  initiation  results.  The  oscillations  in  stress  could  be 
due  to  noise  in  the  computation,  but  are  more  likely  due  to  the  dynamic  response  of  the  material 
near  the  crack  tip.  The  stress  oscillations  are  very  similar  for  the  two  \L  =  100.0  m/s  analyses, 
and  for  the  two  =  1.0  m/s  analyses,  and  are  therefore  not  due  to  random  noise.  The  viscoplastic 
effects  again  raise  the  magnitude  of  the  stresses  appreciably.  The  notion  of  Lo  [11]  that 


Fig  9  Equivalent  stress  vs  distance  from  crack  tip  at  fracture  initiation 
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Fig.  10.  Equivalent  stress  vs  distance  from  crack  tip  after  4  mm  of  crack  growth. 

viscoplasticity  makes  the  near  tip  stress  fields  more  singular  cannot  be  addressed  on  the  bases  of 
Fig.  9-11  since  this  would  require  a  crack  tip  mesh  at  least  an  order  of  magnitude  finer  than  the 
one  used  here. 

Figures  11  is  a  plot  of  the  stresses  after  8  mm  of  crack  growth.  Results  were  only  available 
for  the  iL  =  100.0  m/s  case.  Again  significant  elevation  in  stress  is  apparent  due  to  viscoplastic 
effects. 

Figures  12-14  are  plots  of  the  strain  ahead  of  the  propagating  crack  tip  and  correspond  to 
the  stress  plots  in  Figs.  9-11.  The  inertia  of  the  materia)  near  the  crack  tip  does  not  appear  to 
influence  the  strain.  However,  the  rate  sensitivity  increases  the  strain,  a  small  amount  at  initiation, 
a  little  more  after  4  mm  of  crack  growth,  and  more  still  after  8  mm  of  crack  growth,  as  shown  in 
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Fig  11  Equivalent  stress  vs  distance  from  crack  Up  after  8  mm  of  crack  growth 
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Fig.  ]2.  Plastic  strain  vs  distance  from  crack  tip  at  fracture  initiation. 


Figs  12-14.  This  is  contrary  to  the  results  of  creep  studies  by  Little  et  al  [32]  which  indicated 
smaller  plastic  strains  due  to  increasing  strain  rate. 

The  results  presented  in  Figs  9-14  do  not  represent  the  crack  tip  stress  or  strain  fields  on  a 
sufficiently  fine  size  scale  to  permit  conclusions  to  be  drawn  about  what  effects  inertia  and  rale 
sensitivity  have  on  the  critical  strain  to  fracture  in  the  fracture  process  zone.  The  results  do  give 
a  good  representation  of  stress  and  strain  in  the  region  around  the  crack  tip  where  the  7-field 
solutions  would  be  applicable  if  the  cracks  were  stationary. 

The  plastic  strain  rates  are  plotted  as  a  function  of  distance  from  the  crack  tip  in  Figs  15- 
17.  In  general,  the  iL  =  100.0  m  s  strain  rates  are  2  orders  of  magnitude  higher  than  the  \L  =  1.0  ms 
strain  rates.  The  strain  rates  at  initiation  are  well-behaved,  whereas,  the  strain  rates  after  4  or 


Fig  13  Plastic  strain  \s  distance  from  crack  tip  after  4  mm  of  crack  growth 
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Fig.  14.  Plastic  strain  vs  distance  from  crack  tip  after  8  mm  of  crack  growth. 


8  mm  of  crack  growth  vary  dramatically.  Small  errors  in  strain  tend  to  be  magnified  when 
calculating  strain  rates;  this  is  partly  the  cause  of  the  large  strain  rate  fluctuation.  It  is  therefore 
difficult  to  generalize  these  results.  (Note  that  data  points  which  would  fall  below  the  axes  shown 
in  these  figures  have  been  omitted.  Their  omission  is  indicated  by  joining  the  adjacent  data  points 
with  dotted  lines.) 

Figure  18  is  a  plot  of  the  average  plastic  strain  rate,  defined  by 
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Fig  15  Plastic  strain  rate  vs  distance  from  crack  tip  at  fracture  initiation 
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Fig  16.  Plastic  strain  rale  vs  distance  from  crack  tip  after  4  mm  of  crack  growth 


as  a  function  of  crack  growth.  The  average  plastic  strain  rates  are  about  1/s  and  1 00  s  for  the 
vL  =  1.0  ms  analysis  and  the  \L  =  100.0  m/s  analysis  respectively.  Substituting  these  values  into 
the  Perzyna  equation  (3)  implies  that  the  stresses  increase  by  22%  and  41%,  respectively,  due  to 
viscoplastic  effects.  This  is  substantiated  by  Fig.  9  at  initiation,  and  to  a  lesser  degree  by  Figs  10 
and  1 1  after  4  or  8  mm  of  crack  growth.  This  supports  the  observation  made  previously  [23]  that 
rate  effects  in  the  bulk  of  the  specimen  can  be  estimated  based  on  the  knowledge  of  an  average 
strain  rate. 

The  deformed  crack  profiles  are  shown  in  Fig.  19  for  the  \L  =  100.0  m  s  viscoplastic  analysis. 
A  startling  feature  is  the  large  peak  located  at  .v  =  0,  which  develops  after  crack  initiation.  This 
must  be  due  to  viscoplastic  effects,  since  it  was  not  observed  in  the  rate  insensitive  analyses  or 


Fig  1 1  Pljstic  strain  rate  vs  distance  from  crack  up  after  8  mm  of  crack  growth 
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Fig.  18.  Average  plastic  strain  rate  vs  crack  growth. 


the  stable  crack  growth  analysis  [22],  Due  to  the  very  large  strain  rates  at  initiation  (ep  as  12,000/s), 
the  stress  is  elevated  about  80%  over  the  rate  insensitive  value,  and  therefore  considerably  more 
unloading  takes  place  at  the  initial  crack  tip  than  at  the  surrounding  material. 


CONCLUSIONS 


The  technique  for  modelling  crack  propagation  using  spring  and  gap  elements,  described  in 
a  companion  paper  [22],  was  successfully  applied  to  rapid  crack  propagation  in  viscoplastic 
materials. 


Fig  19  Deformed  crack  profiles  for  100  m  s  load  point  velocity,  viscoplastic  analys 


aaafflsaiaiaa&^^ 


A 


r 


mimi  languor  mi  mi  eji  rs*.  r^K  rv*  rs#.  n^jvuvi  ~^rv^-ry  jv^jwmrv.^rv 


m  rwwvwvi 


| 


■ 

i 

V 


460 


R.  HOFF  ei  al 


Due  to  the  lack  of  an  established  fracture  criterion  for  rapid  ductile  fracture,  the  3-Aa  history 
from  a  stable  crack  growth  and  analysis  [22]  was  employed  for  the  first  4  mm  of  crack  growth, 
whereafter  the  crack  tip  opening  angle  was  used  as  the  fracture  criterion.  This  assumption  implies 
that  the  initial  strain  to  fracture,  as  represented  in  an  integral  sense  by  5,  is  insensitive  to  strain 
rate  and  inertia.  The  effect  of  loading  rate  as  the  calculated  resistance  curves  agree  qualitatively 
with  experimental  results  by  Joyce  [26]  (who  tested  at  a  different  temperature).  However,  the  fact 
that  the  measurements  seem  to  show  a  greater  rate  sensitivity  than  the  calculations  may  assign 
that  S-Aa  and  ac  display  a  positive  rate  sensitivity. 

Computations  for  crack  velocities  around  1.0  m/s  and  lOOm/s  reveal  that  the  viscoplastic 
materials  absorb  considerably  more  energy  during  fracture  than  rate  insensitive  materials.  This 
difference  increases  at  higheT  crack  velocities.  This  is  evidenced  by  the  higher  J-Ao  curve  at  higher 
loading  rates. 

For  the  four  cases  studies,  the  crack  velocity  normalized  with  respect  to  the  load  point  velocity, 
a ;yl,  is  a  linear  function  of  crack  growth  and  insensitive  to  inertia  or  strain  rate.  However,  the 
normalized  crack  velocity  is  likely  geometry  and  load  history  dependent. 

A  detailed  examination  of  the  stresses  in  the  5  mm  region  ahead  of  the  crack  tip  revealed 
that  stresses  were  increased  due  to  strain  rate  effects.  Contrary  to  slow  stable  crack  growth  case, 
where  a  fixed  stress  distribution  appears  to  translate  with  the  crack  tip  after  an  initial  transient, 
it  was  not  observed  that  the  stresses  reached  a  “steady  state"  distribution  after  4  or  8  mm  or  crack 
growth.  The  dynamic  response  of  the  material  near  the  crack  tip  seems  to  preclude  the  determination 
of  such  a  solution. 
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Abstract — This  paper  presents  a  new  technique  for  simulating  crack  extension  in  conjunction 
with  the  finite-element  method.  The  technique  uses  spring  and  gap  elements  to  control  the 
motion  of  nodes  on  the  crack  plane.  These  elements  are  available  in  many  proprietary  finite- 
element  codes,  thereby  obviating  the  need  for  a  user-written  finite-element  code.  Numerical 
results  for  stable  crack  growth  are  in  excellent  agreement  with  corresponding  experiments.  The 
technique  is  also  applied  to  rapid  fracture  in  ductile  materials,  as  discussed  in  a  companion 
paper. 


'  NOMENCLATURE 

a  crack  length 
a  crack  velocity 

>  A  a  crack  growth 

E  Young’s  modulus 
h  element  size 
J\  Rice's  contour  integral 
n  strain-hardening  exponent 

u.  v,  w  displacements  in  the  X,  Y  and  Z  directions,  respectively 
x,  y,  z  coordinate  directions 
a  crack-tip  opening  angle 
8  crack-tip  opening  displacement 
«  strain 
to  yield  strain 
a  stress 

t  (To  yield  stress 

INTRODUCTION 

i 

The  analysis  of  stable  crack  growth  in  ductile  materials  using  finite  elements  is  a  necessary 
precursor  to  the  finite-element  analysis  of  rapid  crack  growth  in  these  same  materials.  Not  only 
is  stable  crack  growth  a  special  case  of  rapid  crack  propagation  (i.e.  as  a  —*  0),  but  many  of 
the  techniques  developed  for  modelling  stable  crack  growth  can  be  extended  to  the  rapid- 
propagation  case.  After  a  review  of  previous  work  in  this  field,  a  new  technique  for  modelling 
stable  crack  growth  will  be  described  in  this  paper.  Finally,  stable-crack-growth  results  will  be 
presented  based  on  an  analysis  using  this  new  technique.  A  companion  paper[l]  describes  the 
application  of  this  technique  to  rapid  viscoplastic  fracture. 

Other  investigators  have  taken  basically  two  approaches  to  modelling  stable  crack  growth 
with  finite  elements;  we  shall  call  them  the  microscopic  approach  and  the  macroscopic  ap¬ 
proach.  The  microscopic  approach  models  a  region  in  the  immediate  vicinity  of  the  crack  tip. 
The  elements  at  the  crack  tip  are  much  smaller  than  the  crack-tip  opening  displacement  at 
initiation  and  the  entire  mesh  is  usually  much  smaller  than  the  crack  length.  With  this  fine- 
mesh  gradation,  it  is  possible  to  model  the  large  strains  in  the  fracture  process  zone  (especially 
if  a  large-strain  finite-element  formulation  is  used).  It  is  also  feasible  to  use  microstructurally 
significant  fracture  criteria  such  as  the  “critical  plastic  strain  at  a  critical  distance”  criterion 
of  McClintock  and  Irwin[2]  (also  discussed  in  Rice[3]).  The  microscopic  approach  can  give 
accurate  representations  of  the  near-tip  stress  and  strain  singularities  for  limited  amounts  of 
crack  growth.  Examples  of  research  involving  the  microscopic-crack-growth  approach  are  An- 
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derson[4],  Rice  et  al.[ 5].  Rice  and  Sorensen[6],  and  Sorensen[7,  8].  Unfortunately,  the  near¬ 
tip  stress  fields  elucidated  using  this  approach  are  often  difficult  to  relate  to  experimentally 
measureable  quantities  such  as  load,  load-point  displacement,  or  even  the  J  integral  (where 
applicable). 

The  macroscopic  approach  to  modelling  stable  crack  growth  provides  data  which  are  easily 
related  to  experimentally  measured  parameters.  The  finite-element  mesh  typically  models  an 
entire  laboratory-sized  specimen  or  a  relatively  large  portion  of  a  structural  component.  The 
size  of  the  elements  near  the  crack  tip  is  of  the  order  of  the  crack-tip  opening  displacement  at 
initiation,  or  larger.  Since  all  dimensions  of  the  model  are  very  large  compared  to  microstruc- 
turally  significant  distances,  the  fracture  criterion  must  also  be  meaningful  and  calculable  on 
that  size  scale.  Currently  the  most  frequently  used  are  the  J  integral,  the  crack-tip  opening 
angle  a,  and  a  combination  of  the  two.  Since  the  critical  values  of  J  and  a  are  initially  determined 
experimentally,  it  is  therefore  not  surprising  that  crack-growth  simulations  using  the  macro¬ 
scopic  approach  and  these  fracture  criteria  correspond  well  with  experimental  measurements. 

The  macroscopic-crack-growth  approach  has  received  far  more  attention  than  the  microscopic  Jyfli 
approach  and  representative  studies  include  work  by  Andersson[9],  Du  and  Lee[  10,1 1  ] ,  Hellan 
and  Lotsberg[12],  Kanninen  et  a/.[13],  and  Shih  et  al. [14]. 

A  finite-element  crack-growth  simulation  must  employ  some  mechanism  for  altering  the 
boundary  conditions  (once  the  fracture  criteron  is  satisfied)  to  simulate  crack  extension.  By 
far  the  most  commonly  used  method  is  the  “node-release”  technique.  Conceptually  speaking, 
when  the  fracture  criterion  is  satisfied  the  displacement  constraint  at  the  crack-tip  node  is 
replaced  by  an  equivalent  force.  This  nodal  force  is  then  gradually  reduced  to  zero  in  subsequent 
increments  of  the  analysis,  allowing  the  node  in  question  to  displace,  thereby  making  the  next 
node  the  new  crack  tip.  This  procedure  is  repeated  for  as  many  subsequent  nodes  as  required 
for  the  desired  amount  of  crack  growth.  Andersson[4.  9,  15],  Du  and  Lee[10,  11],  Hellan  and 
Lotsberg[12],  Kanninen  et  al. [13],  Light  et  al. [16,  17],  Rice  et  al.[ 5],  Rice  and  Sorensen[6], 
and  Sorensen[7,  8]  have  all  used  variations  of  this  technique. 

An  alternative  technique,  used  successfully  by  Shih  et  al. [14,  18],  is  to  shift  the  crack-tip 
node  in  the  direction  of  crack  growth,  thereby  modelling  the  lengthening  crack.  However,  this 
technique  is  still  unproven  in  a  dynamic  context  so  it  is  not  used  in  this  study.  Atluri,  Nishioka 
and  coworkers[  19-21]  have  developed  special  crack-tip  elements  which  contain  embedded 
singularities,  and  translate  with  the  moving  crack  tip.  Although  this  technique  proved  very 
successful  in  elastodynamic  problems[20,  22-28],  it  has  some  theoretical  shortcomings  in  stable- 
crack-growth  problems,  since  the  order  of  the  crack-tip  strain  singularity  must  be  known  a 
priori  in  order  to  formulate  the  element.  Since  the  order  of  the  strain  singularity  is  velocity 
dependent  in  rapid  crack  propagation  (see  Achenbach  et  al. [29]).  the  application  of  the  moving- 
singularity  element  to  elastoplastic  dynamic  problems  is  even  less  accurate. 

The  finite-element  model  and  methodology  discussed  in  this  paper  were  developed  in  view 
of  two  constraints:  (i)  the  methodology  developed  for  modelling  stable  crack  growth  must  be 
useable  for  elastoplastic  dynamic  crack  growth,  and  (ii)  these  computations  must  be  performed 
using  a  commercially  available  finite-element  program.  Stable  crack  growth  is  the  zero-velocity 
special  case  of  rapid  crack  propagation  in  ductile  materials.  Since  the  computations  are  more 
easily  performed  in  the  absence  of  inertia  effects,  stable  crack  growth  was  chosen  as  the  starting 
point  for  the  finite-element  analysis  of  rapid  crack  growth  in  ductile  materials.  Most  of  the 
finite-element  work  performed  in  a  production  environment  is  done  with  commercially  available 
finite-element  codes.  However,  no  results  of  crack-growth  simulations  have  been  reported 
which  were  obtained  using  any  of  these  commercially  available  programs  (Shih  et  al. [14]  used 
ADINA[30];  however,  the  program  was  specially  modified  for  this  purpose).  It  was  felt  that 
the  development  of  crack-growth  methodology,  which  could  be  used  with  a  proprietary  finite- 
element  program,  would  be  of  more  utility  than  a  methodology  which  also  required  a  user- 
developed  finite-element  code.  The  proprietary  code  ABAQUS[31]  was  chosen  for  this  work 
because  of  its  proven  capabilities  in  the  solution  of  nonlinear  problems.  The  installation  at 
Vanderbilt  University  featured  ABAQUS  version  3-12-78,  running  on  a  DEC-1099  computer, 
with  the  TOPS-10  operating  system.  The  technique  described  in  this  paper  should  also  work 
equally  well  with  other  proprietary  finite-element  programs. 
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A  finite-element  model  was  developed  to  permit  the  investigation  of  stable  and  rapid  crack 
growth  in  ductile  steels.  Figure  I  is  a  plot  of  the  undeformed  mesh.  The  cutout  in  the  top  figure 
is  shown  enlarged  below.  In  both  cases  the  initial  crack  tip  is  indicated  by  an  arrow.  The 
dimensions  of  the  model  are  identical  to  specimen  T52,  from  a  recent  study  by  Shih  ei  al. [14], 
at  General  Electric.  The  dashed  lines  show  the  contours  where  the  J  integral  will  later  be 
calculated. 

The  material  is  modelled  as  linear  elastic  up  to  the  yield  point,  with  power-law  hardening 
beyond  first  yield: 
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Values  for  the  material  constants  were  chosen  appropriate  for  A533B  steel  at  93°C:  a0  =  382.866 
MPa,  E  =  197.62  GPa.  e0  =  cro/E.  and  n  =  10.  Since  ABAQUS  requires  a  multilinear  stress- 
strain  curve,  the  power-law  portion  was  discretized.  The  stress-strain  curve  was  discretized 
in  an  optimal  sense  so  that  deviations  were  held  to  within  r  0.5%  (see  Hoff  et  a/. [32]). 

A  series  of  38  springs  is  defined  from  points  in  space  to  ground.  The  purpose  of  these 
springs  is  to  associate  a  stiffness  with  degree-of-freedom  number  2  (i.e.  y  displacement)  for 
nodes  60005  to  60042  inclusive  (see  Fig.  2).  In  fact,  the  v  coordinate  of  nodes  60005  to  60042 
is  initially  zero,  so  that  those  nodes  begin  on  the  crack  plane.  Nodes  60005  to  60042  are  shown 
elevated  from  the  crack  plane  in  Fig.  2  for  clarity.  Gap  elements  are  defined  between  nodes 
60005  and  50005.  60006  and  50006.  etc.  up  to  60042  and  50042.  The  gap  elements  are  formulated 
so  that  they  can  transmit  only  compressive  loads,  i.e.  no  tensile  or  shear  loads.  Since  the 
tendency  of  nodes  50005-50042  is  to  move  upwards,  due  to  tensile  loading  of  the  specimen, 
all  of  the  gaps  will  close.  A  displacement  boundary  condition  is  applied  to  nodes  60005-60042. 
constraining  these  nodes  from  moving  in  the  v  direction.  With  all  of  the  gaps  closed,  nodes 


Fig.  1.  Finite-element  model  of  compact  tension  specimen  used  in  the  investigation  of  propa¬ 
gating  cracks. 


50005-50042  are  also  constrained  to  the  crack  plane,  thereby  locating  the  crack  tip  at  node 
50042.  It  would  not  be  possible  to  grow  the  crack  with  displacement  boundary  conditions  applied 
to  nodes  50005-50042.  since  the  motion  of  these  nodes  after  crack  growth  cannot  be  determined 
a  priori. 

A  static  elastic  analysis  was  performed  to  verify  the  adequacy  of  the  spring-and-gap  method 
of  specifying  the  boundary  conditions.  The  y  displacements  of  the  nodes  constrained  by  the 
gaps  were  all  of  order  10  ~20.  Since  the  displacements  of  the  remainder  of  the  nodes  in  the 
mesh  were  of  order  I0~".  one  can  conclude  that  the  gaps  only  allow  negligible  displacements. 
By  opening  and  closing  the  appropriate  gaps,  the  crack  length  was  varied  from  1 17.2  to  129.2 
mm.  For  each  different  crack  length,  values  of  the  stress-intensity  factor  A'i  were  determined 
from  the  finite-element  results  using  two  different  methods  and  compared  with  the  polynomial 
equation  of  Srawley[33]: 

(2  +  fl/u)[0.866  +  4.64«/ir  -  I3.32(«At)-  +  14.72(«/W')-  -  5.6(fl/W)4] 

(1  -  a!W)y-  '  <3) 

Equation  (3)  is  accurate  to  +0.5%  in  the  range  0.2  <  a/W  <  1.  The  values  of  Aj  determined 
using  the  stiffness-derivative  method  of  Parks[34]  agree  with  eqn  (3)  to  within  0.2%.  This  result 
verifies  that  the  spring-and-gap  method  is  capable  of  modelling  displacement  constraints  ac¬ 
curately.  Unfortunately  the  stiffness-derivative  method  has  not  been  generalized  to  elastoplastic 
dynamic  problems  and  will  not  be  used  henceforth.  The  J  integral  of  Rice[35]  is.  of  course, 
suitable  to  elastoplastic  problems,  and  variations  such  as  yj  [36.  37]  are  appropriate  for  dynamic 
problems.  Values  of  Aj  were  also  determined  by  performing  the  contour  integration  for  Jil38]. 
but  agree  less  well  with  the  polynomial  equation  ( +  5.2%).  This  is,  however,  sufficient  accuracy 
for  the  crack-growth  simulations  to  be  performed  in  this  study. 

Values  of  J i  or  K\  could  have  been  determined  more  accurately  using  contour  integration 
if  singular  elements  were  used  at  the  crack  tip.  However,  the  presence  of  singular  elements 
caused  convergence  difficulties  with  the  gaps  when  this  was  attempted.  Moreover,  the  standard 
singular  element,  created  by  collapsing  one  side  of  an  eight-noded  isoparametric  finite  element 
and  shifting  the  adjacent  nodes  to  the  quarter  points{39,  40],  did  not  offer  any  advantage  in 
crack-growth  studies. 

Numerical  experimentation  revealed  that  the  choice  of  spring  constant  affected  the  con¬ 
ditioning  of  the  stiffness  matrix[4j.  It  was  determined  that  if  the  spring  constant  was  numerically 
equal  to  Young’s  modulus  then  the  best  performance  was  obtained.  For  instance,  if  E  =  200 
x  109  Pa.  then  the  spring  constant  should  be  set  to  200  x  I09  N/m.  Very  stiff  springs  ti  e. 
12-18  orders  of  magnitude  stiffer  than  the  recommended  value)  did  not  cause  any  detectable 
degradation  in  solution  accuracy,  whereas  very  soft  springs  (i.e.  the  stiffness  is  12  or  more 
orders  of  magnitude  less  than  the  recommended  value)  significantly  degraded  the  results.  It  is 
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difficult  to  determine  exact  limits  on  acceptable  spring  stiffnesses  since  the  conditioning  of  a 
stiffness  matrix  depends  on  the  number  of  degrees  of  freedom  in  the  analysis  and  the  number 
of  springs.  For  the  problems  worked  in  this  paper,  spring  stiffnesses  within  6  orders  of  magnitude 
of  the  recommended  value  were  employed  and  produced  no  deleterious  effect  on  the  results. 

IMPLEMENTATION  OF  THE  FRACTURE  CRITERION 

Several  studies!  13,  14]  have  focused  attention  on  which  fracture  parameters  are  suitable 
for  describing  stable  crack  growth.  Candidate  parameters  such  as  .A,  a,  8,  cL//da,  average  crack¬ 
opening  angle,  generalized  energy  release  rate,  etc.  were  assessed  on  the  basis  of  a  large  number 
of  requirements:  not  only  must  the  ideal  parameter  be  physically  significant,  but  it  must  also 
be  readily  determinable  using  finite  elements.  It  was  found[13.  14]  that,  for  a  given  material 
and  small  amounts  of  crack  growth,  the  J  integral  has  a  unique  relationship  to  the  amount  of 
crack  growth,  A  a.  This  relationship  is  only  mildly  specimen  dependent.  The  crack-tip  opening 
angle  a  initially  has  a  very  high  value,  but  after  a  small  amount  of  crack  growth  reduces  to  a 
constant  value.  In  this  study  the  7|-Aa  relationship  determined  from  experiment!  14]  was  used 
for  the  first  4  mm  of  crack  growth,  whereafter  a  constant  crack-tip  opening  angle  was  used  for 
the  remaining  crack  growth.  This  two  parameter  fracture  criterion  is  depicted  in  Fig.  3.  where 
the  first  regime  is  termed  ./-controlled  and  the  second  regime  is  termed  a-controlled. 

The  major  difficulty  encountered  in  modelling  ./-controlled  crack  growth  centers  on  this 
paradox:  the  boundary  conditions,  which  are  a  function  of  A  a  which  in  turn  is  only  a  function 
of  J i,  must  be  specified  at  the  beginning  of  an  increment  in  the  solution  but  Jt  can  only  be 
determined  at  the  end  of  that  particular  increment.  This  would  not  constitute  a  major  problem 
if  boundary  conditions  could  be  handled  implicit y  at  each  increment  in  a  finite-element  analysis 
(i.e.  the  boundary  conditions  are  determined  as  part  of  the  solution  process).  However,  bound¬ 
ary  conditions  must  always  be  stated  explicitly  in  a  finite-element  analysis.  This  necessitates 
the  use  of  a  rather  sophisticated  extrapolation  scheme  to  predict  what  the  value  of  J\  will  be 
at  the  end  of  a  particular  increment.  The  following  paragraphs  will  describe  the  evolution  of 
the  Ji  extrapolation  technique. 
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Fig.  3.  Comparison  of  regimes  of  crack  growth  where  J \  and  a  are  appropriate  characterizing 

parameters. 
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The  specimen  was  loaded  so  that  the  load-point  displacement  varied  linearly  with  a  non- 
dimensional  “time."  Since  Kj  varies  linearly  with  load-point  displacement  in  the  elastic  case, 
and  J\  varies  quadratically  with  Kx,  therefore  J\  should  vary  quadratically  with  time.  Referring 
to  Fig.  4(a),  if  the  solution  increment  just  completed  is  at  then  a  parabola  was  fit  through 
the  values  of  Jx  at  times  /,_2.  /,-i,  and  The  upcoming  increment  would  end  at  t,- ,.  and  J, 
could  be  predicted  by  extending  the  parabola  to  /,■♦  i .  This  technique  proved  to  be  very  unstable. 
Generally  the  predicted  value  of  J\  tended  to  be  much  too  great  (Fig.  4(a))  or  much  too  small 
(Fig.  4(b)).  The  latter  case  was  more  serious  since  the  desired  crack  length  is  determined  from 
the  predicted  Jx  value,  and  if  the  predicted  Jx  value  was  less  than  the  previous  J\  value  then 
the  crack  would  become  shorter  instead  of  longer. 

It  turns  out  that  in  an  elastic-plastic  analysis  J i  varies  almost  linearly  with  load-point 
displacement  (with  some  perturbation  due  to  the  changing  stiffness  of  the  model  as  the  crack 
grows),  so  the  assumed  quadratic  behaviour  was  not  a  good  model.  A  linear-extrapolation 
scheme  was  therefore  attempted,  whereby  J\  at  /,■*  i  was  determined  based  on  data  taken  at 
ti- i  and  As  shown  in  Fig.  4(c),  linear  extrapolation  was  more  stable  than  quadratic  extrap¬ 
olation  but  minor  perturbations  in  Jx  were  still  amplified  in  time. 

The  technique  finally  implemented  was  a  modification  of  the  linear-extrapolation  method. 
Instead  of  using  data  at  i  and  the  predicted  value  of  Jx  was  determined  by  extrapolating 
data  at  2  and  In  addition,  the  slope  of  the  J\  vs  t  curve  was  not  allowed  to  change  by 
more  than  20 9c  between  successive  increments.  This  produced  the  most  stable  results  of  all 
and  is  depicted  in  Fig.  4(d).  The  Jx-t  history  is  shown  in  Fig.  5  for  ./-controlled  crack  grow  th, 
and  it  is  evident  that  this  is  a  very  stable  technique. 

Once  the  value  at  J\  was  predicted,  the  value  of  A  a  was  determined  from  a  parabolic  fit 
to  experimental  data  from  a  General  Electric  study[14).  A  a  was  the  smaller  root  of 

J\  =  -2.562  x  10lo(Aa)2  +  3.985  x  108  Aa  +  2.000  x  105,  (4) 


•  CALCULATED  VALUE 
+  PREDICTED  VALUE 


(b) 


Fig.  4.  Comparison  of  J\  prediction  schemes:  (a)  Parabolic  extrapolation  overestimating  J,  (. 
(b)  parabolic  extrapolation  underestimating  J,  *  i ,  (c)  linear  extrapolation  from  points  at  t,  ,  and 
t,.  and  (d)  linear  extrapolation  from  points  at  and 
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Fig.  5.  J i  history  for  the  ./-controlled  portion  of  stable  crack  growth. 


where  A  a  is  in  units  of  meters  and  J\  is  in  units  of  J/m2.  The  value  of  A  a.  determined  from 
eqn  (4).  was  never  equal  to  an  integer  multiple  of  the  element  length. 

This  raises  the  question  of  how  to  model  these  crack  lengths  which  do  not  coincide  with 
element  boundaries.  Suppose  that  the  “effective"  crack  length  A  a  would  put  the  crack  tip  at 
a  position  as  shown  in  Fig.  6.  Clearly  nodes  n,  n  -  1,  n  -  2.  etc.  must  be  constrained  in  the 
y  direction.  This  is  done  by  leaving  nodes  m.  m  -  1.  m  -  2,  etc.  on  the  crack  plane  as  they 
have  been  since  the  beginning  of  the  analysis.  The  gaps  between  n  and  m,  n  -  1  and  m  -  1. 
etc.  will  remain  closed,  thereby  constraining  nodes  n,  n  -  1,  n  -  2.  etc.  to  remain  on  the 
crack  plane.  Nodes  n  +  3 .  n  +  4.  n  +  5.  etc.,  which  previously  were  used  to  constrain  n  + 
3,  n  +  4.  etc.,  are  now  displaced  upwards  so  that  the  gaps  are  open.  Spring  node  m  +  2  is 
displaced  upwards  in  proportion  to  how  far  the  effective  crack  tip  has  advanced  through  element 
/:  e.g.  if  the  crack  tip  has  advanced  25%  of  the  way  between  nodes  n  +  2  and  n,  then  node  m 
+  2  is  raised  25%  of  its  maximum  value.  The  maximum  displacement  of  m  +  2  is  almost  twice 
what  the  Final  displacement  of  node  n  +  2  will  be.  Therefore,  about  halfway  through  the 
incremental  gap-opening  procedure,  the  node  n  +  2  will  separate  from  the  gap  (i.e.  the  gap 
opens)  and  node  n  +  2  will  Find  its  own  equilibrium  position.  In  general,  the  gaps  were  closed 


Fig.  6.  Details  of  crack-extension  procedure. 
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for  three  increments  and  then  open  for  another  two  increments  during  the  incremental  opening 
procedure. 

When  the  crack  tip  moved  past  node  n,  then  nodes  m  -  1  and  m  were  moved  upwards 
incrementally.  This  procedure  repeated  itself  until  4  mm  of  crack  growth  had  occurred.  At  that 
time  the  critical  value  of  the  crack-tip  opening  angle.  0.24  rad.  was  used  as  the  fracture  criterion 
for  further  crack  growth.  Opening  of  the  crack-tip  element  was  handled  using  springs  and  gaps, 
as  just  described,  except  that  the  determination  of  the  effective  crack  length  was  slightly  dif¬ 
ferent.  The  motion  of  the  spring  nodes.  .  .  .  m  -  1 ,  m.  m  +  1  .  .  .  .  are  specified  in  a  user- 
developed  subroutine  which  is  attached  to  the  ABAQUS  program  at  link  time.  A  listing  of  the 
subroutines  used  in  this  study  appears  in  Ref.  [38]. 

RESULTS  AND  DISCUSSION 

The  results  of  a  simulation  of  stable  crack  growth  are  presented  in  this  section.  An  ex¬ 
haustive  presentation  is  not  made;  however,  sufficient  evidence  is  given  to  demonstrate  that 
the  spring-and-gap  technique  is  appropriate  for  modelling  crack  extension.  In  addition,  the 
stresses  and  strains  ahead  of  the  stably  growing  crack  are  contrasted  with  those  ahead  of  a 
stationary  crack. 

The  variation  of  J\  with  crack  growth  is  plotted  in  Fig.  7.  The  values  of7i  are  determined 
by  integration  along  the  contour  shown  in  Fig.  1.  Contours  much  closer  to  the  crack  tip  were 
generally  unacceptable  for  J\  evaluation  since  material  unloading  took  place  behind  the  crack 
tip  as  crack  extension  progressed.  Unloading  violates  one  of  the  conditions  of  using  J2  flow 
theory  in  the  finite-element  analysis  to  evaluate  the  J integral,  which  is  based  on  the  deformation 
theory  of  plasticity.  Besides,  as  pointed  out  by  McMeeking[42]  on  the  basis  of  finite-strain 
studies,  the  J  integral  is  not  path  independent  when  calculated  along  contours  closer  than  about 
58  from  the  crack  tip. 

The  variation  of  J\  with  A  a  corresponds  very  well  with  the  experimental  results  of  Shih 
et  <z/.[14],  Since  the  crack  growth  is  determined  during  the  analysis  from  eqn  (4).  it  is  not 
surprising  that  there  is  almost  an  exact  correspondence  to  experiment  during  the  first  4  mm  of 
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Fig.  7.  J  integral  i  j  crack  growth. 
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considerable  noise  in  a  at  the  start,  however,  and  this  may  be  due  to  a  number  of  factors.  First, 
it  can  be  shown  that  at  initiation  the  7-field  solution  predicts  -it  rad  as  the  tip  opening  angle. 
Since  a  is  determined  according  to 


a  =  tan 


it  is  not  uniquely  defined  at  initiation.  The  element  size  h  and  the  v  displacement  of  node  n  ~ 
2,  vn  —  2 ,  are  illustrated  in  Fig.  6.  Second,  due  to  the  mesh  refinement  near  the  initiation  site, 
the  element  size  h  keeps  changing  for  the  first  5  elements  of  crack  growth. 

The  steady-state  value  of  the  crack-tip  opening  angle.  0.24  rad,  was  obtained  graphically 
by  determining  the  asymptote  of  the  a  vs  t  variation  during  the  first  4  mm  of  crack  growth 
(although  this  plot  is  not  shown).  This  value  is  slightly  higher  than  the  experimental  steady- 
state  value  of  0.21  rad,  which  is  not  unexpected  since  the  numerical  values  of  a  are  almost 
all  higher  than  the  experimental  values  at  less  than  4  mm  of  crack  growth. 

The  variation  of  8  with  An  is  plotted  in  Fig.  9.  In  this  instance,  8  is  twice  the  v  displacement 
of  the  node  1  mm  behind  the  original  crack  tip.  In  the  7-controlled  region,  the  calculated  values 
deviate  by  less  than  2.2%  from  a  curve  defined  by 

8  =  - 33.58(An)2  +  0.47082An  +  2.5332  x  10~4.  (6) 


where  8  and  A  a  are  in  units  of  meters.  In  the  a-controlled  region,  the  calculated  values  deviate 
by  less  than  1.3%  from  a  line  defined  by 

8  =  0.2D673Afl  +  7.6159  x  I0'4.  (7) 


Although  the  Jrla  relationship  in  Fig.  7  is  similar  to  this,  7|  vs  8  does  not  plot  as  a  straight 
line  over  the  entire  range  of  crack  growth. 
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Fig.  II.  Load  vs  crack  growth. 


The  plot  of  load-point  displacement  vL  vs  da  shown  in  Fig.  10,  exhibits  a  similar  bimoda! 
(parabolic/linear)  characteristic.  The  load-point  displacement  is  defined  as  twice  the  y  displace¬ 
ment  of  the  node  at  the  top  of  the  load  hole.  Since  the  finite-element  model  is  only  of  the  top 
half  of  the  specimen,  this  factor  of  2  is  necessary  for  correspondence  with  experiment  where 
vL  is  defined  as  the  relative  displacement  of  the  top  loading  pin  with  respect  to  the  bottom  one. 
The  slope  of  the  linear  portion  in  Fig.  10  is  significant  because  its  value,  0.834,  can  be  used  to 
estimate  the  crack  velocity  in  a-controlled  rapid  crack  extension.  That  is,  if  d(A a)ldvL  =  1/ 
0.834  for  stable  growth,  then  the  crack  velocity  a  can  be  estimated  from  the  load-point  velocity 
i>L  according  to 


t'L 

0.834 


1.20  i'L. 


(8) 


Departures  from  the  behaviour  predicted  by  eqn  (8)  would  be  due  to  inertia  and  strain-rate 
effects. 

The  variation  of  load  with  da  is  plotted  in  Fig.  1 1.  The  finite-element  results  were  deter¬ 
mined  by  assuming  a  1-m-thick  plane-strain  specimen  with  25%  side  grooves.  The  experimental 
results  were  therefore  multiplied  by  a  suitable  constant  to  allow  comparison  with  the  numerical 
results.  The  two  results  unquestionably  show  the  same  trends.  Any  quantitative  differences 
are  due  to  the  computational  model  being  (intentionally)  slightly  different  from  the  experimental 
specimen  in  terms  of  material  properties,  fracture  criteria,  etc. 

The  deformed-crack  profiles  are  shown  in  Fig.  12  for  the  first  4  mm  of  crack  growth.  Note 
that  the  scales  are  different  on  the  abscissa  and  ordinate  so  that  the  crack  opening  is  exaggerated 
in  the  figure.  Several  features  of  this  plot  are  noteworthy.  The  character  of  deformation  changes 
abruptly  since  the  strain  fields  are  quite  different  at  initiation  from  what  they  are  after  finite 
amounts  of  growth.  Sorensen(7]  and  Miller  and  Kfouri[43]  noted  similar  behaviour  in  the  de¬ 
formed-crack  profiles.  A  more  puzzling  phenomenon  is  the  jagged  nature  of  the  deformed 
profiles.  This  is  perhaps  due  to  the  nature  of  the  gap-release  mechanism. 
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The  equivalent  stresses  and  plastic  strains  at  45°  to  the  tip  of  the  growing  crack  are  plotted 
in  Figs.  13  and  14.  The  stresses  and  strains  are  virtually  the  same  after  4  mm  of  crack  growth 
as  they  are  after  10  mm  of  crack  growth,  indicating  that  the  4  mm  point  gives  indicative  results 
for  stable  crack  growth.  The  stresses  and  strains  are  generally  higher  after  stable  growth  than 


Fig.  13.  Equivalent  stress  along  a  line  at  45°  to  the  propagating  crack  tip. 
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Fig.  14.  Plastic  strain  along  a  line  at  45°  to  the  propagation  crack  tip. 


they  are  at  initiation,  but  are  much  less  singular.  These  smaller  stress  and  strain  singularities 
for  stable  growth  were  anticipated  on  the  basis  work  by  Amazigo  and  Hutchinson[44],  Although 
the  data  in  Figs.  13  and  14  do  not  prove  the  existence  of  a  weaker  crack-tip  singularity,  they 
do  indicate  the  anticipated  trend. 


CONCLUSIONS 

A  new  technique  was  described,  employing  spring  and  gap  finite  elements,  which  can  be 
used  to  simulate  crack  growth.  An  advantage  of  this  technique  is  that  it  does  not  require  a 
user-developed  finite-element  code;  rather,  it  can  be  used  with  commercially  available  codes 
such  as  ABAQUS. 

This  technique  was  applied  to  the  modelling  of  stable  crack  growth  in  a  compact  tension 
specimen  of  A533B  steel.  The  J  integral  and  the  crack-tip  opening  angle  were  used  as  a  two- 
parameter  fracture  criterion.  Numerical  results  agree  well  with  corresponding  experiments 
previously  performed  by  Shih  et  al. 
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Parametric  Representation  of  the 
Yield  Strength-Strain  Rate 
Behavior  of  A533B  Steel 

A  nonlinear  programming  technique  to  estimate  the  parameters  for  a  temperature 
normalized  version  of  the  Malvern  equation  is  presented.  The  general  conjugate 
gradient  algorithm  is  used  to  minimize  a  least-squares  function  formulated  tram 
data  for  A533B  steel.  A  constrained  solution  (which  provides  preferential  treat¬ 
ment  for  a  particular  data  point)  is  also  given.  The  results  satisfactorily  represent 
experimental  data  over  a  temperature  range  from  -  73C  to  260C.  The  approach  v 
well  suited  to  design  applications  where  data  at  a  specific  temperature  of  interest  arc 
limited  or  nonexistent. 


Introduction 

The  yield  strength  of  a  metal  is  a  function  of  a  number  of 
variables,  notably  temperature  and  strain  rate.  For  body- 
centered-cubic  metals,  such  as  steel  above  the  ductile  to  brittle 
transition  temperature,  the  yield  strength  typically  increases 
with  increasing  strain  rate.  Large  increases  in  yield  strength 
are  observed  at  high  strain  rates,  which  may  occur  in  ap¬ 
plications  involving  rapidly  propagating  cracks,  or  projectiles 
impacting  on  a  surface.  In  order  to  take  advantage  of  the 
enhanced  material  properties  at  high  strain  rates,  the  strain 
rate  dependence  of  yield  strength  must  be  characterized 
quantitatively.  With  recent  advances  in  computational  stress 
analysis  techniques  (such  as  the  finite  element  method), 
analyses  employing  strain  rate  dependent  yield  strength  have 
become  tractable,  thereby  emphasizing  the  need  for  an  ac¬ 
curate  yield  strength/strain  rate  relationship. 

A  number  of  stress  strain  strain  rate  temperature 
relationships  have  been  proposed  [1],  One  of  these  is  the 
Malvern  relationship  (1,2)  given  in  equation  (1 ). 

Got  =  al/E  +  D[(a/o„  )  -  I)''  (1) 

If  the  elastic  component  is  eliminated  from  ( 1 ),  we  obtain 

r  =  D[(o/o„  )  -  l]'1  (2) 

where  <  =  the  plastic  strain  rate 

a  =  the  instantaneous  yield  strength 
o„  =  the  "static"  yield  strength  (at  <=  0) 

D.p  -  temperature-dependent  parameters. 

It  is  common  to  estimate  a„,  D.  and  p  by  writing  equation  (2) 
for  three  selected  data  points  and  solving  the  resulting  system 
of  nonlinear  equations.  Such  approaches  are  often  adequate; 
however,  which  three  data  points  to  choose  is  a  subjective 
decision  and  the  quality  of  the  outcome  is  dependent  on  the 
experience  and  insight  of  the  engineer. 

In  this  paper  we  present  an  alternative  method  for  the 


estimation  of  these  parameters  based  on  the  well-known 
principle  of  least  squares.  A  nonlinear  programming  strategy 
is  used  to  accommodate  the  nonlinear  nature  of  the  Malvern 
equation.  Results  are  given  for  a  temperature  normalized 
version  of  equation  (2)  based  on  data  for  A533B  steel.  These 
are  now  being  used  in  FEM  analyses  of  dynamic  crack-tip 
stress  fields  at  Vanderbilt  University . 

Our  purpose  here  is  to  present  an  objective  approach  to  this 
curve-fitting  problem.  The  method  reported  could  be  used  in 
conjunction  with  functional  relations  different  from  equation 
(I).  The  Malvern  model  was  chosen  for  illustration  of  the 
technique  because  it  is  often  included  in  commercially 
available  FEM  codes  and  it  appears  to  be  widely  used.  Ii  is 
importani  to  note  that  the  use  of  a  least-squares  approach  for 
the  estimation  of  the  model  parameters  assures  that  the  model 
will  fully  utilize  the  available  experimental  data.  Presumably, 
this  will  result  in  the  best  possible  curve  fit  within  the 
limitations  of  the  mathematical  form.  Since  there  is  virtually 
always  scatter  in  experimental  data,  any  mathematical 
description  of  the  physical  phenomenon  should  benefit  by 
such  treatment. 

There  are  several  solution  techniques  for  the  nonlinear 
least-squares  problem,  among  them  the  well-known 
Marquardt  algorithm  (3)  and  the  Gauss  least-squares  method 
(4).  These  require  the  formulation  and  inversion  of  a  laigc 
Jacobian  matrix.  The  general  conjugate  gradient  method 
utilized  in  this  note  involves  only  vector  operations.  Since  no 
matrices  or  inversion  processes  are  required,  stability 
problems  are  avoided  and  round  off  is  less  significant.  These 
benefits  are  usually  gained  at  the  expense  of  additional 
iterations  for  convergence;  however,  since  each  iteration  is 
less  complex,  computation  times  are  small  (typically  a  lew 
seconds  or  less  on  a  main  frame  computer) 
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A  considerable  amount  of  yield  strength  strain  rate  data  is 
available  for  A533B  steel  (a  pressure  vessel  qualm  steel  used 
extensively  in  the  nuclear  power  industry  I.  Experimental  data 
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Fig.  1  Yield  strength/strsin  rate  relationship  lor  A533B  at  various 
temperatures.  Strain  rate  is  in  s  "  1 . 


have  been  given  by  Steichen  and  William''  [5]  and  by  Oldfield, 
et  al.  |6).  Their  results  (obtained  from  high-speed  tensile  tests) 
are  plotted  in  Fig.  I .  Server  (’)  has  obtained  similar  data  using 
prepacked  Charpy  specimens. 

Figure  I  show-  that  the  yield  strength  of  A533B  is  tem¬ 
perature  dependent.  An  exponential  dependence  of  strain  rate 
on  temperature  is  predicted  bv  chemical  rate  theory. 
Statistical  mechanics  has  been  used  to  formalize  this 
relationship,  as  discussed  in  (8).  The  relationship,  based  on 
the  Arrhenius  rate  equation  is 

t  «  4e\p[  -  Hi  71  (3) 

where 


Estimation  of  the  Parameters 


•I  =  "frequency"  factor 
H  -  stress-dependent  activation  energy  term 
T  =  absolute  temperature. 

The  activation  energy  term  has  been  shown  to  be  a  t unction  ol 
i he  yield  stress  by  Christian  [9];  hence 

o=/(  71n|.4 /f'|)  (9) 

The  value  of  .4  has  been  experimentally  determined  by 


The  general  curve-fitting  problem  requires  that  we  find  a 
functional  relation.  c(y),  to  represent  n  data  points,  (y  .  c  )■ 
The  least-squares  curve  minimizes  the  sum  of  the  deviations 
squared,  i.e., 

minimize  -  [c(.i, ) -c,);  (6) 

If  ;( i )  is  posed  as  a  polynomial  or  some  other  simple  func¬ 
tion.  then  a  solution  can  be  found  by  solving  a  set  of 
simultaneous  linear  equations.  This  technique  is  discussed  in 
numerous  textbooks  on  statistics  and  numerical  methods  [11. 
12).  If  a  more  complex  relationship  exists  between  ;  and  i, 
then  some  other  approach  is  necessary.  Special  purpose  [3] 
and  general  purpose  algorithms  [13)  have  been  reported  in  the 
literature. 


Bennett  and  Sinclair  [I0[  as  1 0s s  '.  and  is  reportedly  quite 
insensitive  to  stress,  microstructure,  temperature  and  strain 
rate.  Combining  (3)  and  (4)  gives 

ilt=A(i.lA\',>  (?) 

which  converts  strain  rates,  f  ,  at  temperature,  T  ,  to  com¬ 
parable  strain  rates.  i„,  at  (he  reference  temperature,  T„. 
Accordingly,  data  at  any  temperature  can  be  converted  to 
data  at  any  desired  (reference)  temperature. 

The  data  of  Fig.  1  were  obtained  at  three  different  tem¬ 
peratures.  Equation  (?)  was  used  to  transform  these  data  to 
equivalent  points  at  TM  -  366. 15  K  (930.  Of  course,  it  would 
be  possible  to  consider  the  behavior  at  each  temperature 
.eparately  The  principal  disadvantage  ol  this  approach  is  that 
it  requires  experimental  data  at  any  temperature  of  interest. 
Such  data  are  not  always  readily  available  or  easily  obtained. 
The  approach  we  present  here  is  well  suited  to  design  ap¬ 
plications  where  data  at  the  specilic  temperature  of  interest 
are  limned  or  nonexistant  Within  the  range  ot  temperatures 
and  strain  rates  considered  in  Fig  I .  reasonable  results  could 
be  expected  Extrapolation  to  obtain  information  about 
behavior  outside  the  range  ol  the  experimental  data  would  be 
difficult  to  detend  and  is  not  rcsommcnded 


Development  of  the  Objective  Function.  A  relationship 
for  the  yield  strength  as  a  function  of  the  plastic  strain  rate 
can  be  derived  from  (2)  as 

o  =  o  ,[((/£>) ' 1  "'  +  1)  C) 

The  least-squares  formulation  is  then 

minimize/! a,,, D ,p)  =  -(  o„ {(«,/£>) ' 1  ’’ '  <■  H  -  a  I*  (8) 

where  the  function/is  the  objective  function.  V  alues  of  a  . .  D. 
and  p  which  minimize/will  be  obtained  by  the  iterative  search 
strategy  outlined  in  the  forthcoming. 


Algorithm  and  Implementation.  The  objective  I  unction 
was  minimized  by  the  General  Conjugate  Gradient  (GCG) 
algorithm  with  an  acceptable  point  line  search,  as  outlined  in 
[14],  The  variables  a  ..  D,  and  p  were  scaled  by  1  O' .  5  x  I0J. 
and  5.0,  respectively.  Analytical  derivatives  were  used.  Reset 
was  employed  every  ten  multivariate  iterations.  Computations 
were  done  in  double  precision  Fortran  on  the  DEC- 1099 
computer  at  Vanderbilt  University.  Convergence  was  assumed 
when  the  m3enitude  of  the  izradient  was  reduced  to  10  '  or 


The  computer  program  was  used  to  minimize  a  wide  range 
of  functions  with  known  minima  (including  least-squares 
objectives  in  the  lorm  ot  equation  (8))  prior  to  conducting  this 
investigation  The  expected  mimmi/ers  were  correctly 
idem 1 1 icd  bv  the  algorithm  in  everv  case 
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Plastic  Strom  Rate,  Logic) 

Fig.  2  Yield  strength/striin  rite  date  for  AS33B  normelized  to  B3C. 
Strain  rate  is  in  a  ” 1 . 


Results 

The  data  normalized  to  93C  by  equation  (5)  (see  Fig.  2) 
were  used  to  formulate  the  objective  function.  Initial 
estimates  of  the  parameters  were 


' 

100  MPa 

D 

= 

50000 

J>  . 

.5.0 

The  final  solution  obtained  by  the  minimization  algorithm 
was 


Oo  ~ 

'437.40  MPa’] 

D 

= 

297930 

J>  - 

5.6338  J 

rates  in  the  development  of  the  model.  One  approach  would 
be  to  include  the  locally  obtained  data  in  the  least-squares 
formulation  of  equation  (8).  This  would  essentially  gi\e  equal 
consideration  to  all  of  the  available  data  in  the  estimation  of 
the  model  parameters.  A  second  approach,  which  favors  the 
local  data  point  over  the  others,  would  be  to  force  the 
mathematical  function  to  contain  the  preferred  point.  This 
would  be  done  by  taking  one  variable  out  of  equation  (7)  w  ith 
considerable  algebraic  manipulation.  Alternatively  an 
equality  constraint  can  be  introduced  in  the  minimization  to 
enforce  the  preferred  relation.  This  procedure  is  illustrated  tn 
the  forthcoming. 

Consider  the  case  where  the  yield  strength  for  a  particular 
specimen  has  been  found  to  be  41 5  MPa  at  a  plastic  strain  rate 
of  10  '  s ' 1  and  a  temperature  of  93C.  The  required  equably 
constraint  for  the  minimization  would  be 

0o[(O.OOl/D)"  ?'  +  1)-415  x  10‘  =0  (II) 


The  resulting  curve  is  given  by  the  solid  line  of  Fig.  2  where  we 
see  that  equation  (7)  maps  the  data  with  maximum  yield 
strength  error  of  less  than  ±  9  percent  and  an  average  error  of 
less  than  4  percent.  The  solid  line  of  Fig.  2  can  be  transformed 
to  other  temperatures  by  manipulating  equations  (5)  and  (7) 
to  obtain 

a  =  o0[Uj/D' )"  p  '  +  1)  (10) 

where  D'  =  A{D/A)T*'h  and  p'  =  TnplTr  Since  A,  D,p. 
T,  and  T„  are  all  known,  the  yield  strength  can  be  determined 
at  any  strain  rate  ir  corresponding  to  Tr  The  curves 
presented  in  Fig.  1  were  obtained  in  this  manner.  Even  though 
these  curves  are  not  formally  “optimal”  in  the  least-squares 
sense,  there  is  excellent  qualitative  agreement  at  each  tem¬ 
perature  and  deviations  are  acceptably  small.  Again,  ex¬ 
trapolation  outside  the  range  of  existing  data  is  not  recom¬ 
mended. 


The  constrained  nonlinear  programming  problem  i'  given 
by  equation  (8)  such  that  equation  (II)  is  satisfied  and  thi* 
problem  can  be  addressed  in  numerous  ways  (e  g.,  reduced 
gradient,  Lagrange  multiplier  based  penalty  function,  etc..  sec 
[13]  for  a  survey  of  methods).  For  simplicity  we  chose  to  use  a 
sequential  unconstrained  minimization  technique  as  outlined 
by  Fiacco  and  McCormick  [15].  The  initial  constraint 
weighting  factor  was  .1  and  this  was  multiplied  by  0.08  for 
each  successive  minimization.  The  solution  obtained  was 


i  "1  1 

D 

JP  . 

62920  | 

7.2455 


Convergence  was  based  on  step  length  reduction  to  10  * 


Forcing  the  Curve  to  Contain  a  Particular  Point 

The  strategy  outlined  in  the  foregoing  will  be  satisfactory 
for  many  applications;  however,  there  may  be  occasions  when 
it  is  desirable  to  give  preferential  treatment  to  a  particular 
point.  For  example,  facilities  may  be  available  to  conduct 
tensile  tests  at  the  single  strain  rate  of  10  's  '.  It  may  be 
desirable  to  take  advantage  of  published  data  at  higher  strain 


rather  than  the  gradient  test. 

The  constrained  relation  is  given  by  the  dashed  line  of  Fig 
2.  This  curve  minimizes  the  objective  function  of  equation  (8) 
subject  to  the  requirement  that  equation  (ID  be  satisfied 
Although  not  formally  least-squares  optima)  because  it  gives 
preferential  treatment  to  a  single  data  point,  there  are  oc¬ 
casions  when  such  an  approach  may  be  considered  ap¬ 
propriate.  The  nonlinear  programming  method  outlined  here 
can  readily  accommodate  the  situation. 
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Conclusions 

In  this  note  we  have  presented  an  objective  technique  (based 
on  the  principle  of  least  squares)  to  estimate  the  parameters  of 
the  Malvern  equation.  The  nonlinear  programming  approach 
provides  values  for  <r„,  D,  and  p  which  exhibit  excellent 
agreement  with  experimental  data  for  A533B  steel  over  a  wide 
range  of  temperatures. 
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ANALYSIS  OF  CRACK  ARREST  UNDER  ELASTIC-PLASTIC  CONDITIONS 


by 

G.  T.  Hahn*,  P.  Bastias**,  A.  Kumar**,  and  C.  A.  Rubin*** 

ABSTRACT 

Elastic-plastic  finite  element  analyses  of  quasistatic  crack  propagation  and 
arrest  in  a  compact  test  specimen  are  described.  The  analyses  evaluate  the 
effect  of  the  prior  plastic  deformation  history  of  the  propagation  event  on 
the  shape  of  the  J-resistance  curve  of  the  arrestor.  In  addition,  experiments 
employing  a  special  device  are  used  to  simulate  fast  fracture  and  arrest  in 
small  laboratory  test  specimens  of  ductile  HY-80  steel  and  tough  7075  aluminum. 
The  results  show  that  large  dynamic  elevations  of  the  J-curve  are  encountered 
during  arrest.  A  relatively  simple  procedure  for  assessing  or  designing  welded 
in,  steel  crack  arrestors,  based  on  these  findings  is  proposed. 

Key  Words:  Fracture,  crack  arrest,  J-resistance  curve,  HY-80  steel, 

7075  aluminum,  crack  arrestors,  design,  history,  rate,  dynamic 
J-curve. 

INTRODUCTION 

There  is  long-standing  interest  in  designs  that  will  stop  a  rapidly 
advancing  crack  before  structural  integrity  is  compromised^"4^ .  Methods  are 
available  for  determining  the  NDT-temperature^3)  or  the  arrest  toughness 
values^-8)  0f  steel  plates  that  will  arrest  a  short  crack,  i.e.,  ~  10  mm- 
long,  emerging  from  a  locally  embrittled  region.  Existing  toughness 

^Professor,  **Graduate  Student,  and  ***Associate  Professor,  Department  of 
Mechanical  and  Materials  Engineering,  Vanderbilt  University,  Nashville,  TN 
37235 


specifications  have  thus  controlled  the  kind  of  "pop-in"  crack  that  originated 
from  an  arc  strike  and  severed  the  hull  of  the  U.S.S.  Ponaganset  in  1947 
However,  there  is  no  general  and  rational  method  for  designing  welded-in 
crack  arrestors  that  can  stop  a  long  fracture,  i.e.,  0.1m  <  o  <  2m  in  large 
structures  such  as  ship  hulls,  pressure  vessels,  line  pipe,  etc. (9)  This 
remains  a  difficult  problem  because  very  high  toughness  levels:  200  MPa/m  < 

KIa  <  600  MPa^m  (200  KJ/m^  <  Jja  <  2  MJ/m*)  are  needed  to  arrest  long  fractures 
in  a  fully  loaded  structure.  Figure  1  describes  a  practical  example:  a 
fracture  initiated  by  a  short  defect  in  a  100  mm-wide,  brittle  weld.  The 
adjacent  steel  plates  probably  must  possess  an  arrest  toughness,  Kja  *  200 
MPa 4  to  stop  the  unstable  crack  from  penetrating  with  a  reasonable  margin  of 
safety.*  If  the  plates  fail  to  stop  the  fracture,  a  tough  welded-in  arrestor 
of  the  same  weakness  located  in  the  path  of  the  crack  1  m  from  the  weld  may 
produce  arrest  if  it  possesses  a  toughness  Kja  *  600  MPa.**  Arrestors 
located  farther  from  the  weld  would  have  to  be  proportionately  thicker  or 
tougher. 

The  LEFM-based  toughness  estimates  given  above  are  not  reliable  because 
the  LEFM  requirements  are  difficult  to  satisfy  at  the  high  Kj-level s.*** 


*The  statically  evaluated  Kj-value  for  the  2a  =  100  mm-long  crack  emerging 
from  the  weld  zone  is  Kj  =  149  MPa/m.  This  assumes  a  total  stress  of  o=  375  MPa 
(Applied  =  250  MPa  plus  Residual  =  125  MPa),  and  that  the  other  dimensions 
are  much  larger  than  2a.  This  calls  for  an  arrest  toughness  of  at  least 
Kja  =  200  MPa^m  to  arrest  the  crack  with  a  margin  of  safety. 

**The  statically  evaluated  Kj-value  for  the  2a  =  2m-long  crack  entering  a 
welded  in  arrestor  of  the  same  thickness  is  Kj  =  443  MPa.  This  assumes  that 
the  weld  residual  stresses  can  now  be  neglected  and  that  the  other  dimensions 
of  the  structure  are  still  much  larger  than  2a.  An  arrest  toughness  of  at 
least  Kja  *  600  MPa/m  is  required  to  stop  the  crack  with  a  margin  of  safety. 

***For  example,  to  satisfy  the  requirement  that  the  crack  length  is  less  than 
1/10  of  the  plastic  zone  extent:  a  >  10  rQ  ■  2  (Kj/cl)  .  When  Kj  =  200  MPa/m 
and  o0  =  500  MPa,  this  means  that  2a  >  0.6m  at  arrest  which  is  6x  the  value 
involved  in  the  preceding  example. 


Figure  1.  Example  of  a  practical  crack  arrest  problem.  A  100  mm-wide  weld: 
fusion  plus  HAZ,  with  a  low  toughness:  Kjq  =  40  MPa  (Jjc  =  7.7  KJ/nr)  contains 
a  2a  =  4mm  long  crack.  The  crack  becomes  unstable  when  the  total  stress 
(°applied  s  250  MPa  Plus  °residual  *  250  MPa)  reaches  o  =  500  MPa.  The 
crack  begins  to  extend  and  becomes  a  2a  =  100  mm-long  rapidly  propagating 
fracture  when  it  emerges  from  the  weld  zone  and  begins  to  penetrate  the 
adjacent  steel  plates.  Question:  How  tough  do  the  plates  have  to  be  to  assure 
that  the  crack  will  be  arrested?  How  tough  does  a  welded-in  arrestor  plate 
located  in  the  path  of  the  fracture  1  m  from  the  weld  have  to  be  to  assure 
arrest? 
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They  are  also  unreliable  because  the  estimates  given  do  not  take  account  of 
the  increase  in  toughness  obtained  when  the  ductile  crack  begins  to  penetrate 
the  arrestor  as  described  by  its  R-curve.  The  J-integral  concept  together 
with  the  J-resistance  curve  are  therefore  better  suited  for  characterizing 
the  response  of  relatively  tough  materials.  With  this  in  mind,  Hahn  and 
Dantam^10^  proposed  that  the  JR-value  evaluated  at  to  =  1.5  mm  can  serve  as  a 
measure  of  the  crack  arrest  toughness:  Ja(2)  =  JR.*  This  measure  is  valid 
for  long  cracks  and  arrestors  that  display  the  ductile  fibrous  mode  of  fracture 
during  the  arrest  process.** 

However,  the  use  of  J-curves  in  arrest  analyses  involves  several  unresolved 
complications  illustrated  in  Figure  2: 

( i )  Deformation  History  Effects.  Figures  2a,  b  and  c  illustrate  that  the 
shape  of  the  J-resistance  curve  depends  on  the  prior  plastic  deformation 
history.  The  conventional  J-curve  (Figure  la),  produced  by  ductile  crack 
extension  with  a  prior  history  of  completely  brittle  extension,  is  different 
from  the  J-curve  (Figure  lb)  for  ductile  crack  extension  with  a  prior  history 
of  ductile  crack  extension.  The  J-curve  for  a  propagation-arrest  event  (in 
Figure  lc)  is  intermediate.  This  raises  questions  whether  the  conventional 
JR-curve  has  the  right  shape  for  a  crack  arrest  analysis. 

(ii)  Geometry  Dependence.  The  shape  of  the  J-curve  for  a  tough  material  is 
geometry  dependent  for  anything  but  small  amounts  of  ductile  crack  extension, 

*Hahn  and  Dantam  identify  two  measures  of  arrest  toughness:  Ja(i)  =  Jic,  and 
Ja(2)  =  JR(Aa  =  1.5  mm).  The  quantity  Ja(l)  is  appropriate  when  the  crack 
extension  is  a  small  fraction  of  the  initial  crack  length;  Ja(2)  when  the 
extension  is  a  large  multiple  of  the  initial  length.  The  Kj -equivalents  of 
these  2  values  are:  K ja ( 1 )  =  Kjc  =  I CE *  and  Kja(2)  =  KJR  =  "UrE* 

**The  transition  temperature  of  a  steel  arrestor  must  be  well  below  the  service 
temperature  to  assure  that  it  will  display  high  toughness  fibrous  fracture 
under  dynamic  conditions  rather  than  low  toughness  cleavage. 
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Figure  2.  Complications  attending  the  use  of  J-resistance  curves  in  the 
analysis  of  crack  arrest:  (a),  (b)  and  (c)  compare  the  shape  of  the  J-curve 
for  ductile  crack  extension  from  a  pre-existing  crack  of  length  a2  with 
different  prior  deformation  histories  in  the  interval  from  aj  to  a2, 

(d)  schematic  drawing  illustrating  the  excess  strain  energy  (shaded  area  Aj) 
part  of  which  may  be  returned  to  make  up  the  energy  deficit  (shaded  area  A2) 
and  permit  the  crack  to  split  the  arrestor,  and  (e)  example  of  dynamic 
elevation  of  the  J-curve. U7.18) 


typically  no  more  than  0.05  <  Aa/(w-a)  <  0.1  where  w  is  the  uncracked  ligament. t11) 
Consequently,  it  is  again  not  clear  that  the  conventional  J-curve  is  appropriate 
when  arrest  is  preceded  by  large  amounts  of  extension,  as  in  the  examples 
cited  earlier. 


(iii)  Kinetic  Energy  Release.  A  part  of  the  excess  strain  energy  (defined  by 
the  difference  between  the  static  J-driving  force  and  J-resi stance  curves  at 
instability  in  Figure  2d)  is  converted  into  kinetic  energy,  and  part  of  this  may 
be  returned  adding  to  the  crack  driving  force  during  arrest  (see  Figure  2).(fi»12) 
The  kinetic  energy  return  is  difficult  to  evaluate  without  a  dynamic  analysis 
and  this  complicates  the  task  of  defining  the  arrest  condition. 

(iv)  Dynamic  Elevation  of  the  J-Curve.  There  is  evidence  from  several  sources 
that  the  J-curve  ductile  crack  extension,  normally  measured  for  very  low  rates 
of  crack  extension,  i.e.,  a  *  10"4  ms”*  is  elevated  by  increasing  the  crack 
speed. (13-15)  This  does  not  fundamentally  alter  the  ultimate  crack  stopping 
ability  of  the  material  because  the  crack  velocity  must  be  reduced  to  zero  even¬ 
tually  to  produce  arrest.  However,  the  transient  elevation  of  the  J-curve 
produced  by  the  high  velocities  of  unstable  cracks  102ms_1  <  a  <  103ms'1,  could 
provide  a  mechanism  for  dissipating  the  returned  kinetic  energy. 

This  paper  presents  finite  element  calculations  of  the  effect  of  the  prior 
deformation  history  on  the  shape  of  the  J-curve.  It  also  presents  measurements 
of  the  arrest  of  fully  brittle  fractures  in  tough,  ductile  aluminum  and  HY-80 
steel.  The  results  show  that  prior  history  can  frequently  be  neglected  and 
that  the  J-curve  retains  its  geometry  independence  when  arrest  is  preceded  by 
relatively  brittle  fracture.  There  are  indications  the  elevation  of  the 
J-curve  with  crack  velocity  can  compensate  for  kinetic  energy  return.  Simplified, 
guidelines  for  designing  crack  arrestors  are  proposed  on  this  basis. 


ANALYSES  AND  PROCEDURES 


Finite  Element  Analysis 

Finite  element  calculations  were  performed  to  examine  the  effects  of  the 
deformation  history  -  -  the  plastic  deformation  produced  by  an  advancing 


fracture  before  it  reaches  the  arrestor  -  -  on  the  J-resistance  curve 


subsequently  displayed  by  the  arrestor.  To  simplify  the  problem,  quasi¬ 
static,  stable  crack  growth  and  arrest  were  simulated  in  the  elastic-plastic 


finite  element  model  of  a  4T  Compact  Test  Specimen  shown  in  Figure  3.  The 
model  is  composed  of  172  quadrilateral  and  8  noded  elements  with  a  total  of 


38  nodes  forming  the  path  of  the  crack.  The  simulations  were  performed  with 
the  finite  element  code  ABAQUS  on  VAX  11/780  and  11/785  machines. 


The  initial  position  of  the  crack  tip  is  within  the  region  labeled 
Material  A,  which  is  endowed  with  a  relatively  low  toughness.  The  model  is 


subjected  to  a  continuously  increasing  load-point  displacement  and  the  nodes 
ahead  of  the  crack  are  released  at  the  appropriate  time  to  simulate  crack  exten 
sion  on  the  symmetry  plane.  Details  of  the  procedure  are  given  el sewhere. 


The  crack  eventually  reaches  the  region  labeled  Material  B  after  a  crack 
extension  of  either  Aa  =  6mm  (shown  in  Figure  3)  or  4a  =  2mm  (configuration 
not  shown).  The  crack  then  arrests  for  a  time  since  Material  B  is  tougher 
than  Material  A,  and  then  continues  to  grow  as  the  conditions  for  extension 
in  Material  B  are  satisfied  as  a  result  of  the  increasing  load  point  displacement. 

The  deformation  behavior  of  the  2  materials  is  the  same.  It  is  given  by 
the  piecewise-1 inear  representation  of  the  stress-strain  curve  of  A533B 
pressure  vessel  steel  at  93° C  (see  Figure  4a),  and  approximates  power  law 
hardening  after  the  initial  yield:  E  =  198  GPa,  oQ  =  383  MPa,  n  =  10. (13,17) 

The  crack  extension  in  Material  A  was  controlled  to  proceed  according  to  a 
predetermined  J-resistance  curve.  The  effects  of  several  J-curves  appropriate 
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Figure  3.  Finite  element  model  of  4T  compact  test  specimen. 
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Figure  4.  Stress-strain  characteristics  and  J-resistance  curve  of  model 
materials:  (a)  piecewise-1 inear  stress  strain  curve  of  A533B  pressure  vessel 
steel  employed  for  Material  A  and  Material  B,  and  (b)  conventional  J-resistance 
curve  of  Material  B.U7.18) 


for  toughness  levels  in  the  range:  50  MPa^m  <  Jc  <  120  MPa/m  were  examined. 

Crack  extension  in  Material  B  was  controlled  by  the  crack  opening  displacement- 
resistance  curve  previously  derived  for  A533B  at  93*C^3,18):  COD  =  2.5332  •  10-4 
+  0.47082  4a  -  33.58  (4a)2.  This  resistance  curve  is  viewed  as  a  local  criterion 
that  is  insensitive  to  the  prior  deformation  history.  For  ductile  crack  exten¬ 
sion,  it  produced  a  toughness  of:  Jjc  *  200  KJ/m2  (KJc  =  203  MPa^m),  and  the 
conventional  J-curve  shown  in  Figure  4b,  when  the  extension  proceeds  from  a 
pre-existing,  perfectly  brittle  crack.  In  the  present  calculations,  the 
different  deformation  histories  prior  to  arrest  produced  different  J-Aa 
relations  for  the  arrestor.  These  J-curves  were  evaluated  with  the  finite 
element  model  as  the  crack  penetrated  into  Material  B. 

Crack  Arrest  Measurements 

Measurements  of  crack  arrest  in  a  ductile  steel  and  aluminum  alloy  were 
carried  out  to  evaluate  the  dynamic  elevation  of  the  J-curve  by  large  crack 
velocities.  These  measurements  employed  a  non-standard,  6  in  x  6  in  (152.4  mm 
x  152.4  mm)  pre-cracked,  compact-type  test  specimen  whose  measuring  capacity 
is  increased  by  attaching  reusable,  hardened  steel  arms.  The  operation  of  the 
system,  referred  to  here  as  a  "fast  fracture  device"  is  shown  schematically  in 
Figures  5a  and  5b.  The  test  specimen  and  other  components  are  described  in 
Figures  5c  and  5d.  To  facilitate  the  analysis  of  the  experiments,  the  variation 
of  compliance  of  the  system  with  crack  length  was  evaluated  with  the  aid  of  a 
finite  element  model.  Details  of  these  calculations  are  presented  elsewhere. (19) 

After  the  bolt  breaks,  the  system  behaves  as  if  split  by  a  zero  toughness 
crack  which  runs  to  the  precrack  and  is  arrested  by  the  compact  specimen,  i.e., 
ductile  arrest  with  a  history  of  completely  brittle  crack  propagation.  The 
displacement  of  the  arms  at  the  bolt  and  wedge  at  the  onset  of  fracture  were 


no 
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measured  with  standard  displacement  gages.  These  displacements  and  the  static 
compliance  values  for  the  system  before  and  after  arrest  were  used  to  estimate 
the  elastic  energy  released  during  the  arrest  event.  The  average  value  of 
the  propagating  crack  toughness,  Gg,  was  then  determined  by  dividing  the 
released  energy  by  the  total  fracture  area,  assuming  that  all  of  the  released 
energy  Is  consumed  as  fracture  energy.  The  quantity  Gq  can  be  viewed  as  a 
measure  of  the  elevation  of  the  J-resI stance  curve  as  a  result  of  crack 
extension  and  crack  velocity. 

Measurements  were  carried  out  on  HY-80  steel  and  the  7075-Aluminum 
alloy,  which  was  solution  treated  at  470*C  for  1/2  hr.  and  aged  at  120°C  at 
different  time  Intervals  to  achieve  HRB  *  45,  78,  and  92.  The  toughness 
values  Jjc  and  Ja( 2 )  *  Jr  are  listed  in  Table  1. 

RESULTS 

Finite  Element  Analyses 

The  results  of  the  finite  element  computations  are  presented  in  Figure  6. 
These  show  the  different  J-curves  which  were  prescribed  for  Material  A 
and  controlled  the  "start"  of  the  relatively  brittle  propagation,  and  the 
corresponding,  calculated  J-curves  for  Material  B  which  served  as  the  arrestor. 
The  results  illustrate  the  J-curve  for  the  arrestor  is  influenced  by  the 
deformation  history,  and  is  unaffected  by  the  deformation  history  when  the 
starter  toughness,  i.e.,  the  Jic-values  for  the  material  in  which  the  crack 
starts  propagating,  are  much  lower  than  the  arrest  toughness: 


Jic  (starter) 

-  <  0.1  or 

(la) 

Jjc  (arrestor) 

K jc  (starter) 

-  <  0.3 

(lb) 

Kjc  (arrestor) 
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Figure  6.  Results  of  finite  element  calculations  showing  the  effects  of  the 
prior  plastic  deformation  history  on  the  shape  of  the  J-resistance  curve  of 
the  arrestor. 


Figure  7.  Summary  of  finite  element  calculations  showing  the  effects  of  the 
toughness  of  Material  A  (starter  toughness)  on  the  calculated  toughness  of 
Material  B  (arrestor  toughness)  for  crack  extension  Increments  of  Aa  =  2mm 
and  Aa  *  6mm. 


These  results  are  summarized  in  Figure  7,  which  also  shows  that  these  conclusions 
are  relatively  insensitive  to  the  extent  of  propagation  prior  to  arrest.  In 
other  words,  it  appears  that  the  arrestor's  response  to  a  relatively  brittle, 
prior  propagation  history  is  the  same  as  its  response  to  a  pre-existing  crack. 

It  follows  that  the  rules  for  the  geometry  independence  of  the  J-curve  for 
extension,  from  a  pre-existing  (brittle)  crack,  will  also  apply.  This  means 
that  the  conventional  J-curve  for  the  arrestor  material  will  describe  resistance 
to  penetration  even  when  arrest  is  preceded  by  large  crack  extensions  and  for 
relatively  small  remaining  ligaments. 

Crack  Arrest  Measurements 

The  results  in  Table  1  show  that  the  tough  steel  and  aluminum  alloy  suffer 
relatively  small  crack  penetrations  prior  to  arrest.  When  the  relatively 
large  energy  releases,  e.g.  ~  400  J,  produced  by  the  fast  fracture  device  are 
combined  with  small  penetrations  and  deep  face  grooves,  large  values  of 
Gq,  the  fast  propagating  crack  toughness,  result.  The  Gp  value  is  18X  and  3.6X 
larger  than  Ja(2)  for  the  HY-80  steel  and  the  tough,  HRB  =  45,  aluminum  alloy, 
respectively.  These  increases  are  qualitatively  consistent  with  the  measurements 
of  Joyce^4»15)  -jn  Figure  2e  which  reveal  a  ~  2X  increase  in  J-values  for  an 
increase  in  crack  velocity  from  a  ~  10”4  ms"1  to  a  ~  0.25  ms”1.  These  increases 
in  toughness  are  related  to  the  rate  sensitivity  of  the  plastic  flow  stress^13) 
and  to  changes  in  the  fracture  process^).* 

The  interpretation  of  the  present  measurements  involves  uncertainties 
because  other  sources  of  energy  dissipation  have  been  neglected.  However, 
even  allowing  for  a  factor  of  2  error  in  the  energy  consumed  by  fracture, 

*The  fracture  surfaces  of  the  HY-80  steel  and  HRB  45,  7075  aluminum  are  ~  5X  and 
~  3X  rougher,  respectively,  than  the  surfaces  produced  by  a  slow  tearing  test. 


Table  1.  Summary  of  Crack  Arrest  Measurements 


7075  A1 


HRB  45 

0.25 

8.5 

192 

740 

2680 

HRB  78 

0.25 

50.8 

52 

292 

670 

HRB  90 

0.25 

>88 

22 

148 

<490 

B^/B  is  the  fraction  of  the  specimen  thickness  remaining  at  the  face  grooves. 


it  appears  that  the  dynamic  elevation  of  the  J-curve  for  arrest  of  fast 
propagating  cracks  by  ductile  steel  could  be  ~  10X.  Such  a  high  transient 
toughness  would  consume  the  kinetic  energy  return  which  is  not  accessible  to 
the  static  analysis. 

DISCUSSION 

The  results  of  this  study  suggest  that  a  relatively  simple  J-integral 
based  procedure  may  be  adequate  when  assessing  or  designing  tough,  welded-in, 
steel  crack  arrestors  for  long  crack  extensions,  i.e.,  Aa  >  10  a0.  The  main 
features  of  such  a  procedure  and  a  tentative  requirement  for  the  arrestor 
width  are  proposed  below: 

(i)  Crack  Driving  Force.  The  crack  driving  force  produced  when  the  crack 
reaches  the  arrestor  interface  is  evaluated  in  terms  of  Jj  using  a  static 
analysis  usch  as  the  one  developed  by  Kumar,  German  and  Shih(2°)  and  Kumar, 
et  al.(21).  when  appropriate,  an  LEFM  analysis  can  be  used. 

( i i )  Crack  Arrest  Toughness.  The  plane  strain,  crack  arrest  toughness: 

Ja(2)  =  Jr  (Aa  =  1.5  mm),  is  derived  from  information  provided  by  a  conventional 
Jic  -determination  (ASTM  E-813)  or  a  Resistance  Curve  measurement  (ASTM  E-561). 
To  assure  ductile,  fibrous  crack  extension  in  the  arrestor  material,  DT 
tests  of  the  arrestor  must  show  1002  fibrous  fractures  at  temperatures  from 
25°C  to  50°C  below  the  lowest  service  temperature.  Effects  of  prior  history 
are  neglected.  Plane  strain  conditions  are  assumed  since  these  are  more 
easily  obtained  with  a  rapidly  propagating  crack,  and  effects  of  the  deformation 
history  on  the  0-curve  of  the  arrestor  are  neglected. 

(Hi )  Arrest  Criterion.  The  criterion  for  arrest  is  Jj  <  Ja(2).  Complications 
attending  kinetic  energy  return  are  neglected  for  steel  arrestors. 
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(iv)  Arrestor  Thickness.  The  arrestor  thickness  determines  the  level  of 
tensile  stress  and  the  value  of  Jj.  Increasing  the  thickness  will  reduce  Oj 
and  improve  the  crack  stopping  capability. 

(v )  Arrestor  Width.  The  arrest  criterion  employed  insures  that  in  the  absence 
of  kinetic  energy  return,  the  maximum  crack  penetration  into  the  arrestor  Aa 

<  1.5  mm.  With  kinetic  energy  return  (i.e.,  total  crack  extensions  prior  to 
arrest  that  are  large  compared  to  the  overall  dimensions  of  the  structure), 
arrestor  penetrations  are  likely  to  be  Aa  <  50  mm.  To  assure  that  the  arrestor 
can  both  survive  such  a  penetration  and  contain  the  plastic  zone  of  the  crack 
to  reduce  the  risk  of  weld-failure,  an  arrestor  width  dimension,  w,  related 
to  the  plastic  zone  size  at  arrest  and  the  Upvalue  of  par  1  is  proposed. 

0.4  JjE 

W  =  0.05m  +  - - —  (2) 

V 

where  E  and  aQ  are  the  elastic  modulus  and  yield  strength  of  the  arrestor. 

A  similar  procedure  for  aluminum  arrestors  may  require  an  added  provision 
for  kinetic  energy  return  since  the  dynamic  elevation  of  the  J-curve  may  be 
too  small  to  accommodate  it.  The  proposed  procedure  is  intended  to  serve  as 
a  starting  point  for  further  discussion.  Clearly,  as  a  design  approach,  it 
requires  careful  scrutiny  and  should  be  critically  tested  before  adoption. 

Values  of  Ja(2)  and  corresponding  estimates  of  the  crack  stopping  capabilities 
of  several  steels  derived  from  the  proposed  procedure  are  listed  in  Table  2. 

CONCLUSIONS 

1.  Effects  of  the  prior  plastic  deformation  history  on  the  shape  of  the  J- 

resi stance  curve  have  been  examined  with  a  finite  element  model.  The  calcu¬ 
lations  reveal  that  history  effects  can  be  neglected  when  the  toughness  rati 
Jjq  ( starter) /J i c  (arrestor)  <  0.1  or  Kjc  (starterJ/Kp  (arrestor)  <  0.3 
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2.  Experiments  employing  a  special  fast  fracture  device  show  that  very 
large  dynamic  elevations  of  the  J-resistance  curve  are  encountered  during 
the  arrest  of  fast  propagating  cracks  in  tough,  ductile  arrestors.  The 
measurements  point  to  an  18-fold  and  3.6-fold  increase  in  the  level  of  the 
J-curve  for  steel  and  an  aluminum  alloy,  respectively. 

3.  The  dynamic  elevation  of  the  J-curve  for  fibrous  crack  extension  in  a 
tough  steel  may  be  large  enough  to  consume  any  kinetic  energy  which 
is  returned  to  the  crack  during  arrest. 

4.  The  findings  suggest  that  a  relatively  simple  J-based  procedure  may  be 
adequate  for  assessing  or  designing  tough,  welded-in  steel  crack 
arrestors. 
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Figure  1.  Example  of  a  practical  crack  arrest  problem.  A  100  mm-wide  weld: 
fusion  plus  HAZ,  with  a  low  toughness:  KjC  =  40  MPa  (Jjc  =  7.7  KJ/m2)  contains 
a  2a  =  4mm  long  crack.  The  crack  becomes  unstable  when  the  total  stress 
(Applied  =  250  MPa  plus  <Vesidual  =  250  MPa)  reaches  a  =  500  MPa.  The 
crack  begins  to  extend  and  becomes  a  2a  =  100  mm-long  rapidly  propagating 
fracture  when  it  emerges  from  the  weld  zone  and  begins  to  penetrate  the 
adjacent  steel  plates.  Question:  How  tough  do  the  plates  have  to  be  to  assure 
that  the  crack  will  be  arrested?  How  tough  does  a  welded-in  arrestor  plate 
located  in  the  path  of  the  fracture  1  m  from  the  weld  have  to  be  to  assure 
arrest? 

Figure  2.  Complications  attending  the  use  of  J-resistance  curves  in  the 
analysis  of  crack  arrest:  (a),  (b)  and  (c)  compare  the  shape  of  the  J-curve 
for  ductile  crack  extension  from  a  pre-existing  crack  of  length  a£  with 
different  prior  deformation  histories  in  the  interval  from  ai  to  a2» 

(d)  schematic  drawing  illustrating  the  excess  strain  energy  (shaded  area  A^ ) 
part  of  which  may  be  returned  to  make  up  the  energy  deficit  (shaded  area  A£ ) 
and  permit  the  crack  to  sol  it  the  arrestor,  and  (e)  example  of  dynamic 
elevation  of  the  J-curve. 117,18) 

Figure  3.  Finite  element  model  of  4T  compact  test  specimen. 

Figure  4.  Stress-strain  characteristics  and  J-resistance  curve  of  model 
materials:  (a)  piecewise-linear  stress  strain  curve  of  A533B  pressure  vessel 
steel  employed  for  Material  A  and  Material  B,  and  (b)  conventional  J-resistance 
curve  of  Material 

Figure  5.  Fast  fracture  device:  (a)  and  (b)  schematic  drawing  of  operation, 

(c)  test  specimen,  and  (d)  details  of  assembly. 

Figure  6.  Results  of  finite  element  calculations  showing  the  effects  of  the 
prior  plastic  deformation  history  on  the  shape  of  the  J-resistance  curve  of 
the  arrestor. 

Figure  7.  Summary  of  finite  element  calculations  showing  the  effects  of  the 
toughness  of  Material  A  (starter  toughness)  on  the  calculated  toughness  of 
Material  B  (arrestor  toughness)  for  crack  extension  increments  of  Aa  =  2mm 
and  Aa  *  6mm. 
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ABSTRACT 

This  paper  reports  a  preliminary  study  of  a  new  laboratory  test  method 
designed  to  measure  crack  arrest  toughness  for  tough,  arrester  grade  steels 
The  design  Is  basically  focused  on  producing  large  enough  driving  force  to 
produce  last  fracture  and  arrest  in  a  small  test  piece  of  a  very  tough, 
ductile  material.  Tests  conducted  on  HY-80  steel  and  7075  aluminum  specimens 
are  reported.  The  analysis  of  arrest  measurements  presented  assumes  that  all 
the  energy  released  Is  conserved  or  converted  Into  fracture  energy.  A  static, 
lower-bound  arrest  toughness  Is  evaluated  from  compliance  analysis  using  an 
elastic-static  finite  element  simulation  of  the  loading  device.  A  photo gramme  trie 
aspect  of  the  fracture  surface  Is  performed  to  quantify  and  relate  the 
micros tructural  features  with  fracture  toughness.  The  results  are  consistent 
with  the  view  that  arrest  toughness  value  derived  from  the  J-Reslstance  curve 
Is  a  valid  measure  of  arrest  capabilities  of  a  tough  material.  It  also 
indicates  that  rapid  crack  extension  In  a  tough,  ductile  material  can  consume 
large  amounts  of  kinetic  energy  released  from  the  process  of  fracture. 
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INTRODUCTION 


Initial  attempts  to  measure  crack  arrest  capabilities  of  steel  by 
Robertson! 1 1 ,  were  followed  by  the  development  of  the  Esso  test!?!,  and  the 
double  tension  test  described  by  Yoshlkl,  et  al.!®].  These  tests  provided 
measures  of  applied  stress,  crack  length  and  temperature  for  which  arrest  is 
possible.  Crosley  and  Ripllng  applied  fracture  mechanics  concepts  to  measure 
the  arrest  toughness,  Kja,  with  a  double  cantilever  beam  (DCB)  specimen^]  . 
This  was  followed  by  the  developed  ot  a  wedge  loaded  'duplex'  specimen  with  a 
hardened  steel  starter  section  by  Hahn,  et  al.f5l,  and  a  similarly  loaded 
specimen  with  a  hard  weld  starter  by  Crosley  and  Rlpllng!®].  Several  methods 
have  been  proposed  for  analyzing  measurements  of  crack  arrest  In  stiff, 
wedge  loaded  specimens,  whose  static  elastic  energy  release  rate  parameter, 

G,  decreases  continuously  with  crack  growth.  Kraft  and  Irwin!7],  and  Crosley 
and  ripllng!®!  proposed  a  static,  nonconservative  interpretation  of  the  crack 
arrest  process  In  which  they  assumed  that  all  the  kinematic  energy  of  tne 
running  crack  Is  dissipated  before  arrest.  Hoagland,  Hahn  and  Rosenf leld!^] 
proposed  an  alternative  static  analysis  in  which  the  klneatlc  energy  was 
assumed  to  be  completely  conserved.  In  this  case  the  same  experiment  Is 
Interpreted  differently.  Accordingly,  in  the  latter  stages  of  the  propagation 
event,  the  stored  kinetic  energy  supplied  part  of  the  fracture  energy 
requirements.  The  same  authors  proposed  a  dynamic  analysis  assuming  the 
conservation  of  kinetic  energy  and  accounting  for  the  dependence  of  the 
propagating  crack's  material  toughness,  Gp  on  crack  velocity!*®].  The  latter 
theory  was  confirmed  by  experiments  conducted  by  Kalthoff,  Beinert  and 
Winkler! **1 . 

The  above  concepts  provide  a  reasonable  estimate  of  the  crack  arrest 
behavior  of  relatively  brittle  materials.  However,  with  very  tough  materials, 


neither  the  apparent  critical  energy  release  rate,  Gc,  nor  Jc  remain  constant 
but  Increase  with  crack  length.  Recent  work  by  Hoff,  Hahn  and  Rublnt^]  on 
strain  rate  effects  using  finite  element  analysis  revealed  that  the  J- 
resl stance  curve  Is  elevated  when  a  crack  extends  rapidly  In  a  material  with 
rate  sensitive  yield  and  flow  properties.  This  result  agrees  with  the 
experimental  evidence  of  Joycet^]^  jn  real  structures,  crack  extentlon  Is 
also  likely  to  proceed  very  rapidly,  producing  large  plastic  strain  rates 
(lO^/s  to  10^/s).  As  a  consequence,  there  will  be  large  elevations  in  the 
yield  stress,  resulting  In  higher  toughness  values  while  the  crack  Is 
propagating.  With  this  In  mind,  Dantam  and  Hahnl^]  have  proposed  two  points 
on  the  J-resistance  curve  (or  their  corresponding  K  equivalents)  as  lower 
bound  measures  of  crack  arrest  toughness  that  can  be  used  to  design  arresters: 

Ja(i)s  JIc,  KJa ( 1 )  -  UlcE/(l-  v2))1/2  (1) 

Ja(2)s  JR(  a  -  1.5  mm),  KJa(2)  -  URE/(1-  v2)]1/2  (2) 

The  quantities  Ja(l)  and  KJa(l)  are  appropriate  tor  short  arrester  cracks, 
l.e.,  a  <  100  mm,  where  the  instability  condition  Is  determined  by  the 
Initial  part  of  the  Jg-curve;  Ja(2)  and  Kja(2)  are  Indicated  tor  longer 
arrested  cracks,  l.e.,  a  >  500  mm,  where  the  instability  condition  Is  determined 
by  the  latter  part  ot  the  JR-curve .  While  these  quantities  are  relatively 
easy  to  obtain  for  tough  materials  from  JR-meaurements ,  their  use  as  arrest 
toughness  parameters  has  not  been  critically  tested. 

The  proposed  ASTM  crack  arrest  method  relies  on  the  non-conservative, 

LEFM  analysis  that  places  lower  limits  on  the  size  of  the  test  piece.  The 
maximum  allowable  Initiation  stress  Intensity,  Kq,  for  a  .200  m  x  .200  m  x 
050  m  specimen  Is  about  120  MPa  m1/0.  Therefore  a  test  specimen  of  the 
design  4x  to  ~  16x  larger  would  be  required  to  evaluate  very  tough  materials 
like  ship  hull  steels  with  arrest  toughness  values  In  the  range  200  MPa  m^/2 


to  500  MPa  ml/2#[15]  in  other  words,  the  ASTM  specimen  size  required  becomes 
prohibitive. 

There  are  several  ways  of  Increasing  the  capacity  of  crack  arrest  tests. 
One  Is  to  insert  deeper  face  grooves  on  the  specimens,  thus  reducing  the 
energy  consumed  by  the  crack,  proportionately.  Although  questions  are  raised 
about  the  3- dime ns tonality  of  the  crack  tip  stress  field  Introduced  by  the 
face  grooves,  experimental  evidence! 16 ,17]  shows  that  face  grooves  occupying 
80%  of  the  cross  section  do  not  interfere  with  static  measurements.  Another 
way  of  Increasing  the  capacity  Is  to  attach  reusable  hardened  steel  arms. 
Initiation  can  be  facilitated  by  an  external  mechanical  restraint  which  Is 
suddenly  removed.  The  concept  that  evolved  in  the  present  work  utilizes  a 
combination  of  the  above  techniques  to  increase  capacity. 

Figure  1  schematically  Illustrates  the  operation  of  the  modified  test. 
Load  is  applied  at  point  D  by  means  of  a  wedge  and  split  pin,  which  offers  a 
8 tiff  loading.  The  hardened  steel  arms  (B)  attached  to  the  pre-cracked 
specimen  (A)  are  wedged  open  against  the  restraint  of  the  notched  bolt  (C). 

As  the  wedge  load  Is  Increased,  the  notched  bolt  breaks  and  releases  the 
elastic  energy  stored  In  the  arms,  causing  the  crack  to  grow  under  a  high 
strain  rate  condition.  Since  stiff  loading  produces  a  decreasing  stress 
intensity,  K,  with  crack  length,  the  crack  arrests  at  a  point  F  when  the 
stress  Intensity  factor  drops  below  the  fracture  resistance  of  the  material 
at  arrest. 

As  part  of  the  present  study,  an  elastic-static  finite  element  analysis  . 
has  been  performed  to  evaluate  the  compliance  of  the  fast  fracture  loading 
device  with  the  bolt  Intact  and  after  It  fractures,  using  the  general  purpose 
flntle  element  code,  ABAQUs! 18 ,19] .  Measurements  of  fast  fracture  and  arrest 
were  also  carried  out  with  the  device  on  7075  aluminum  alloys  and  HY-80  steel. 


Estimates  of  the  propagating  crack  toughness,  and  the  arrest  toughness, 

Kja,  are  derived  from  the  compliance  values  assuming  that  all  of  the  kinetic 
energy  is  conserved  and  is  converted  into  fracture  energy. 

EXPERIMENTAL  PROCEDURE  AND  ANALYSIS 
Materials 

Experimental  materials  include  HY-80  steel  and  7075-T651  aluminum  alloys. 
The  aluminum  alloys  are  solution  treated  at  470°C  for  1/2  hr.,  quenched  In 
water  and  aged  at  120°C  for  different  time  intervals  to  achieve  different 
hardness  levels  After  the  heat  treatment,  the  hardness  values  measured  were 
45,  78,  and  92  on  Rockwell  B  Scale  for  the  three  aluminum  specimens, 
respectively.  Variation  of  the  toughness  parameters  Jjc  and  Ja(2),  as 
functions  of  hardness  is  shown  in  Figure  2.  The  specimen  geometry,  shown  In 
Figure  3,  Is  a  simple  compact  type  specimen  measuring  .150  m  x  .150  m  x  .025  m 
Sharp  45°  face  grooves  are  used  on  the  specimen  occupying  25%  and  80%  of  the 
thickness  tor  the  aluminum  and  steel  specimens,  respectively.  A  fatigue 
crack  is  introduced  in  the  specimens,  following  the  procedures  listed  In  the 
standard  Klc  test^^l  ,  in  order  to  reduce  the  initial  stress  intensity,  Kq. 

Loading  Device 

The  fast  fracture  device  consists  of  reusable,  hardened  steel  arms  which 
are  attached  to  the  specimen  by  means  of  a  key-llke  coupling  and  a  set  of 
bolts  (refer  to  Figure  4).  The  loading  device,  designed  by  Xu^lJ,  also 
consists  of  a  wedge  and  split  pin  arrangement  located  towards  one  end  of  the 
steel  arms.  The  wedge,  the  spilt  pins,  and  the  arms  are  made  of  hardened 
4140  steel  (HRC  50).  Mechanical  restraint  Is  provided  by  a  notched  bolt  made 
of  4340  steel,  quenched  and  tempered  to  a  hardness  of  HRC  38.  The  bolt  In 
this  condition  is  capable  ot  withstanding  about  275  kN  before  fracture. 


Loading  and  Instrumentation 

The  tests  are  conducted  using  a  Tlnlus-Oisen  closed  loop  servohydraullc 
testing  machine.  The  load  applied  on  the  wedge  is  measured  by  means  o£  a 
load  cell  and  is  input  directly  into  an  XY-recorder.  Displacements  are 
measured  by  means  of  clip  gages,  one  attached  to  the  edge  of  steel  arms  (near 
the  wedge)  and  the  other  located  near  the  notched  bolt.  This  enables 
simultaneous  recording  of  the  end  and  the  bolt  displacements.  Alter  lubricating 
the  contact  surfaces  in  the  wedge  and  split  pin  assembly,  the  wedge  is  inserted 
at  a  rate  of  0.05  meters  per  minute.  Loading  is  continued  till  the  bolt 
fractures  and  pre-existing  crack  in  the  specimen  extends  and  gets  arrested. 

The  wedge  is  then  extracted  from  the  device,  which  is  held  down  by  means  ol 
stoppers . 


Analysis  ol  Arrest  Measurements 

At  the  instant  when  the  pin  breaks,  the  elastic  energy  stored  in  the 
arms  is  released  and  the  crack  extends  with  the  aid  of  this  energy.  In  terms 
of  compliance,  this  energy  can  be  expressed  as 


^vend^ 
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Here,  Ua  stands  for  the  energy  stored  in  the  arms  before  the  pin  breaks, 
o 

CD  a  denotes  the  compliance  of  the  device  that  is  pinned  by  the  bolt,  and 
o 

Vend  ls  the  displacement  measured  at  the  load  point.  After  the  pin  breaks 
and  the  crack  extends  from  aQ  to  aa,  the  compliance  of  the  device  increases 
to  a  value,  CUj8  .  The  elastic  energy  stored  drops  down  to  a  value,  Ua^  , 


given  by 


Here,  it  is  assumed  that  the  crack  propagates  under  a  constant  displace¬ 
ment  condition.  The  difference  between  the  two  energies  is  the  total  energy 
available  for  fracture,  and  should  give  an  estimate  of  the  elastic  energy 
release  rate,  Gr  as  follows: 


Gr  “  -  (5) 
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where  B  is  the  effective  specimen  thickness.  Known  Gr,  fracture  resistance, 
can  be  evaluated  using  the  equation 


(Kir)2 


EGr 


(1  -  v2) 


(KiD)2 


(6) 


Finite  Element  Analysis 

A  linear,  elastic-static  finite  element  simulation  of  the  loading  device 
was  performed  using  eight  noded,  quadrilateral,  isoparametric  elements.  Tne 
mesh  shown  In  Figure  5  contains  228  plane  strain  elements,  with  649  nodes  and 
1947  degrees  of  freedom.  The  region  labeled  Zone  1  represents  two  half  Inch 
thick  arms  enclosing  air  gap  an  Inch  tnlck.  The  specimen  sandwiched  between 
these  arms  Is  represented  by  Zone  2.  Since  a  difference  in  elastic  stiffness 
arises  due  to  this  arrangement,  a  weighted  average  of  the  elastic  modulus 
value  was  used  for  each  of  the  two  zones.  A  concentrated  load  was  applied  at 
D  to  simulate  tne  wedge  opening  force  and  a  distributed  load  was  applied  at  C 
to  simulate  the  constraint  due  to  the  notch  bolt.  Points  A  and  B  represent 
the  positions  of  the  two  clip  gages  used  In  the  test  to  measure  the  displacements. 
The  location  of  the  fatigue  crack  before  the  test  is  represented  by  the  point  E. 


After  the  bolt  fractures,  load  applied  at  C  is  removed  and  the  nodes  are  released 
one  by  one  starting  from  E  to  simulate  static  crack  growth. 


Compliance  values  obtained  from  the  finite  element  computations  agreed 


reasonably  well  with  those  from  the  experiment  (refer  to  Figure  6).  For 
constant  displacement,  the  change  in  stress  intensity  with  crack  length  can 
be  calculated  knowing  the  change  In  compliance  values  using  the  equation 
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The  quantity  gives  a  static,  lower-bound  estimate  of  arrest  tougnness. 
Since  the  above  function  In  compliance  Is  very  sensitive  to  small  variations 
In  the  compliance  values,  the  following  polynomial  fit  (with  a  maximum  error 
of  5%)  Is  used,  to  smooth  out  any  scatter: 


cu  a  =  ma' 
U  >a  l 


(8) 


where,  the  value  of  ra  Is  calculated  to  be  3.4596  x  10“^  N~^m~2.  Using 
equations  (7)  and  (8),  the  stress  intensity  change  as  a  function  of  crack 
length  is  calculated  (refer  to  Table  1). 


Fractographlc  Analysis 

Since  the  fracture  appearance  In  the  case  of  fast  fracture  is  significantly 
different  from  that  of  slow  fracture,  an  attempt  Is  made  to  quantify  the 
differences.  A  photogramme  trie  technique [ 22]  ls  modified  to  measure  the 
(x,y,z)  coordinates  of  various  points  on  the  fractograph  with  reference  to  an 
Imaginary  reference  coordinate  systme.  Assuming  that  the  difference  between 
the  adjacent  peaks  and  valleys  on  the  fracture  surface  gives  a  measure  of  the 
surface  roughness,  the  average  difference  between  the  z-coordlnates  is  measured 
using  the  above  technique.  This  difference  between  the  z-coordinates  of 
adjacent  peaks  and  valleys  can  be  taken  as  a  measure  of  surface  roughness 
only  If  the  tilt  axis,  the  Imaginary  x-axls,  and  the  edge  of  the  SEM  specimen 
coincide  with  each  other.  These  measurements  are  done  on  both  slow  and  fast 
fracture  surfaces  and  are  summarized  In  Table  2. 


EXPERIMENTAL  RESULTS 


Figure  7  snows  the  load-displacement  record  for  the  specimens  of  steel 
and  aluminum.  The  displacement  shown  Is  recorded  by  the  clip  gage  attached 
to  the  device  near  the  notched  bar  and  the  load  measured  by  the  load  cell  Is 
the  vertical  load  applied  on  the  wedge.  Displacement  measured  at  the  end  of 
the  device  close  to  tne  wedge  Is  recorded  versus  time  and  Is  shown  In  Figure  8. 

The  load-displacement  values  are  listed  In  Table  3  with  those  obtained 
from  finite  element  calculations  tor  the  difference  specimens  tested.  These 
values  represent  the  load-displacement  situation  at  an  Instant  just  before 
the  fracture  of  the  notched  bolt.  After  the  bolt  fractured,  the  crack 
extensions  measured  are  listed  In  the  last  column  of  Table  3. 

Results  of  Fractographlc  Observations 

Some  of  the  micrographs  of  fracture  surfaces  produced  by  the  fast  fracture 
loading  device  are  shown  here.  For  purposes  of  comparison,  fractographs  of 
the  corresponding  slowly  fractured  surfaces  are  also  shown.  Figures  9  and  10 
depict  the  typical  fast  and  slow  fracture  surfaces,  respectively,  for  steel 
(50X  magnification).  Similarly,  Figures  11  and  12  show  the  typical  fractographs 
of  aluminum  surfaces  (at  the  same  magnification).  On  comparing  the  corresponding 
pairs  of  surfaces,  It  Is  apparent  that  the  appearance  of  fast  fracture  surface 
Is  significantly  coarser  than  that  of  slow  fracture.  Mlcrostructural  features, 
such  as  the  aspect  ratio  (defined  as  the  ratio  of  average  "height"  of  a  dimple 
on  Its  average  width  or  diameter)  seem  to  be  higher  In  the  former  case. 

DISCUSSION 

The  analysis  of  the  arrest  measurements  described  earlier  yields  estimates 
of  Gd  and  KjD  ■  ( EGq/( 1- v^) ] ^ ^ ,  the  average  value  of  the  energy  consumed  by 
the  propagating  crack  and  the  corresponding  toughness.  The  values  of  Gq  and 


Ga  (and  Kq  and  KIa)  tne  statically  evaluated  elastic  energy  release  rates  at 
the  onset  ot  crack  extension  and  arrest  (and  the  corresponding  stress  intensity 
values)  can  also  be  derived  trom  the  compliance  calculations.  These  values  are 
identified  In  Figure  13,  which  also  shows  schematically  the  variation  ot  Gayn 
and  Gstat  with  crack  length. 

At  the  instant  the  bolt  breaks,  elastic  energy  is  conveyed  down  the  arms 
of  the  device  and  test  piece  to  the  pre-existing  crack  tip  at  wave  speeds. 

The  changes  In  the  strain  energy  distribution  are  similar  to  those  produced 
by  a  very  brittle  crack  propagating  from  the  bolt  to  the  test  material.  The 
dynamic  energy  release  rate  Gdyn,  at  the  location  of  the  pre-existing  crack 
tip,  increases  rapidly,  consistent  with  a  very  high  loading  rate  and  rises 
well  above  Gstat  =  Gq,  because  ot  the  kinetic  energy  Inserted  In  the  test 
piece  by  the  fracture  of  the  bolt. 

The  crack  begins  to  extend  when 

Gayn  =  GD(a  =  0)=  JIC(J)  (7) 

where  JIC(J)  signifies  the  J-value  at  the  onset  ot  crack  extension  tor  the 
high  loading  rate.  The  requirement: 

Gdyn  =  GD(a)  (8) 

where  G[)(a)  represents  the  crack  velocity  dependence  of  the  fracture  energy, 
is  satlstled  while  the  crack  Is  propagating  and  determines  Its  speed.  The 
crack  velocity  diminishes  after  the  bulk  ot  the  kinetic  energy  In  the  system 
Is  consumed  as  fracture  energy.  The  crack  arrests  when  a  =  aa  and: 

Gdyn  <  ga  -  GD(a  *  0)  a  Ja(2)  (9) 

It  should  be  noted  that  the  value  of  Gdyn  at  the  onset  ot  crack  extension, 
Equation  (7],  and  at  arrest,  Equation  [8],  are  not  necessarily  the  same  In 
view  ot  differences  In  the  character  of  crack  tip  strain  field  (the  rising 
nature  of  the  J^-curve)  existing  at  these  locations.  The  average  values  of 


value  of  Kjq.  However,  it  should  be  noted  that  even  if  30%  of  the 
released  energy  were  dissipated  in  this  way,  remote  from  the  crack 
tip,  this  would  only  reduce  Kjd  by ~15%. 

(11)  Rate  dependence  of  the  yield  strength.  The  calculations  of  Hoff, 

Hahn  and  Rubinfl^]  snow  that  the  rate  dependence  of  the  yield 
strength  of  steel  can  elevate  the  Jj^-curve  by  an  amount  comparable 
to  Jk;,  independent  of  the  local  resistance  to  fracture. 

(iii)  Rate  dependence  of  the  resistance  to  fracture.  The  changes  in  the 
fracture  appearance  described  in  the  previous  section  reflect  an 
increase  in  the  local  resistance  to  fracture.  This  is  discussed 
more  fully  in  the  next  paragraph. 

As  pointed  out  earlier,  the  fracture  appearance  of  fast  fracture  surface 
is  much  coarser,  and  is  accompanied  by  larger  plastic  deformation,  than  the 
slow  fracture  surface.  Recently,  there  have  been  several  attempts  to  relate 
the  elastic  energy  release  rate  (or  the  corresponding  stress  intensity)  with 
the  micros tructural  features  in  the  case  of  very  tough  and  ductile  materials. 
One  of  the  important  contributions  is  from  Stuwet^]^  whose  work  indicates 
that  fracture  resistance  increases  as  the  magnitude  of  the  microroughness  or 
the  depth  of  an  average  dimple  Increases.  In  more  specific  terms,  doubling 
the  roughness  will  approximately  Increase  the  stress  intensity  by  a  factor  of 
1.3.  Similar  observations  of  Thompson  and  Ashby[24]  also  verify  the  results 
of  Stuwe.  While  these  analyses  are  based  on  empirical  formulations,  it  is 
reasonable  to  state  that  a  fracture  surface  revealing  enormous  plastic  work 
and  large  micros tructural  features  resulting  in  a  coarse  structure,  indicates 
an  Increase  in  the  resistance  to  fracture. 

In  summary,  the  above  effects  add  up  in  such  a  way  as  to  enhance  the 
static  values  of  energy  release  rate  (or  the  corresponding  values  of  stress 


Intensity).  From  the  energy  viewpoint,  the  factors  due  to  the  3  effects  are 
••'1.3,  --<2.0  and  ^1.7,  for  energy  losses,  rate  dependence  of  yield  strength  and 
rate  dependence  of  resistance  to  fracture,  respectively.  These  add  up  to  a 
factor  of~5.  Experimentally  measured  values  (Tables  4  and  5)  Indicate  a 
factor  ~4x  to  ~16x  larger  than  the  static,  rate  Insensitive  values  of  the 
energy  release  rate.  While  the  net  effect  of  the  above  explanations  falls 
somewhat  short,  more  sensitive  instrumentation  and  a  device  of  higher  capacity 
are  required  to  assess  the  dynamic  effects  In  a  thorough  way. 

CONCLUSIONS 

1 •  The  "fast  fracture  device"  Is  able  to  simulate  crack  extension  and  arrest 
In  a  very  tough  material.  However,  the  energy  and  toughness  measuring 
capacity  of  the  present  design  is  marginal  for  tough,  arrester  grade  steels. 

2.  The  relatively  small  crack  extensions  produced  with  relatively  large 
expenditures  of  energy  Indicate  that  the  average  toughness  values  for 
fast  propagating  fractures  are  -^2x  to  ~4x  larger  than  the  values  tor 
slow  crack  extension  In  a  tough  aluminum  alloy  and  steel. 

3.  The  large  values  of  the  toughness  Kjq,  derived  from  the  measurements  may 
have  3  sources  (1)  the  rate  dependence  of  the  resistance  to  plastic  flow, 
(ll)  the  rate  dependence  of  the  resistance  to  fracture  and  (111)  unaccounted 
for  energy  losses  In  tne  device. 

4.  Fracture  surfaces  produced  by  fast  fracture  are  substantially  rougher  than 
those  produced  by  slow  crack  extension.  This  change  Is  consistent  with  an 
Increase  in  the  resistance  to  fracture. 

5.  Rapid  crack  extension  In  a  tough,  ductile  material  can  consume  substantial 
amounts  of  kinetic  energy  released  In  the  process  of  fracture. 


6.  The  Kja(2)  toughness  values  are  intermediate  with  Kid  an<1  ^la*  While  this 
Is  consistent  with  the  view  that  Kja(2)  Is  a  valid  measure  of  the  arrest 
capabilities  of  a  tough  material,  more  direct  evidence  ot  a  close  relation 
between  Kja(2)  and  Kja  Is  called  for  to  support  the  use  of  Kja(2)  in  design. 

7.  Measured  values  of  Gq  Indicate  a  factor  4x  to  16x  the  energy  of  slowly 
extending  fractures  and  the  net  factor  due  to  various  effects  accompanying 
fast  fracture  adds  up  to  ~5x. 


REFERENCES 

1.  Robertson,  T.  S.,  J.  Iron  Steel  Inst.,  1953  Vol.  175,  p.  361. 

2.  Akita,  Y.  and  Ikeda,  K.,  Transporatlonal  Tech.  Rest.  Inst.,  Report 
56,  Tokyo,  December,  1962. 

3.  Yoshlkl ,  M. ,  Kanazawa,  T.,  and  Machida,  S.,  Proc.  Seventh  Japanese  Cong. 
Test  Materials,  1964,  p.  71. 

4.  Crosley,  P.  B. ,  and  Rlpling,  E.  J.,  "Towards  Development  of  a  Standard 
Test  for  Measuring  Kia,"  ASTM  STP  627,  Philadelphia,  1979. 

5.  Hahn,  G.  T.,  Hoagland,  R.  G.,  Rosenfield,  A.  R. ,  and  Barnes,  C.  R. , 

"A  Cooperative  Program  for  Evaluating  Crack  Arrest  Testing  Methods," 

ASTM  STP  711,  Philadelphia,  1980. 

6.  Crosley,  P.  B. ,  and  Rlpling,  E.  J. ,  "Significance  of  Crack  Arrest 
Toughness  (Kia)  Testing,"  ASTM  STP  711,  Philadelphia,  1980. 

7.  Krafft,  J.  M.  and  Irwin,  G.  R.,  "Crack  Velocity  Considerations," 

ASTM  STP  581,  Philadelphia,  1965. 

8.  Crosley,  P.  B. ,  and  Rlpling,  E  J.,  "Crack  Arrest  Studies,"  Second  Annual 
Report  for  EPRI,  Project  RP303-1,  September,  1976. 

9.  Hahn,  G.  T.,  Hoagland,  R.  G.,  and  Rosenfield,  A.  R.,  Met.  Trans., 

Vol.  7A,  1976,  pp.  49-54. 

10.  Hahn,  g.  T.,  Hoagland,  R.  G.,  Rosenfield,  A.  R.,  and  Kanninen,  M.  F., 
"Dynamic  Crack  Propagation,"  G.  C.  Shi,  editor,  Noordhoff  Int.  Publishing, 
Leyden,  The  Netherlands,  1973. 

11.  Kaltholf,  J.  F.,  Belnert,  J.,  Winkler,  S.,  and  Klemm,  W.,  "Experimental 
Analysis  of  Dynamic  Effects  In  Different  Crack  Arrest  Test  Specimens," 

ASTM  STP  711 ,  Phil. ,  1980. 

12.  Hahn,  G.  T.,  Hoff,  R.,  and  Rubin,  C.  A.,  "Strain  Rate  Dependence  of  the 
Deformation  at  the  Tip  of  a  Stationary  Crack,"  Sixteenth  National  Symposium 
on  Fracture  Mech.,  Columbus,  August,  1983. 


■uMtMWJuwwwvrewyvwwwyvv*  I’wwvv*  \s%v%  rewn  yr  v*.  vrv*^)n^j^jrvnjrKX*jrRii7jr*it7Ji'*jrmjrmj<  *jt  ’■^  wj  rjum  ,’JLVLV-V.7\PL,\nmw< 


13.  Joyce,  J.  A.,  "Static  and  Dynamic  J-R  Curve  Testing  of  A533  B  Steel 

Using  the  Key  Curve  Analysis  Technique,"  No.  NUREG/CR-2274 ,  Washington, 
1981. 


14.  Dantam,  V.  K.,  and  Hahn,  G.  T.  “Definition  of  Crack  Arrest  Performance 
of  Tough  Alloys,"  Proc.  of  U.S. -Japan  Coop.  Seminar,  Fracture  Tolerance 
Evaluation,  Honolulu,  December,  1971. 


15. 


Broek,  D. ,  "Elementary  Engineering  Fracture  Mechanics,"  3rd  Rev.  Edition, 
Martlnus  Nljhoff,  Boston,  1982. 


16. 
17  . 


Dantam,  V.  K. ,  M.S.  Thesis,  Vanderbilt  University,  1982, 


Mejias,  H.  D. ,  and  L.  A.  Devldla,  "Influence  of  Side  Grooving  on  Crack 
Arrest  Toughness  of  C-Mn  Steels,”  submitted  for  publication  to  Engr. 
Fracutre  Mechanics,  April,  1985. 


18.  Hlbbltt,  Karlsson  and  Sorensen,  Inc.,  "ABAQUS  -  User's  Manual," 
Providence,  RI ,  1982. 


19. 


Hlbbltt,  Karlsson  and  Sorensen,  Inc.,  "ABAQUS  -  Tneory  Manual,' 
Providence,  RI,  1982. 


20.  ASTM  Standard  E399-78A,  "Standard  Test  Method  for  Plane  Strain  Fracture 
Toughness  of  Metallic  Materials,"  Annual  Book  of  ASTM  Standards,  Part  10, 
Philadelphia,  PA  1979. 


21. 


Xu,  N . ,  "Towards  Development  of  Fast  Fracture  Loading  Device,”  Vanderbilt 
University,  Private  Communication,  1983. 


22. 


Nelson,  J.  L.,  and  Beavers,  J.  A.,  "Application  of  a  Photogramme  trie 
Technique  to  the  Determination  of  the  Orientation  of  Stress  Corrosion 
Fractures,"  Met.  Trans.,  May  1979,  pp.  658-662. 


23. 


Stuwe,  H.  P.  ,  "The  Work  Necessary  to  Form  a  Ductile  Fracture  Surface,' 
Eng.  Frac.  Mech.,  Vol.  13,  1980,  pp.  231-236. 


24. 


Thompson,  A.  W. ,  and  Ashby,  M.  F. ,  "Fracture  Surface  Microroughness,' 
Script  Met.,  Vol.  18,  1984,  pp.  127-130. 


.v 


Ja(2)  (KJ/m2 


FIGURE  3.  Test  specimen  geometry  (drawing  is  not  to  scale) 


o  A1  7075  (HRB  45) 
O  A1  7075  (HRB  78) 
O  A 7  7075  (HRB  92) 
A  Steel  HY-80 


FIGURE  12.  Slow  fracture  surface  of  7075  aluminium 
(HRB45) . 


TABLE  1 


SUMMARY  FOR  FINITE  ELEMENT  ANALYSIS  FOR  STATIC  CRACK  GROWTH 


Aa(m) 

C(mN-l) 

dC/da 

(N-l) 

G(  Jn)-2) 

K(MPa  ml/2) 

0.0 

1.1395x10-7 

1.0664x10-6 

1.2279x105 

164.3 

0.39x10-3 

1.1437x10-7 

1.0690x10-6 

1.2219x10s 

163.9 

0.71x10-3 

1.1471x10-7 

1.0712x10-6 

1.2171x10s 

163.5 

1.42x10-3 

1.1547x10-7 

1.0759x10-6 

1.2064x10s 

162.8 

2.85x10-3 

1.1712x10-7 

1.0855x10-6 

1.1852x10s 

161.4 

5.67x10-3 

1.2010X10"7 

1.1045x10-6 

1.1449x10s 

158.6 

8.48x10-3 

1.2323x10-7 

1.1236x10-6 

1.1062x10s 

155.9 

11.3x10-3 

1.2643x10-7 

1.1430x10-6 

1.0691x10s 

153.3 

14.12x10-3 

1.2968x10-7 

1.1625x10-6 

1.035x10s 

150.7 

16.94x10-3 

1.2399x10-7 

1.1821x10-6 

9.9928xl04 

148.2 

19.76x10*3 

1.3635x10-7 

1.2020x10-6 

9.6664xl04 

145.8 

TABLE 
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TABLE  3 

SUMMARY  OF  LOAD-DISPLACEMENT  DATA 

I 


Material 

Data 

Source 

vend(m) 

vload(m) 

vp  1  n  ( m ) 

Vertical 
Load  (kN) 
on  wedge 

HY80 

Steel 

Finite  Element 
Calculations 

0.0122 

0.0099 

3.81X10-4 

144 

Experiment 

0.0122 

— 

3.74xl0“4 

151 

7075 

Aluminum 

Finite  Element 
Calculations 

0.0124 

0.0101 

3.05x10-4 

158 

Specimen  1 
(HRB  45) 

Experiment  1 

0.0124 

— 

3.03x10-4 

160 

Specimen  2 
(HRB  78) 

Experiment  2 

0.0125 

— 

3.62xl0'4 

161 

Specimen  3  Experiment  3 
(HRB  92) 


0.0125 


3.73x10-4  166 


TABLE  4 


SUMMARY  OF  ARREST  MEASUREMENTS 


Material 

®n/b 

CP,a0 

cu,ax 

U*0 

Ua 

al 

gD 

(mN"1) 

(mN"l) 

(J) 

(J) 

(KJ/m2) 

Steel 

80S 

6.875x10“® 

1.18X10"7 

712.8 

416.0 

5499 

Aluminum 

7075 

HRB  45 

252 

6.32x10"® 

1.43xl0"7 

808.0 

357.0 

2439 

Alumi num 
7075 

HRB  78 

252 

6.32x10"® 

2.29x10-7 

808.0 

223 

609.6 

Aluminum 

7075 

HRB  92* 

252 

6.32x10"® 

>2.6x10-7 

808.0 

<200 

<446.5 

♦Crack  did  not  arrest  at  this  hardness  level  and  therefore  the 
figures  calculated  represent  the  upper  limits. 


TABLE  5 


SUMMARY  OF  ARREST  MEASUREMENTS  AND  CALCULATIONS 


Specimen 

bn/b 

Aa 

mm 

M?a/m 

KJC  . 
MPa^m 

kID 

MPa/m 

*Ja(2) 

MPa™ 

KIa 

MPa^m 

HY-80 

Steel 

80X 

9.0 

402 

214 

1099.3 

259 

385 

7075  A1 
(HRB  45) 

25X 

8.5 

189 

115 

430 

226 

174 

7075  A1 
(HRB  78) 

25X 

50.8 

193 

60 

215 

142 

126 

7075  A1 
(HRB  92) 

25X 

No 

Arrest 

190 

39 

<184 

101 

<104 

